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Chapter  1 

INTRODUCTION 


1.1  Objective 

This  study  is  concerned  with  the  processing  of  radar  signals  using  higher 
order  spectra.  The  purpose  of  radar  signal  processing,  in  general,  is  to 
extract  useful  information  about  target  scattering  mechanisms.  Such  infor¬ 
mation  may  be  used  for  classification  of  unknown  non-cooperative  aircraft 
targets. 

Two  main  factors  can  be  used  to  evaluate  any  radar  signal  processing 
technique:  1)  the  significance  of  the  extracted  features  in  terms  of  providing 
a  better  understanding  of  the  scattering  mechanisms  that  can  lead  to  more 
accurate  description  of  the  target  geometry,  and  2)  the  robustness  of  the 
extracted  features  with  respect  to  noise  contamination,  changes  in  target 
position,  minor  changes  in  the  target  physical  structure  and  so  on. 

The  goal  in  this  study  is  to  define  the  advantages  and  limitations  of 
the  processing  of  radar  signals  using  higher  order  spectra  and  to  evaluate, 
when  possible,  the  utility  of  such  techniques.  Although,  the  final  goal 
is  to  detect  and  classify  unknown  radar  signals,  the  use  of  the  proposed 
processing  technique  in  identifying  scattering  mechanisms  in  the  laboratory 
for  canonical  objects  and  model  targets  is  also  addressed  in  this  study.  The 
focus  is  on  the  interpretation  of  the  responses  that  result  from  higher  order 


1 


spectral  processing  in  terms  of  the  scattering  features  of  radar  targets. 

A  radar  target  recognition  system  that  uses  higher  order  spectral  re¬ 
sponses  as  features  is  shown  in  Figure  1.1.  The  catalog  consists  of  noise- 
free  backscatter  data  for  several  aircraft  at  different  aspect  angles.  The 
“impulse  response”  and  the  “bispectral  response”  of  the  unknown  target 
are  computed  and  then  compared  with  those  of  the  catalog  targets  and  the 
decision  regarding  the  identity  of  the  unknown  target  is  based  on  the  com¬ 
parison.  This  report  is  concerned  with  analyzing  the  relevant  features  and 
performances  of  the  target  recognition  system  shown  in  Figure  1.1.  The 
proposed  radar  signal  processing  technique  is  tested  over  a  wide  spectrum 
of  target  models  and  frequency  bands. 

Higher  order  spectral  processing  of  radar  signals  belongs  to  the  class  of 
time-domain  mechanism  extraction  techniques  [1,2].  Further,  this  research 
effort  can  be  considered  as  a  part  of  the  automatic  target  recognition  prob¬ 
lem.  A  summary  of  the  numerous  contributions  to  the  radar  target  iden¬ 
tification  problem  is  reported  in  [75].  These  contributions  include  studies 
in  the  VHF  and  UHF  frequency  bands  of  a  variety  of  scattering  objects. 
Although,  the  emphasis  in  [75]  is  on  time-domain  analysis  (including  po- 
larimetric  studies),  numerous  efforts  on  feature  selection  and  statistical 
analysis  and  evaluation  in  the  frequency-domain  are  also  summarized.  Key 
results  in  these  investigations  relevant  to  the  problem  at  hand  are  given 
in  [75]. 


1.2  Overview  of  The  Radar  Target  Recogni¬ 
tion  Problem 


Automatic  recognition  of  radar  targets  has  been  extensively  investigated 
during  the  last  three  decades  [24].  In  a  large  number  of  these  investigations, 
and  in  many  current  and  proposed  target  identification  radar  systems,  the 
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Decision 


Figure  1.1:  Block  diagram  of  proposed  target  identification  system. 
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radar  operates  at  stepped  frequencies.  The  recorded  data  are  either  cali¬ 
brated  scattering  coefficients  (complex  valued)  or  normalized  radar  cross- 
section  (RCS)  at  each  of  the  stepped  frequencies.  The  radar  cross  section  of 
a  target  is  proportional  to  the  magnitude  squared  of  the  ratio  of  received  to 
incident  complex  electric  (or  magnetic)  fields  and  is  independent  of  target 
range.  In  [65],  the  radar  cross-section  is  defined  as 

(47r). (power  per  unit  solid  angle  scattered  towards  the  receiver) 
power  per  unit  area  of  the  incident  wave  at  the  target 

(1.1) 

The  recorded  radar  cross-section  data  (amplitude  and  phase),  form  the 
frequency  response  of  a  radar  target.  The  impulse  response  is  the  inverse 
Fourier  transform  of  the  measured  frequency  response  if  an  infinite  number 
of  samples  is  available.  When  only  band-limited  frequency  response  data 
are  available,  which  is  always  the  case  in  practice,  the  corresponding  time- 
domain  representation  is  known  as  the  band  limited  transient  response  of 
a  radar  target  [1,2,3].  Scattering  mechanisms  can  be  determined  from  the 
transient  response  [1,2]. 

The  primary  objective  for  any  radar  target  recognition  system  is  to  re¬ 
liably  and  efficiently  discriminate  between  classes  of  targets  using  features 
extracted  from  the  target  responses  in  either  the  frequency  or  time  do¬ 
main.  Classification  of  radar  signals  can  be  entirely  based  on  the  raw  data 
in  the  frequency  domain.  Numerous  approaches  have  been  followed  for  the 
purpose  of  classification  of  unknown  radar  targets  based  on  the  measured 
band-limited  frequency  response  [26,37,39,68].  Furthermore,  a  limited  num¬ 
ber  of  responses  at  discrete  frequencies  can  be  properly  selected  (feature 
selection  as  in  [26])  in  order  to  reduce  the  dimension  of  the  problem  without 
sacrificing  considerable  loss  of  information  about  the  target. 

Classification  techniques  in  the  frequency  domain  are  divided  into  two 
categories:  parametric  and  nonparametric.  If  prior  information  about  the 
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probability  of  occurence  of  each  target  and  the  likelihood  function  of  the 
measured  noisy  frequency-domain  samples  are  known,  then  it  is  possible 
to  use  an  optimal  parametric  classification  technique  that  minimizes  Bayes 
risk  [38].  In  some  cases,  classification  is  based  on  more  than  one  frequency 
sample,  thus  improving  the  percentage  of  correct  recognition.  The  num¬ 
ber  of  samples  may  be  fixed  (Bayesian  tests  with  fixed  sample  size  [18]) 
or  may  vary  depending  on  the  probability  of  misclassification  (sequential 
probability  ratio  tests  [37]).  If  the  parameters  of  the  underlying  distribu¬ 
tion  of  random  observations  are  unknown,  then  they  can  be  estimated  using 
learning  techniques  [38]. 

An  alternative  to  estimating  the  parameters  of  the  underlying  distribu¬ 
tion  is  to  adopt  a  nonparametric  classification  approach  [38].  Nonparamet- 
ric  classification  is  often  based  on  linear  discriminant  functions  [20]  or  on 
nonlinear  functions  as  in  neural  network  algorithms  [25].  Also,  in  [39],  a 
language  theoretic  approach  to  classification  of  radar  targets  is  developed 
that  uses  syntactic  recognition  algorithms  where  each  target  is  represented 
by  a  finite  set  of  strings.  Relational  graph  matching  for  the  purpose  of 
recognizing  radar  signatures  is  reported  in  [76]. 

Classification  of  radar  signals  can  also  performed  in  the  time  domain 
using  both  parametric  and  nonparametric  decision-theoretic  techniques. 
Time  domain  features  include  transient  response  (using  classical  Fourier 
algorithms),  high  resolution  spectral  processing  (using  autoregressive  mod¬ 
eling  [8]),  and  transient  polarization  responses  (using  both  Fourier  [75]  and 
AR  modeling  [66].  Nonparametric  classification  of  time  domain  features  is 
more  common  simply  because  it  is  difficult  to  parametrize  the  underlying 
distribution,  especially  if  nonlinear  processing  (such  as  autoregressive  mod¬ 
eling)  is  used.  Time-domain  features,  however,  can  be  associated  with  the 
geometry  and  orientation  of  the  target  and  may  provide  valuable  informa- 
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tion  about  the  scattering  mechanisms  of  the  target. 

The  performance  of  any  target  recognition  system  is  dependent  on  the 
type  of  features  used.  The  decision  to  adopt  a  certain  set  of  features  is 
based  on  the  noise  distribution,  knowledge  of  the  target  azimuth  position, 
and  computation. 

I. 3  Organization  of  the  Report 

The  remainder  of  this  report  is  organized  as  follows.  A  background  on  some 
of  the  time-domain  mechanism  extraction  techniques  is  provided  in  Chapter 

II.  Also,  a  review  of  some  of  the  existing  radar  signal  processing  algorithms 
and  feature  extraction  techniques  is  presented  there.  A  summary  of  the 
key  features  of  higher  order  spectra  (HOS),  particularly  the  bispectrum,  is 
also  given  in  Chapter  II. 

Chapter  III  introduces  the  birange  profile  as  an  additional  way  of  dis¬ 
playing  radar  signatures.  Properties  of  the  birange  profile,  estimation  pro¬ 
cedures,  and  interpretations  in  terms  of  scattering  mechanisms  are  dis¬ 
cussed.  Examples  of  birange  profiles  for  hypothetical  targets,  canonical 
objects  and  scaled  aircraft  models  are  presented.  The  validity  of  some  of 
the  assumptions  needed  to  compute  the  bispectrum  is  tested.  Also,  Chap¬ 
ter  III  includes  a  study  on  the  effect  of  unwanted  signals  (such  as  noise)  on 
the  birange  profile.  Finally,  a  characterization  of  radar  targets  in  terms  of 
a  finite  set  of  scattering  centers  and  birange  responses  is  presented. 

The  results  reported  in  Chapter  III  are  based  on  a  classical  bispectral  es¬ 
timation  technique  that  use  a  Fourier  transform  methods.  It  is  known  [60] 
that  classical  signal  processing  of  band-limited  data  using  Fourier  tech¬ 
niques,  especially  in  two  dimensions,  yields  responses  with  limited  resolu¬ 
tion.  In  Chapter  IV,  a  parametric  approach  to  birange  profile  is  investi¬ 
gated.  The  focus  in  Chapter  IV  is  on  estimating  the  birange  profile  via 
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autoregressive  modeling  of  radar  returns.  The  advantages  and  disadvan¬ 
tages  of  this  approach  are  discussed. 

Some  of  the  statistical  properties  of  the  birange  profile  are  discussed 
in  Chapter  V.  Further,  nonparametric  classification  of  radar  targets  based 
on  the  birange  profile  is  presented.  Different  nonparametric  classification 
algorithms  are  compared  in  terms  of  performance  and  efficiency. 
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Chapter  2 

BACKGROUND 


The  research  subject  to  be  investigated  in  the  following  chapters  is  con¬ 
cerned  with  time-domain  representations  of  radar  target  signatures  using 
higher  order  spectral  analysis.  The  scattering  data  (usually  noisy)  cor¬ 
responds  to  calibrated  radar  returns  at  uniformly  spaced  frequencies.  A 
general  description  of  the  radar  system  that  provides  the  raw  data,  and  a 
description  of  the  nature  of  the  measured  data,  are  given  below.  Higher 
order  spectral  processing  of  radar  signatures,  similar  to  Classical  spectral 
analysis  of  radar  data,  produces  signatures  in  the  time  domain  that  are 
related  to  the  geometrical  shape  of  the  target.  A  background  on  spectral 
processing  of  radar  signatures  with  emphasis  on  time-domain  techniques  is 
also  given.  Finally,  an  introduction  to  bispectral  analysis  of  real-valued  sig¬ 
nals,  properties,  and  estimation  techniques  are  summarized  in  this  chapter. 
We  extend  some  of  the  properties  of  the  bispectrum  of  real-valued  signals 
to  complex-valued  signals. 

2.1  Stepped  Frequency  Radar  Measurements 

The  radar  data  considered  in  this  study  (obtained  from  the  Ohio  State 
University  compact  range)  consist  of  coherent  radar  cross  section  measure¬ 
ments  at  a  number  of  uniformly  stepped  frequencies.  These  measurements 
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are  in  units  of  dBm2,  which  is  the  radar  cross  section  of  the  target,  relative 
to  1  square  meter,  and  correspond  to  multipolarization  modes  HH  (hori¬ 
zontally  polarized  transmit  and  horizontally  polarized  receive),  HV,  VH, 
and  VV  [64].  The  data  is  calibrated  so  that  all  system  related  parameters 
are  removed.  A  detailed  description  of  the  measurement  process  and  data 
base  can  be  found  in  [68]. 

The  radar  system  simulated  in  this  study  operates  mostly  in  the  8- 
58  MHz  frequency  range.  This  range  of  radar  frequencies  corresponds  to 
the  case  where  the  size  of  the  target  is  about  1-10  times  the  wavelength 
(resonant  region).  It  was  shown  in  [22]  that  this  range  of  frequencies  con¬ 
tains  sufficient  information  about  the  shape  of  the  target  and  can  be  used 
efficiently  for  target  identification  purposes.  Furthermore,  the  target  signa¬ 
tures  over  this  frequency  range  are  relatively  robust  with  respect  to  changes 
in  aspect  angle.  Other  frequency  ranges  (S-band)  are  also  considered  in  this 
study. 

2.1.1  Scattering  Regions  Suited  for  This  Study 

The  application  of  the  bispectrum  in  radar  signature  analysis  studies  re¬ 
quires  that  the  target  scattering  behavior  includes  significant  multiple  in¬ 
teractions.  This  requirement  imposes  some  constraints  on  the  frequency 
band  to  belong  to  a  specific  scattering  region.  A  brief  description  of  the 
three  major  scattering  regions  as  a  function  of  frequency  is  given  below. 

The  ratio  of  the  target  size  to  the  radar  operating  wavelength  can  be 
used  to  categorize  the  scattering  properties  of  the  target.  The  scattering 
region  targeted  in  this  study  is  mainly  the  resonance  region  where  the 
scattered  returns  include  multiple  interactions. 

The  range  of  frequencies  where  the  wavelength  is  large  compared  to 
the  object  dimensions  is  denoted  by  the  Rayleigh  region.  The  radar  cross- 
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section  in  this  region  is  proportional  to  the  square  of  the  frequency.  The 
range  of  frequencies  where  the  size  of  the  target  is  of  the  order  of  a  wave¬ 
length  is  the  resonance  region. 

The  range  of  frequencies  where  the  wavelength  is  small  compared  to 
target  size  is  the  called  the  high  frequency  or  optical  scattering  region.  The 
coupling  between  currents  in  different  parts  of  the  body  is  insignificant  and 
the  scattering  is  more  localized  in  this  region.  Scattering  in  this  region  is 
mainly  due  to  specular  returns  in  addition  to  diffracted  returns  from  edges 
and  other  surface  discontinuities. 


2.2  Processing  of  Radar  Signals 


Radars  transmit  short  pulses  of  electromagnetic  waves  and  receive  echoes 
from  any  reflector  in  the  path  of  the  transmitted  waves.  The  range  R  to  the 
electromagnetic  waves  reflector  (target)  is  related  to  the  time  A t  elapsed 
between  transmission  and  reception  of  echo. 


(2.1) 


where  c  is  the  speed  of  light  (the  2  accounts  for  the  two  way  path  of  radar 
signals).  The  range  resolution  SR,  which  is  the  minimum  distance  for  re¬ 
solving  two  scatterers,  when  the  transient  response  is  obtained  using  Fourier 
techniques,  is  inversely  proportional  to  bandwidth  B. 


(2.2) 


For  stepped  frequency  radars  with  an  increment  A /  (distance  between  fre¬ 
quency  samples),  an  echo  is  expected  within  a  time  interval  (—T/2,  T/2] 


where  T  is  given  by 


(2.3) 
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which  corresponds,  in  units  of  distance,  to  a  range  r,r  =  cT/2.  Thus,  for 
a  target  of  length  L ,  the  maximum  frequency  increment  A fmax  (maximum 
distance  between  samples  in  frequency  domain)  should  satisfy 

A/™,  <  ^  (2.4) 

This  is  Nyquist  sampling  criterion  to  avoid  aliasing  of  target  responses.  The 
same  criterion  is  sufficient  to  obtain  unaliased  bispectral  responses  [50]. 

The  measured  scattering  coefficients  at  frequencies  w*  to  u>h  are  samples 
of  the  target  transfer  function  (frequency  response)  H(u).  which,  for  a 
single  point  target,  can  be  written  as 

H(u)  =  A(u)  exp  -  wTP)}  (2.5) 

where  u  is  the  radian  frequency  and  TP  is  the  time  required  for  the  radar 
signal  to  propagate  to  the  target  and  back.  The  amplitude  A(u;)  of  the 
returned  signal  may  depend  on  u,  and  similarly  for  the  electric  phase  ip 
(this  is  a  phase  change  introduced  by  target  scattering).  The  range  to  the 
target  is  simply  R  =  It  is  convenient  to  measure  the  range  from  a 
fixed  point  at  the  center  of  the  target.  Then,  R  =  R0  +  r  where  R0  is  the 
distance  to  the  fixed  point  on  the  target.  Similarly,  the  propagation  time 
delay  is  Tr  =  To  +  tT  where  To  =  Then  the  actual  transfer  function  is 

H'(w)  =  exp  (—ju>To)H(w).  (2.6) 

A  target  may  consist  of  P  scattering  points  (or  can  be  modeled  as  such).  If 
there  are  no  interactions  between  various  scattering  points  then  the  transfer 
function  is 

p 

H{v)  =  2  Mu)  exP  ~  uTi)  (2.7) 

i=i 

where  Ai,tpi,Ti  depend  only  on  the  ith  scattering  point.  Multiple  inter¬ 
actions  between  scattering  points  introduce  extra  terms  that  characterize 
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the  coupling  between  the  two  scatterers  and  do  not  correspond  to  a  true 
physical  scattering  object.  Therefore  a  more  accurate  model  of  scattering 
from  a  radar  target  is 

H(u)  =  S^(w)exp{-i(^i(c4>) -wTi)} 

i=i 

+  E  E  4m(w)  exP  {-iCMw)  -  u(th))} 

*:=i  i=i 

+  ....  (2.8) 

Usually  the  magnitude  Am  of  the  scattered  signals  due  to  an  interaction  be¬ 
tween  two  specular  scattering  centers  is  small  compared  to  either  Ak  or  At. 
However,  this  is  not  always  the  case,  especially  for  certain  classes  of  canoni¬ 
cal  targets  [71].  Conventional  radar  targets  such  as  aircraft,  ships,  or  tanks 
are  complicated  to  the  extent  that  one  can  no  longer  exclude  the  possibility 
of  a  significant  number  of  interactions  with  relatively  large  magnitudes.  It 
is  likely,  but  not  necessary,  that  the  interactions  between  scatterers  become 
smaller  in  magnitude  as  the  number  of  interacting  scatterers  is  increased. 
One  focus  of  this  study  is  to  investigate  the  use  of  the  bispectrum  for 
discriminating  between  specular  scattering  and  multiple  interaction  scat¬ 
tering.  The  use  of  the  trispectrum  and  even  higher  order  spectral  analysis 
can  also  be  related  to  more  complex  scattering  mechanisms.  In  this  study 
we  consider  only  the  bispectrum. 

Signal  processing  techniques  can  be  applied  in  either  the  time  or  fre¬ 
quency  domain.  The  measured  frequency  response  has  important  intrinsic 
characteristics  that  permit  the  extraction  of  important  features  such  as 
complex  natural  resonances  (CNR)  [70].  Such  features  are  extracted  by 
applying  spectral  processing  techniques  to  the  recorded  data  over  finite 
uanuvYiuMi.  AiuiiOugu  cue  emphasis  in  this  report  is  on  time  domain  radar 
signal  processing,  a  summary  of  spectral  analysis  studies  when  applied  to 
frequency-domain  data  is  given  below.  It  will  be  shown  later  that  the  same 


spectral  estimation  techniques  can  be  applied  to  either  time  or  frequency 
domain  representations  of  target  responses  and  the  difference  is  only  in 
the  characteristics  of  the  extracted  features.  Therefore,  the  application  of 
higher  order  spectra  is  not  limited  to  time-domain  waveforms  but  may  have 
important  applications  in  the  frequency  domain  as  well. 


2.2.1  Signal  Processing  of  Radar  Signatures:  Frequency- 
Domain 


The  transient  response  properties  of  scattering  objects  are  characterized 
by  their  complex  natural  resonance.  The  number  of  natural  resonances  of 
scatterers  is  infinite.  Kennaugh  [2]  suggested  that  it  is  possible  to  approx¬ 
imately  characterize  the  transient  responses  of  scatterers  at  longer  times 
with  a  finite  number  of  low-frequency  natural  resonances.  The  complex 
natural  resonances  can  be  determined  using  an  iteration  process  which 
searches  for  the  zeros  of  the  impedance  determinant  of  the  object  in  the 
frequency  domain  [72].  Chuang  and  Moffatt  [72]  used  the  Prony  method  to 
extract  complex  natural  resonances  from  approximate  backscattered  ramp 
responses  of  radar  targets. 

If  a  plane  wave  Ei  =  E0exp{j(iot  —  kr)}  is  transmitted  to  illuminate 
a  target,  then  the  backscattered  field  is  E,  =  E0c^f^  exp  {j(iot  —  2 kr)} 
where  H (o>)  is  the  transfer  function  seen  by  the  radar  and  r  is  the  range 
to  the  target.  If  a  ramp  signal  tu(t)  (where  u(t)  =  1  t  >  0,0  otherwise)  is 
transmitted,  then  the  ramp  response  is  related  to  the  frequency  response 
as  follows 


(-»*) 

HP  h  fy  Hi  —/  nnrt  A  »»  i . 

JL  XX^  xtixxip  XV.O|vOllOV>  tvflyvj  t*U>XX  VXllslX  Uv  WXXltVWXj  LlCMXXg  JL'JUXX 


(2.9) 


-  .  .  ...  .  n  «« 
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M*)  =  E 
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cos  r*  +  (2.10) 
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where  is  the  phase  of  H(ju) )  at  u)  =  ivn/T  [72].  If  one  considers  the 
first  few  terms  in  the  above  expansion  (which  is  a  reasonable  assumption 
because  as  the  coefficients  of  the  Fourier  series  are  weighted  by  n2)  then  it 
is  easy  to  notice  the  resemblance  to  the  function 

/(t)  =  53A-exp{7iO  (2.11) 

i=i 

where  the  parameters  .At-,7i  can  be  estimated  using  the  Prony  method  [72]. 
The  poles  7,-  extracted  using  Prony  method  represent  the  dominant  natu¬ 
ral  resonances  of  the  target.  The  complex  natural  resonances  (CNR’s)  of 
a  target  are  aspect  independent  and  therefore  have  been  proposed  as  an 
aspect  insensitive  method  for  discriminating  between  targets.  The  Prony 
algoiithm  is  sensitive  to  noise  contamination  [72]  and  may  produce  inac¬ 
curate  CNR’s.  The  extraction  of  complex  natural  resonances  from  noisy 
radar  signals  remains  an  active  research  area. 

For  electrically  small  objects,  the  CNR’s  can  be  obtained  from  an  in¬ 
tegral  equation  formulation  and  numerical  search  [73].  For  targets  with 
complicated  geometry  the  extraction  of  the  CNR’s  is  from  either  measured 
scattering  data  or  from  measured  surface  currents.  The  Prony  method, 
like  most  high  resolution  techniques,  requires  prior  knowledge  of  the  model 
order  (N  in  equation  (2.11)). 

Frequency  domain  target  identification  studies  focus  on  the  extraction 
of  complex  natural  resonances  that  are  aspect  independent  features.  No 
information,  however,  about  the  scattering  mechanisms  of  a  radar  target 
are  directly  evident  from  the  frequency  response  data. 

2.2.2  The  Impulse  Response 

The  information  about  target  shape,  size,  and  orientation  resides  in  the 
measured  backscattered  frequency  response.  Such  information  can  be  ex- 
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traded  by  examining  the  target  scattering  mechanisms  that  are  displayed 
in  the  target  impulse  response. 

In  general,  electromagnetic  scattering  from  most  objects  tends  to  be 
localized  at  specific  areas.  In  [70],  the  individual  scatterers  are  represented 
as  specular  reflection,  edge  diffraction  and  creeping  wave  mechanisms.  It 
is  often  convenient  to  separately  quantify  the  scattering  cross  section  of 
each  mechanism,  particularly  for  experimental  purposes.  This  can  only 
be  accomplished  if  the  response  to  each  individual  scatterer  is  separable 
and  is  often  done  by  processing  windowed  scattering  spectrum  using  the 
fast  Fourier  transform  (FFT)  [5].  The  FFT  of  the  band-limited  frequency 
response  H (u)  is  the  bandlimited  transient  response.  A  thorough  discussion 
on  the  application  of  transient  and  impulse  response  in  scattering  studies 
can  be  found  in  [2], 

The  relationship  between  the  impulse  response  and  radar  cross  section 
can  be  described  as  follows.  If  an  impulsive  plane  wave  E{  is  reflected  by  a 
linear  scatterer  [70]  then  the  scattered  field  h(t ),  is  defined  as  the  impulse 
response  and  its  Fourier  transform  H(u)  is  the  transfer  function  of  the 
target  as  seen  by  the  radar.  If  a  monochromatic  plane  wave  (exp  {jwi}) 
strikes  an  object,  then  the  phasor  response  E(w)  is  related  to  a  transverse 
component  of  the  scattered  electric  field,  E‘(t  —  *)  at  large  r  by  [70] 

»[j7(u;)exp(jw(i-^))  =^E‘{t--c)  (2.12) 

where  r  is  the  range  from  scatterer  to  radar  and  c  is  the  speed  of  light,  and 
denotes  the  real  part  operation.  The  radar  cross  section  of  the  target  <r 
is  related  to  H{w)  as  [70] 

o-  =  lim  47tr2  ~|  (2.13) 

r-*oo 

=  C7r  \H(io)\2  (2.14) 

A  more  concise  definition  of  the  impulse  response,  [70]  is:  The  im- 
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pulse  response  waveform  is  the  time- dependent  electromagnetic  field  inten¬ 
sity  produced  at  a  fixed  point  in  space  when  a  plane  electromagnetic  im¬ 
pulsive  wave  strikes  an  object .  If  the  impulse  response  is  known,  then  the 
scattered  signal  for  any  incident  pulse  shape  can  be  predicted  [70].  The 
impulse  response  concept  (or  the  time- domain  profile)  is  widely  used  for 
scattering  analysis  [2,5,8,10]  because  it  adds  insight  to  the  complex  rela¬ 
tion  between  the  shape  of  the  target  and  its  reflectivity  [70].  To  further 
emphasize  the  importance  of  the  impulse  response  (or  transient  response) 
in  electromagnetics  research,  it  is  claimed  in  [70]  that  such  techniques  make 
it  possible  to  integrate  all  existing  information  concerning  an  object  into  a 
single  characteristic  function,  thus  combining  the  best  features  of  several 
approximate  methods  such  as  Rayleigh  theory,  physical  optics,  geometrical 
theory  of  diffraction,  etc.  into  a  single  waveform  given  a  specific  aspect  and 
polarization  [70]. 

2.2.3  Signal  Processing  of  Radar  Signatures:  Time- 
Domain 

The  problem  concerned  with  the  extraction  of  scattering  mechanisms  of  a 
radar  target  in  order  to  acquire  some  information  about  the  target  size  and 
geometry  has  been  investigated  by  numerous  researchers  during  the  last 
three  decades.  There  are  two  aspects  to  this  problem: 

•  Resolving  two  closely  spaced  peaks  in  the  impulse  response  or  discrim¬ 
inating  two  distinct  mechanisms  that  are  spaced  by  small  increment 
in  time. 

•  Associating  each  mechanism  in  the  transient  response  with  a  physical 
scattering  phenomenon.  That  is,  resolving  the  integration  of  all  scat¬ 
tering  mechanisms  into  classes  of  mechanisms  (specular,  interactions, 
creeping  waves,  etc...). 
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The  first  aspect  of  the  problem  is  a  signal  processing  problem  where  the 
resolution  of  two  closely  spaced  peaks  is  limited  by  the  bandwidth  of  the 
measured  frequency  response  if  the  transient  response  was  obtained  using 
Fourier  transform  technique.  Some  studies  suggest  using  high  resolution 
spectral  estimation  techniques  in  order  to  overcome  the  bandwidth  lim¬ 
itation  constraint.  The  second  aspect  of  the  problem  is  related  to  the 
electromagnetic  scattering  behaviour  and  depends,  not  only  on  the  band¬ 
width,  but  on  the  range  of  frequencies  used  for  measurements  (Rayleigh, 
resonance,  optical,  etc..). 

Improving  the  resolution  will  generally  have  positive  impact  on  the 
mechanism  identification  problem.  One  should  be  careful,  however,  that 
some  of  the  techniques  used  to  improve  the  resolution  are  prediction  based 
techniques  that  implicitly  extrapolate  the  spectrum  outside  the  measure¬ 
ment  window.  The  extrapolated  spectrum  using  signal  processing  tech¬ 
niques  may  differ  from  the  true  spectrum.  This  could  result  in  false  scatter¬ 
ing  mechanisms  that  cannot  be  associated  with  the  geometry  of  the  target. 
Thus,  a  compromise  between  achieving  a  high  resolution  transient  response 
without  grossly  changing  the  target  scattering  features  is  required. 

This  study  is  concerned  with  both  sides  of  the  time-domain  signal  pro¬ 
cessing  problem.  We  will  use  higher  order  spectra  (HOS)  to  help  in  solving 
the  mechanism  identification  problem  and  we  will  also  study  the  use  of 
parametric  HOS  to  help  solve  the  bandwidth  limitation  problem.  Several 
approaches  to  obtain  high  resolution  transient  responses  are  summarized 
below.  These  techniques,  slightly  modified,  may  be  used  to  obtain  high 
resolution  bispectral  responses. 
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Autoregressive  modeling  of  radar  signatures 

The  key  to  the  estimation  of  high  resolution  transient  responses  is  to  as¬ 
sume  that  the  scattering  data  satisfies  a  certain  model.  Autoregressive 
modeling  of  frequency  response  is  one  possibility.  In  [7]  it  is  assumed  that 
the  measured  scattering  coefficients  yn  at  the  nth  frequency  satisfy  the  au¬ 
toregressive  model 

p 

Vn  =  ]Cai2/n-t  +  en  (2.15) 

»=i 

where  a;  is  the  ith  AR  coefficient  {en}  is  a  zero  mean  Gaussian  distributed 
random  signal  with  variance  p ,  and  p  is  the  model  order.  Autoregressive 
modeling  of  frequency-domain  target  signatures  permits  the  estimation  of 
the  transient  response  scattering  centers  in  terms  of  the  zeros  of  the  above 
difference  equation  satisfied  by  the  data.  The  model  order  is  assumed  to  be 
known  a  priori  and  this  means  (when  translated  to  electromagnetic  terms) 
a  prior  knowledge  of  the  number  of  scattering  mechanisms.  The  autoregres¬ 
sive  model  for  radar  data  is  physically  intuitive  because  responses  at  one 
radar  frequency  can  be  approximately  expressed  as  a  linear  combination  of 
responses  at  other  frequencies.  Autoregressive  modeling  can  be  sensitive  to 
noise  contamination  because  of  implicit  nonlinearities.  In  [7],  this  model 
was  used  for  identification  and  classification  of  commercial  aircraft  tran¬ 
sient  responses.  The  results  in  [7]  indicate  that  the  above  model  produces 
satisfactory  performance  for  SNR  above  10  dB. 

Autoregressive  moving  average  modeling  (ARMA) 

Autoregressive  (AR)  modeling  assumes  that  the  transient  response  of  a 
scattering  object  consists  only  of  peaks  at  different  points  in  time,  and  no 
nulls.  More  accurate  modeling  of  radar  signatures  would  include  both  peaks 
and  nulls  (poles  and  zeros).  Thus,  autoregressive  moving  average  modeling 
is  an  alternative  choice. 


18 


In  [8],  the  measured  scattering  coefficients  at  the  nth  frequency  is  as¬ 
sumed  to  satisfy  the  following  ARMA  model 

y(fn)  =  -  E  aiV{fn-i )  +  S  CAfn-j)  (2-16) 

i=l  j=0 

where  e(*)  is  zero  mean  white  noise  and  p,q  are  the  model  orders.  The 
target  impulse  response  is  given  as 

@  (exp  {j'tt(1  -  2 t/T)}) 

A  (exp  {j*7r(l  —  2t/T)}) 

where  C(*),A(*)  are  polynomials  given  as 

A(®)  =  1  +  aix  +  . . .  +  apxp  (2.18) 

C(x)  =  1  +  Ci®  +  . . .  +  cqx q 

The  orders  p  and  q  may  be  chosen  a  priori ,  or  can  be  estimated  from  the 

data.  One  of  the  findings  in  [8]  is  that  only  20  samples  are  needed  to 

generate  an  impulse  response  that  includes  the  mechanisms  extracted  by 
Fourier  techniques  using  150  samples.  Noise  is  a  problem,  especially  at  low 
SNR’s,  but  this  method  may  be  useful  in  a  low  noise  environment  where 
the  main  objective  is  to  study  scattering  mechanisms.  Although,  both  p 
and  q  have  to  be  specified  (while  in  the  AR  modeling  case,  only  p  is  needed) 
it  is  claimed  in  [8]  that  ARMA  modeling  produces  more  accurate  results 
than  AR  modeling. 

Impulse  response  estimation  via  maximum  entropy  method 

The  maximum  entropy  method  (MEM)  can  be  used  to  obtain  high  reso¬ 
lution  transient  responses  [5,9].  The  MEM  estimation  technique  does  not 
require  the  data  to  be  stationary  or  the  samples  to  be  uniformly  spaced  but 
it  requires  the  exact  values  of  the  autocorrelation  function  of  the  data  [62]. 
According  to  Jaynes  [62],  the  estimation  of  the  autocorrelation  function 
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does  not  yield  an  accurate  transient  response,  especially  if  the  data  is  noisy, 
and  it  does  not  necessarily  maximize  the  entropy  of  the  data.  Unfortu¬ 
nately,  in  the  target  identification  problem,  only  estimates  of  the  autocor¬ 
relation  of  the  frequency  response  are  available. 

In  [5],  the  MEM  transient  response  estimation  technique  is  compared 
to  Fourier  methods.  It  is  shown  that  closely  spaced  scatterers  of  sonic 
canonical  objects  are  easily  resolved  using  maximum  entropy  method  even 
though  the  data  is  nonslationary.  As  in  the  AR  case,  the  model  order  should 
bo  known  a  priori  in  order  to  apply  the  MEM  technique.  Furthermore,  it  is 
important  to  avoid  the  line  splitting  problem  in  MEM  which  indicates  the 
existence  of  two  closely  spaced  scatterers  where  there  is  only  one  scattering 
mechanism. 

Transient,  response  using  MUSIC  algorithm 

The  MUlli  Signal  Classification  algorithm  (MUSIC)  is  used  in  [7]  to  obtain 
high  resolution  estimates  of  scattering  mechanisms.  The  MUSIC  algorithm 
produces  satisfactory  results  when  the  measured  scattering  coefficients  are 
noisy.  The  scattering  data  is  assumed  to  be  sinusoidal  with  a  uniformly 
distributed  random  phase  (the  electric  phase  ifi).  The  covariance  matrix  of 
the  noisy  data  is  given  as 


=  £{/r(u>)J/'*(u;)} 

(2.19) 

=  H{u)H‘{u>)  T  <r2J 

(2.20) 

= 

(2.21) 

where  H'(w)  =  H(u)  -b  n(u;),  <r2  is  the  noise  variance,  and  I  is  the  identity 
matrix.  The  eigenvalues  and  eigenvectors  of  the  covariance  matrix  $  are 
denoted  by  A;  and  e,  respectively.  The  matrix  $  is  Hermitian  and  the  space 
spanned  by  its  eigenvectors  can  be  divided  into  two  subspaces,  one  spanned 
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by  noise  and  the  other  subspace  is  spanned  by  eigenvectors  associated  with 
scattering  mechanisms.  Details  of  the  algorithm  can  be  found  in  [7]. 

This  algorithm  assumes  that  ^  is  uniformly  distributed  over  [0,27r].  It 
is  possible  to  avoid  this  assumption  by  manipulating  the  data  in  H'(io) 
and  .thus  randomizing  ij).  This  is  done  by  dividing  the  measured  data  into 
subbands  with  subband  covariance  matrices.  Then  the  covariance  matrix 
$  is  simply  the  average  of  all  subband  covariance  matrices  [7], 

In  addition  to  high  resolution  estimation  of  scattering  centers,  the  MU¬ 
SIC  algorithm  also  requires  less  data  than  conventional  Fourier  techniques. 
The  main  disadvantage  of  the  MUSIC  algorithm  as  well  as  AR,  ARM  A, 
MEM  techniques,  is  the  loss  of  information  about  the  amplitude  of  scatter¬ 
ing  mechanisms.  These  algorithms  are  mainly  used  to  estimate  the  location 
of  scattering  mechanisms  in  time. 

Extensions  of  the  MUSIC  algorithm  to  higher  order  spectra  are  derived 
in  [30].  Unfortunately,  the  computational  burden  becomes  significant  and 
the  dimensions  of  the  covariance  matrix  increase  rapidly. 

Transient  response  estimation  using  Prony  based  models 

One  may  argue  that  all  the  above  techniques  are  based  on  statistical  spec¬ 
tral  estimation  theory;  however,  in  the  RTI  problem  the  data  may  consist 
of  a  deterministic  signal  corrupted  by  additive  noise.  Furthermore,  the 
stationarity  of  the  radar  scattering  data  is  not  a  strictly  valid  assumption. 
Therefore,  a  more  appropriate  radar  data  model  should  be  based  on  a  deter¬ 
ministic  modeling  approach.  The  Prony  method  is  a  deterministic  approach 
to  modeling  of  radar  data.  In  [67],  the  transient  response  is  estimated  using 
an  improved  version  of  the  Prony  method. 

In  the  Prony  method  the  scattering  yn  at  the  nth  frequency  is  modeled 
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as 

m 

!/„  =  £<«  <2'22) 

i=  1 

where  m  is  the  number  of  scatterers  (known  a  priori).  This  data  satisfies 
the  AR  model 

m 

Vn  f  Yj  aiVn+i  =  0  (2.23) 

»=1 

The  poles  pi  are  estimated  using  the  same  equations  as  AR  parameter 
estimation  methods  where  the  covariance  matrices  are  replaced  with  data 
matrices  and  the  amplitudes  di  are  then  estimated  directly  from  the  data 
in  terms  of  the  estimates  of  the  poles  p;.  Singular  value  decomposition 
of  an  associated  data  matrix  results  in  improved  estimation  of  the  poles 
especially  for  noisy  backscattered  signals.  Although  it  is  necessary  to  know 
the  model  order  a  priori ,  it  seems  that  this  method  is  less  sensitive  than 
other  parametric  techniques  to  changes  in  model  parameters  and  frequency 
band  employed  [67]. 

2.3  Introduction  to  The  Bispectrum 

The  focus  in  this  report  is  on  the  processing  of  radar  data  (treated  as  a  time 
series)  using  higher  order  spectral  analysis  techniques  and  the  bispectrum 
in  particular.  The  application  of  HOS  in  radar  scattering  analysis  is  novel 
and  the  goal  of  this  study  is  to  investigate  the  benefits,  if  any,  of  applying 
HOS  to  target  signatures.  There  is  no  doubt  that  HOS  processing  of  any 
set  of  data  provides  additional  knowledge  about  the  hidden  characteristics 
of  the  data.  Higher  order  spectra,  and  the  bispectrum  in  particular,  are  not 
commonly  used  in  the  radar  community.  For  this  reason,  a  brief  summary 
of  the  bispectrum  in  a  signal  processing  framework  is  given  below.  In  this 
summary,  we  investigate  some  properties  of  the  bispectrum  of  complex¬ 
valued  signals. 
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Signal  processing  using  higher  order  spectra  was  suggested  by  Brillinger 
and  Rosenblatt  in  [12].  The  bispectrum  has  long  been  considered  a  valu¬ 
able  tool  for  analyzing  nonlinear  processes,  and  there  are  numerous  ap¬ 
plications  of  the  bispectrum  in  a  wide  variety  of  nonlinear  problems  and 
disciplines  [11].  But  the  interest  in  the  bispectrum  as  a  general  signal  pro¬ 
cessing  technique  is  recent.  The  reason  is  that  the  bispectrum  is  computa¬ 
tionally  expensive  technique  with  limited  rewards.  That  is,  the  additional 
knowledge  provided  by  the  bispectrum  has  not  justified  the  increased  com¬ 
putational  burden.  Also,  the  mathematics  of  the  bispectrum  are  far  more 
complicated  than  for  other  spectral  analysis  techniques,  and  the  physical 
insight  of  the  bispectrum  is  not  fully  understood.  Today,  with  the  develop¬ 
ment  of  fast  computing  facilities,  more  attention  is  given  to  the  bispectrum 
and  its  applications.  The  main  problem  that  remains  unsolved  is  the  in¬ 
tuitive  interpretation  and  understanding  of  bispectral  processing  of  a  time 
series. 

Currently,  there  are  three  main  applications  for  higher  order  spectra  in 
signal  processing.  These  applications  are  stated  in  [11]  as  follows:  1)  to 
extract  information  about  deviation  from  Gaussianness  (normality)  [11], 
2)  to  estimate  the  phase  and  amplitude  of  a  non-Gaussian  parametric  sig¬ 
nals  [11],  and  3)  to  detect  and  characterize  the  nonlinearities  of  mechanisms 
in  a  time  series  [11].  The  first  application  is  based  on  the  fact  that  for  a 
Gaussian  signal  all  polyspectra  of  order  greater  than  two  are  identically 
zero.  The  second  application  is  based  on  the  fact  that  the  bispectrum  pre¬ 
serves  the  phase  information  of  non-Gaussian  parametric  signal  [11].  The 
third  application  is  due  to  the  fact  that  HOS  plays  a  key  role  in  detect¬ 
ing  and  characterizing  nonlinearities  from  the  output  data  of  a  nonlinear 
system  [11]. 

Recently,  the  interest  in  the  bispectrum  as  a  signal  processing  tool  has 
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moved  into  different  research  areas  and  disciplines.  There  are  numerous 
papers  on  the  bispectrum  in  analyzing  data  structures  such  as  electroen¬ 
cephalogram  (EEG)  data  [30],  seismic  data  [49],  ocean  gravity  waves  [51], 
sonar  data  [11],  underwater  acoustic  transients  [50],  and  bioelectric  data. 
There  are  also  several  applications  of  the  bispectrum  in  texture  discrim¬ 
ination,  reconstruction  of  shift  and  invariant  objects,  tomography,  image 
reconstruction,  non-Gaussian  image  modeling  [59].  Cumulants  and  bispec¬ 
tra  are  also  used  to  separate  stochastic  processes  and  in  training  of  neural 
networks  and  adaptive  antenna  arrays  and  detection  problems  [59].  In  [44] 
Bartelt  et.  al.  propose  an  algorithm  to  estimate  the  phase  and  amplitude 
of  the  original  signal  from  its  triple  correlation  using  bispectral  process¬ 
ing.  Signal  reconstruction  using  the  bispectrum  is  also  investigated  in  [49]. 
Furthermore,  two  sub-optimal  detection  schemes  using  both  spectra  and 
bispectra  are  proposed  in  [56,57]. 

To  our  knowledge,  there  does  not  exist  any  published  literature  on  the 
bispectrum  of  radar  signals  other  than  our  recent  work  in  [21,52].  The 
following  chapters  address  this  novel  problem  of  bispectral  processing  of 
radar  signals  and  its  impact  on  identification  of  scattering  mechanism  and 
target  classification. 


2.3.1  Definitions 


Real-Valued  Random  signals: 


Higher  order  spectra  are  defined  in  terms  of  cumulants  (often  called  cumu- 
lant  spectra)  [11].  Given  a  set  of  n  real  random  variables  {®i,®2,.  •  •  »®u}, 
the  joint  cumulants  of  order  r  =  ki  +  k2  4- . . .  +  kn  are  defined  as 


Ci 


ki...kn 


( /  -\r  ^  in  •  •  •  >^n)\ 

V  3  du )Wi=...=Un 


(2.24) 
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where 


$(w,, . . .  ,wn)  =  E{ex p +  . . .  +  u>n®n)]} 


The  joint  moments  of  order  r  of  {$„}  is  given  as 

“  ((  j)  “&>}■...  flu* 


uii=...=:wn=0 


(2.25) 


(2.26) 

(2.27) 


The  joint  cumulants  are  related  to  the  joint  moments.  For  example,  if  the 
first  moment  is  zero  then  assuming  that  we  have  a  single  random  variable 

*i  (11] 


(2.28) 


C20...0  =  m.2o.  ..0 


C30...0 


m30...0 


C40...0  —  7^40,.  .0  3to20  0 

For  a  zero  mean  process  the  third  order  moment  is  equal  to  the  third 
order  joint  cumulant.  A  random  process  {®i}"=1  is  stationary  to  order 
m  =  ki  +  . . .  +  kn,  if 

If  the  random  process  {*fc}5J=i  is  stationary  with  zero  mean  then  the  second 
and  third  order  joint  cumulant  (identical  to  second  and  third  order  joint 
moment)  are  given  as 


CW..o(fc,  l) 

C2{m) 

Cmo...o(i> 0 
<73(m,n) 


E{x(k)x(l)}  or 
E{x(k)x(k  +  to)} 
E{x(j)x(k)x(l)}  or 
E{x(k)x(k  +  m)X(k  -1-  n)} 


(2.29) 


(2.30) 
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Notice  that  the  fourth  order  cumulant  is  not  identical  to  fourth  order 
moment.  If  the  process  has  non-zero  mean  then  the  cumulant  Cno...o  = 
Covariance{s(A:)»(i)}.  The  general  relation  between  moments  and  cumu- 
'  lants  have  been  derived  by  Leonov  and  Shiryaev  [14]  as 

£{*,...*„}  =  Y,  c(7.)c(ii)...c(7») 

71 

where  the  summation  extends  over  all  partitions  7 ;  of  the  integers  (1, . . .  ,n); 
and  C(.)  are  the  cumulants  of  {®n}.  The  following  examples  (given  in  [14]) 
are  based  on  this  theorem 

E{xk}  =  C(k) 

E{xkxi}  =  G(k)C(l)  +  G(k,l) 

E{XjxkXl}  =  C{j)C{k,l)  +  C{k)C{j^)  +  C{l)C{j,k)  +  C{j,k,l) 

If  we  assume  E{xi}  =  0,V  i,  then  we  obtain  the  relation  [14] 

E{xjXkxixm}  =  C(j,k)C[l,m)  +  C(j,l)G(k,m) 
-\-C(j,m)C(k,l)  4- 

If  xk  is  a  Gaussian  random  variable  then  C{j,  k,l,m)  =  0. 

The  following  are  few  important  properties  of  cumulants  (taken  from 
[13,14]), 

•  Cumulant(o1*i,a2®2  •  •  •  j°n®n)  =  cna2  •  •  •  On  Cumulant(®1, . . .  ,*n). 

•  If  the  random  process  can  be  partitioned  into  two  disjoint  sets  that 
are  independent  from  each  other,  then  Cumulant(si, . . . , En)  =  0. 

•  Cumulant(»i  +yi,...ixn  +  yn)  =Cumulant(xi, . . .,*„)+ 
Cumulant(yi  +  ...,j/n),  where  {x’n}  and  {yn}  are  two  independent 
random  processes.  This  relation  does  not  hold  for  moments. 


26 


The  spectrum  and  bispectrum  are  defined  as  the  Fourier  transform  of  the 
second  and  third  order  joint  cumulants  respectively  [11]. 


OO 


S(w)  = 

Y  ^(rn)  exp  (-j(wm)} 

m=-oo 

(2.31) 

j5(wi,u>2)  = 

oo  00 

Y  Y  C73(jn,n)exp{-i(a;lm  +  a>2n)} 

(2.32) 

m——oo  n— oo 


Unlike  the  spectral  density,  the  bispectrum  is  generally  complex.  A  suffi¬ 
cient  condition  for  the  existence  of  the  bispectrum  is  that  the  third  order 
cumulant  is  absolutely  summable. 

The  bispectrum  is  also  defined  for  a  continuous  time  signal  *(<)  as  [11] 


Cs{ti}T2)  =  E{x(t)x(t  +  rl)®(t  +  r2)}  (2.33) 

/oo  roo 

/  C3(T1,T2)exp{-j(cj1r1+iV2T2)}dT1  dr2(  2.34) 

■OO  «/—  oo 


Little  insight  is  provided  from  the  above  definition  of  the  bispectrum. 
Intuitively,  using  equation  (2.32),  a  spiky  bispectrum  is  the  result  of  flat 
triple  correlation  and  therefore  the  time  series  sc(t)  changes  slowly  over  time 
intervals.  A  flat  bispectrum,  however,  implies  that  the  time  series  changes 
rapidly  over  short  intervals  of  time. 

An  alternative  definition  of  the  bispectrum  is  based  on  the  Cramer  (or 
Fourier-Stieltjes)  spectral  representation  dZ{y) )  [11],  where  the  time  series 
{»(&)}  can  be  written  as 

x(k)  =  7T  /  exp  {jwfc}  dZ(u>)  (2.35) 

where  dZ(u>)  is  random  process  with 


E{dZ{u>)}  =  0  (2.36) 

E{dZ(u>i)dZ~(w2)}  =  2ttS(w)cLi;,  =  w2  =  w 
2?{dZ(u>i)d,Z(w2)dZ*(w3)}  :=  jB(wi,u>2)dwidw2,  wj  +  w2 -- a>^(2.37) 
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then 


uu 

®(ft)  =  £  -^(w)exp  {jwA:}  then  (2.38) 

u>=  —  oo 

S(u>)  =  (X(w)X»)  and 
j5(u>i,c»;2)  =  (JC(wi)jr(a;2)^C*(wi  +  W2)).  (2.39) 

where  ()  denotes  the  ensemble  average.  The  above  definition  of  the  bis¬ 
pectrum  is  more  insightful  than  the  cumulants  based  definition.  This  def¬ 
inition  implies  that  the  bispectrum  at  a  particular  frequency  pair  (u>;,u>j) 
(or  bifrequency)  is  nonzero  if  at  least  one  of  the  three  spectral  responses 
Ar(u>j),  and  (w,-  +  uij)  is  correlated  with  the  other  two. 

Real- Valued  Deterministic  signals 

The  definition  of  the  bispectrum  given  in  the  previous  section  is  based  on 
the  expectation  of  the  triple  product  of  a  random  signal.  In  many  practi¬ 
cal  deconvolution  and  filtering  problems  the  signal  is  deterministic.  Also, 
the  addition  of  noise  to  the  deterministic  signal  produces  a  nonstationary 
random  signal.  The  expectation  based  definition  of  the  bispectrum  is  not 
applicable  for  time  varying  deterministic  signals  even  if  these  signals  are 
embedded  in  noise.  The  following  definition  of  the  bispectrum  is  applicable 
if  the  measured  time  series  has  a  deterministic  time  varying  component 
added  to  a  random  component  such  as  noise.  The  distinction  between  de¬ 
terministic  and  time  varying  signals,  however,  is  only  in  the  statement  of 
the  problem  and  not  in  the  estimation  procedure. 

Bispectrum  of  Real-Valued  Energy  signals  [79]: 

If  {x(k)}  is  a  real  energy  signal  i.e., 

CO 

£*=  £  |*(fc)|2  <  oo  (2.40) 

fc=- 00 
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(2.41) 


then,  the  spectrum  of  {s(fc)}  is  given  as 

X(u>)  =  x(k)exp{—ju)k} 

k—- oo 

and  the  second  and  third  moments  are  given  as 

M^(m)  =  ^2  x(k)x(k  +  m) 

k=—oo 

OO 

M3  (m,n)  =  yt  x(k)x(k  +  m)x(k  +  n) 

-OO 

the  spectral  and  bispectral  densities  are  respectively 


Sx{io)  =  X{u)X'(u>) 

, ^2)  =  X (wj)A  (ti^A"  (uq  +  W2) 


(2.42) 


(2.43) 


Notice  that  the  bispectrum  for  deterministic  signals  is  defined  in  terms  of 
moments  and  that  of  random  signals  is  defined  in  terms  of  cumulants.  Also 
notice  that  the  estimation  procedures  for  both  third  order  cumulants  and 
third  order  moments  are  identical. 

Bispectrum  of  Real- Valued  Periodic  signals  [79]: 

If  the  sequence  {*(&)}  is  periodic,  i.e.,  x(k  +  N)  =  x(k)  then  x(k)  is  a  power 
signal  with  finite  average  power  given  by  Sx  =  |x(A;)|2.  The  line 

spectrum  of  {*(&)}  is  obtained  using  discrete  Fourier  transform  defined  as 


x(k ) 
X(n) 


4  E  ^ (n)  exP  {j2irfcn/ JV} 

n=0 

N-l 

y  x(k)  exp  {— j2irkn/N} 

k—0 


(2.44) 


where  X(n)  is  also  periodic  with  period  N.  The  second  and  third  order 
moments  are  given  as 


R*(m)  ~ 
R^(m,n)  = 


N 

N 


N-l 

\  *!*(  l*  I  w) 

/  .  I  III/) 

k—0 


N-l 

y  x(k)x(k  -f  m)x(k  4-  n) 
k=  0 
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and  the  spectral  density  and  bispectrum  are,  respectively 

S-(P)  =  jX(p)X'(p)  (2.46) 

B’(p,q)  =  ±X(p)X(q)X-(p  +  q) 

The  estimation  of  polyspectra  of  periodic  signals  is  similar  to  that  of  energy 
signals  with  finite  data  lengths. 

2.3.2  Properties  of  The  Bispectrum  of  Real- Valued 
Signals 

The  cumulant  C(m,n)  of  a  real  stationary  random  sequence  {*(&)}  has  six 
symmetry  regions  in  the  two-dimensional  time  domain.  These  symmetry 
regions  are  defined  by  the  following  equations  [11] 

C7(m,n)  =  C(n,m)  (2.47) 

=  C(—n,m  —  n) 

—  C(n  —  m,  — m) 

=  C(m  —  n,  —n) 

=  C(—m,n  —  m) 

Therefore  it  suffices  to  compute  third  order  cumulant  in  one  region,  and 
the  other  sectors  can  be  generated  using  the  above  relations.  Figure  2.1 
shows  how  these  symmetry  regions  are  distributed  in  the  (m,  n)  plane. 
Notice  that  with  this  type  of  skewed  symmetry,  the  correspondence  between 
data  points  in  one  sector  is  not  a  twofold  symmetry  in  the  other  sector. 
Figure  2.2  shows  how  the  data  is  arranged  in  different  sectors  of  the  (m,n) 
plane.  For  example,  samples  that  are  distributed  horizontally  in  the  sector 
Tn ,  7x  ^  0,  77i  ^  7i  arc  the  ITU rror  image  of  the  sector  m ,  ti  ^  0,  tti  n  tl 
but  they  are  arranged  vertically  in  the  sector  m  <  0,n  >  0.  The  arrows  in 
Figure  2.1  indicate  the  arrangement  of  cumulant  samples  (from  c  to  b)  in 
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each  symmetry  sector  of  the  third  order  cumulant  plane.  All  sides  of  the 
symmetry  triangles  are  labeled  accordingly  for  all  symmetry  regions  (see 
Figure  2.1).  Furthermore,  the  data,  samples  on  the  axis  m  >  0  has  a  mirror 
image  on  the  axis  n  >  0  but  the  third  image  is  on  the  diagonal  defined  by 
m,n  <  0,  and  m  —  n.  Therefore  a  set  of  responses  on  the  axes  m>  0,  or 
n  >  0  may  look  similar  but  less  dense  on  the  third  axis  m  =  n,m,n  <  0. 

The  bispectrum  i?(uq,u;2)  of  a  real  stationary  sequence  x(k)  has  twelve 
symmetry  regions.  The  twelve  sectors  in  the  two-dimensional  frequency 
domain  are  related  as  [11] 

B{wi,u2)  =  i?(w2,uq)  (2.48) 

=  -B*(-W2,-Uq) 

=  B*(-u>  i,-w2) 

=  B(-uft  -wa,u>a) 

=  J3(uq,-uq  -  u>2) 

=  B(u> 2,  -uq  -  0>2) 


The  bispectrum  is  generally  complex  with  an  amplitude  and  phase.  Also, 
for  discrete  signal  processing,  the  bispectrum  is  doubly  periodic  with  period 
27T. 


B(o)i,io  2)  =  B(u)i  +  2fc7r,o;2  +  2i7r).  (2.49) 


Figure  2.3  shows  the  twelve  symmetry  regions  of  the  bispectrum.  The 
principle  triangle  in  the  first  quadrant  of  the  (aq,u>2)  plane  is  defined  by 
uq  >  u)2  and  uq  <  7 r.  There  is  an  additional  triangle  that  is  adjacent  to  the 
principle  triangle  and  is  defined  by  7r/2  <  uq,u>2  <  2/37T.  The  two  triangles 
combined  define  the  first  of  the  twelve  sectors  of  the  bispectrum.  As  in  the 
cumulant  case,  the  knowledge  of  the  bispectrum  in  one  sector  is  sufficient 
to  determine  the  other  sectors  using  the  above  relations.  Also,  beca.use 
of  the  skewness  of  symmetry  regions,  one  should  be  careful  when  trying 
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to  read  the  bispectral  responses.  Figure  2.4  shows  a  data  arrangement 
in  the  bispectral  plane  where  the  solid  lines  represent  the  boundaries  of 
the  symmetry  regions.  Notice  that  the  data  is  arranged  parallel  to  these 
boundaries  as  indicated  by  the  arrows  (from  side  c  to  side  b )  in  each  of 
the  six  major  symmetry  regions.  Each  of  these  six  symmetry  regions  has 
twofold  symmetry  along  an  axis  S  where  S  €  {u> 2  =  o>i,o>2  =  — Wi,u>2  = 
— 2u>j,u>2  =  — |wj}  (see  Figure  2.4). 

Due  to  the  hexagonal  shape  of  the  bis^ectrum,  it  is  a  bit  difficult  to 
follow  its  periodicity  in  the  u;1,u;2  plane  every  increment  of  27r  in  each 
frequency  axis.  Figure  2.5  shows  how  the  (u>i,u>2)  plane  is  subdivided  into 
hexagons  of  bispectra  with  period  27T. 

Aliasing  in  the  bispectral  domain  occurs  in  the  triangle  abf  (see  Figure 
2.5)  adjacent  to  the  principle  triangle  oaf  [50].  Notice  that  if  a  rectangular 
display  of  bispectral  signatures  is  used  then  the  responses  in  the  triangles 
{abc,  egb,  bgd ,  bde,  bef ,  aft/}  (see  Figure  2.5)  are  all  symmetric  with  different 
data  arrangements.  The  reason  for  this  symmetry  is  that  the  triangles 
{abc,  abf}  belong  to  the  first  bispectral  hexagon  centered  at  (0,0),  the 
triangles  {egb,  bdg}  belong  to  the  hexagon  centered  at  (27T,  0),  and  the 
triangles  {bde,  bef}  belong  to  the  hexagon  centered  at  (0,27r).  The  same 
symmetry  is  repeated  in  the  third  quadrant  of  the  rectangular  display  of 
bispectral  signatures.  Estimates  in  these  regions  are  unreliable  because 
they  may  be  due  to  an  overlap  of  bispectral  responses  in  the  three  adjacent 
bispectral  hexagons. 

An  important  property  of  the  bispectrum  is  that  it  vanishes  when  {*(&)} 
is  a  zero  mean  stationary  Gaussian  signal.  Also,  if  the  waveform  x(k)  has  a 
bispectrum  jB(u>i,u;2),  then  the  shifted  signal  x(k  — n )  has  the  same  bispec¬ 
trum.  Thus,  the  bispectrum  as  well  as  the  spectral  density  are  unaffected 
by  a  linear  phase  shift  of  the  spectrum  of  (*(n)}.  Further,  if  e(k)  is  white 


32 


non-Gaussian  noise  defined  such  that  E{e(k)}  =  0,  E{e(k)e(k  +  m)}  =  cr2 
and  E{e(k)e(k  +  m)e(k  +  n)}  =  f35(m,n).  The  bispectrum  of  e(k)  is  con¬ 
stant  over  the  (101,11)2)  plane  with  magnitude  /3. 

In  many  practical  situations,  the  waveform  {x(fc)}  is  bandlimited  where 
_X”(u;)  is  defined  over  the  frequency  interval  [<ol^h\-  The  bispectrum  of 
x(k),  Bx((Oi,io2)  is  then  nonzero  only  in  the  triangle  <  u>i,u;2  <  wh,  Wi  + 
u2  <  o)H  [49].  Figure  2.7,  shows  the  bispectrum  of  a  bandlimited  signal. 

2.3.3  The  Bispectrum  of  Complex- Valued  Signals: 

The  above  definitions  and  properties  of  the  bispectrum  are  limited  to  real¬ 
valued  signals.  Little  information  can  be  found  on  extensions  of  these  def¬ 
initions  to  complex-valued  signals  [76,43,44].  However,  for  the  radar  signal 
processing  problem,  the  bispectrum  of  complex  signals  is  needed.  Below 
we  discuss  definitions  and  properties  of  the  bispectrum  of  complex- valued 
signals. 

If  the  data  {x(k)}  is  a  complex  stationary  process,  then  the  second 
order  cumulant  is  defined  as  C2(m)  =  E{x’k(k)x(k  -f  m)}.  It  is  also  pos¬ 
sible  to  alternate  the  conjugate  so  that  C2(m)  =  E{x(k)x“(k  +  m)}.  For 
the  third  joint  cumulant  case,  the  conjugate  can  be  placed  either  on  one 
or  two  entries  of  the  triple  product  in  (2.30).  Only  one  In  fact,  plac¬ 
ing  the  conjugate  in  a  particular  position  defines  which  one  of  the  six 
symmetry  relations  given  in  (2.48)  remains  valid.  In  general,  the  bispec¬ 
trum  of  complex  signals  has  twofold  symmetry  about  an  axis  S  such  that 
S  £  {wj  =  u>2,  (o2  =  —\<Oi,  (Oi  =  -\w2,  <Oi  =  —  u;2}.  Figure  2.6  shows 
the  configuration  of  the  bispectral  response  with  respect  to  its  symmetry 
axis  for  all  possible  positions  of  the  complex  conjugate  in  the  cumulant 
(as  designated  in  each  symmetry  region  in  Figure  2.6).  If  the  first  quad¬ 
rant  u>x,  (o2  >  0  is  preferred  then  the  third  order  cumulant  is  defined  as 
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Figure  2.1:  The  six  symmetry  regions  of  the  third  order  cumulant  function 
c(m,n) 
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Figure  2,2:  Data  arrangement  of  third  order  cumulants. 
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Figure  2.3:  The  twelve  symmetry  regions  of  the  bispectrum  B(ui, u2) 
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Figure  2.5:  Periodicity  of  the  bispectrum  in  the  (u>i, io2)  plane. 
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W1 


Figure  2.6:  Distribution  of  the  responses  in  the  bispectral  domain  depen¬ 
dent  on  the  placement  of  the  conjugates  in  the  third  order  cumulants 
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C3(m,n )  =  E{x*(k)x(k  +  m)s(fc  +  n)}. 

The  following  relations  define  the  bispectrum  of  complex-valued  signals 
for  all  possible  positions  of  the  complex  conjugate. 

E{x’"(k)x(k  +  m)x(k  +  n)}  *-£-*  (X(u>1)X(w2).X’“(u>i  +  u>2)) 
with  the  symmetry  relation  B(u> i,w2)  =  l?(u>2,u>j) 

E{x(k)x*(k +m)x(k +n)}  (Jf’(-o>i)X(a;2)X(-a>i  —  u>2)) 

with  the  symmetry  relation  B( u>i,u>2)  =  .B(u>1}  — a>x  —  a>2) 

E{x(k)x(k  4-  m)x*(k  +  n)}  {X(o>i)^l'“(—  w2)X(—  u>i  —  u>2)) 

with  the  symmetry  relation  5(o>1}w2)  =  f?(— u>i  —  u>2,u>2) 

+  m)x“(k  +  n)}  <-^— >  (X‘(— Wi)X*(— w2)A’(— —  u>2)) 
with  the  symmetry  relation  .6(0;!, u>2)  = 

E{®*(fc)®*(jfc+m)®(Jb  +  n)}  (X*(-o>1)X(o;2)X*(a;; +a;2)) 

with  the  symmetry  relation  B( u>i,u>2)  =  B(-u)2 •  —  u>i,w2) 

E{**(Jfc)*(Jfc  +  m)x*(k  +  n)}  (X(Wl)X*(-a;2)Jr(a>a  +  w2)) 

with  the  symmetry  relation  B(w i,u>2)  =  — wx  —  u>2) 

The  symmetry  properties  for  each  position  of  the  complex  conjugate  can 
be  easily  derived  from  the  above  equations.  Note  that  it  is  possible  to  ob¬ 
tain  one  bispectral  response  from  other  bispectral  responses  using  cumulant 
transformations.  For  example  if  Ri(m,n )  =  E{x*(k)x(k  +  m)x(k  +  n)}  and 
•R2(m,n)  =  E{x(k)x*(k  +  m)x(k  +  n)}  then  i?i(m,n)  =  R2(— m,n  —  m). 
Similarly,  if  i?3(m,n)  =  E{x(k)x(k+m)x*(k+n)}  then  Ri(m,n )  =  R3(m  — 
n,—n ),  and  Jij(m,n)  =  R\(—m,n  —  m)  where  i24(m,n)  =  E{x(k)x*(k  + 
m)x*(k  +  n)}.  Other  relations  can  be  obtained  similarly. 
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Figure  2.7:  Bispectrum  of  bandlimited  signals 


2.3.4  Examples  of  Interpretation  of  The  Bispectrum 

Some  of  the  properties  of  the  bispectrum  are  better  understood  in  the 
light  of  the  following  examples.  For  this  reason,  the  following  examples 
are  common  in  the  bispectrum  literature  [11,48],  and  are  mentioned  here 
because  of  their  relevance  to  bispectral  processing  of  radar  signals. 

Example  1  [11]:  Quadratic  phase  coupling  of  two  harmonic  components 
of  a  random  process  due  to  quadratic  nonlinearities.  The  signal 

x(k)  =  cos(0j k  -b  (f>i)  +  cos(02k  +  <f>2)  +  cos((#i  +  02 )k  +  <f>i  +  <f>2)  (2.50) 

where  <j>i  and  <^2  are  random  variables,  uniformly  distributed  over  [0, 27r]. 
The  bispectrum  of  x(k )  is  a  single  impulse  at  (0i,02), 

Bx{u i,w2)  =  -  0i,u>2  -  02)  (2.51) 

This  example  demonstrates  that  the  bispectrum  is  useful  in  detecting  quadratic 
phase  coupling  between  harmonic  components  [11]. 

Example  2  [48]:  If  {aj(fc)}  is  a  zero  mean  non-Gaussian  stationary  random 
process  which  is  the  input  into  a  Unear  filter  (with  impulse  response  h{k)) 
then  the  bispectrum  of  the  output  signal  y(k)  is  given  as 

By(u>i,u  2)  =  Bx(u>i,w2)H(u)i)H(u)2)H*(u)i  +  u>2)  (2.52) 

The  special  case,  where  x{k)  is  white  non-Gaussian  (j9*(u;i,u>2)  =  1),  is 
useful  for  parametric  estimation  of  the  bispectrum  [42],  where  the  bispec¬ 
trum  of  the  filter  output  is  simply  H(wi)H(w2)Hm(w j  +  u>2). 

Example  3  [11]:  Let 

y(k)  =  x(k)  +  ax2{k )  (2.53) 

where  x(k)  is  a  zero  mean  stationary  Gaussian  random  process.  Then 

By(u  i,uj2)  =  2a{5I(w1)5E(u>2)  +  Sx(w2)Sx(ui,v2)}  (2.54) 
+ar(0){Sx(w1)£(u;2)  +  S'E(w2)£(u;1) 

-\-Sx{w2)8{y)\,  w2)}  +  o(o3)  (2.55) 
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where  Sz(u>)  is  the  power  spectral  density.  This  example  shows  the  role 
of  the  bispectrum  in  analyzing  nonlinear  systems.  Exploring  the  relation 
between  the  bispectrum  and  the  spectrum  of  y(k )  may  indicate  the  type  of 
nonlinearity  applied  to  the  process  x(fe),  particularly  for  a  <  1. 

The  above  examples  are  basic  examples  of  the  major  applications  of 
the  bispectrum  in  signal  processing.  Additional  examples  can  be  found 
in  [11,12,42,48].  An  extensive  study  on  the  analysis  of  bilinear  systems 
using  bispectra  can  be  found  in  [14,48]. 

2.3.5  Bispectral  Estimation  Techniques 

The  problem  addressed  in  this  section  is  the  estimation  of  the  bispectrum 
of  the  finite  complex- valued  data  sequence  ®(1),®(2), . . .  ,  ®(1V).  As  in  the 
spectral  estimation  case,  the  notion  of  expectation  if  x(k)  is  random  and 
the  notion  of  infinite  summation  if  x(k)  is  deterministic  can  be  replaced 
by  finite  summation  over  the  available  data  record  weighted  by  the  factor 
1/N.  The  expectation  of  the  triple  product  can  be  replaced  by  segmenting 
the  data  records  and  then  averaging  the  triple  correlations  of  all  segments. 
Segmentation  also  minimizes  the  estimate  variance,  thus  improving  the 
stability  of  the  bispectral  estimates. 

The  following  techniques  produce  approximate  estimates  of  the  bispec¬ 
trum  and  the  true  bispectral  response  is  obtained  only  if  the  data  record 
is  of  infinite  length  ( N  — »  oo).  The  estimation  techniques  given  below 
are  classical  Fourier  based  bispectral  estimation  procedures  and  can  be  im¬ 
plemented  using  the  fast  Fourier  transform.  There  are  several  parametric 
estimation  techniques  that  are  based  on  modeling  the  data  as  an  AR  or 
ARM  A  process  [11,14,27,28,32]. 


43 


Indirect  method 


The  indirect  bispectral  estimation  method  [11]  can  be  outlined  as  follows, 
(this  method  is  given  in  [11]) 

1)  Segment  the  data  into  K  records  of  M  samples  each. 

2)  Subtract  the  average  value  of  each  record  (to  get  a  zero  mean  record). 

3)  Estimate  the  third  moment  for  each  segment  rO)(m,n),  i  =  1, . . . ,  k  by: 

r«  =  £  *(,T  (l)xW(l  +  m)®W(*  +  n)  (2.56) 

M  l=n 

where  Si  =  max  (0,  -m,  -n),  s2  =  min  {M  -  l,M-l-m,M-l-  n)}, 
and  At  is  the  sampling  interval. 

4)  Average  r^(m,n)  over  all  segments 

1  * 

JR(m,  n)  =  -  r(,)(m,  n)  (2.57) 

k  »= i 

5)  The  bispectrum  is  computed  using  two-dimensional  discrete  Fourier 
transform  (2DFT)  of  i?(m.,n). 

m=L  n—L 

5(wi,u>2)  =  £  R(m’n)W{rn,n)exp{-j(wirn  +  w2n)}  (2.58) 

m=—L n=—L 

where  L  <  M  —  1  is  the  triple  correlation  lag  and  W(m,n )  is  a  two- 
dimensional  window  function  that  can  be  generated  from  a  one  dimensional 
window  d(m)  as  follows 


W{m,n)  =  d(m)d(n)d(n  —  m) 


(2.59) 


where  d(m)  =  d(—  m),  d(0)  =  1,  and  H(w)  >  0  Vu;.  There  are  several 
choices  of  d(m),  [11,14].  One  popular  choice  is  the  Optimum  Bias  Supre- 
mum  window  [11]  defined  as 


7r  m 


sm- 


+ 


(2.60) 
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Figure  2.8  shows  the  two-dimensional  optimum  bias  supremum  window 
function  in  frequency  and  time  domains.  An  impulse  in  the  bispectrum 
before  windowing  is  convolved  with  the  frequency  response  of  the  window, 
thus  producing  an  elliptically  spread  response.  The  window  design  criterion 
is  similar  to  the  one  dimensional  case  [11]  but  has  to  satisfy  the  symmetry 
properties  of  third  order  cumulants. 

The  above  method  is  similar  to  the  Correlogram  spectral  estimation 
method  with  the  slight  difference  that  we  do  not  have  to  worry  about 
nonnegativity  of  the  estimate.  Van  Ness  [35],  has  shown  that  if  the  12th 
moment  of  the  sequence  {®(fc)},  JS{|a;(fc)|12}  is  finite,  and  all  cumulants 
of  order  less  than  or  equal  to  six  are  also  defined,  then  the  estimates  of 
the  bispectrum  using  the  indirect  method  are  asymptotically  ( N  — >  oo) 
complex  normal. 

Direct  method 

The  direct  method  [11]  is  based  on  the  definition  of  the  bispectrum  as 
the  ensemble  average  of  the  product  -V(u>1)A’(a>2)A'*(a;1  +  w2),  and  does 
not  require  using  the  two-dimensional  FFT.  Further,  bispectral  estimation 
as  given  below  is  similar  to  the  periodogram  spectral  estimation  method. 
The  direct  method  consists  of  the  following  steps  (this  method  is  outlined 
in  [11]): 

1)  Segment  the  data  into  K  records  of  M  samples  each. 

2)  Subtract  the  average  of  each  record. 

3)  Generate  the  DFT  of  all  K  records 

i  M- 1 

XW  =  ~M  £  x(k)exp{-j2nk\/M)}  (2.61) 

**—  k=0 
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Figure  2.8:  The  two-dimensional  window  (2.59),  (2.60)  in  the  time  (top) 
and  the  frequency  domain  (bottom) 
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4)  Let  M  =  MyxN  and  My  =  2Ly+l.  Also,  assume  that  the  time  increment 
between  samples  is  At.  Estimate  the  bispectrum  of  each  segment  using 

&(i)(Ai,A2)  =  (NAt)2  £  E  Y^iXy  +  ky)  (2.62) 

ki=-Li  ki=~Li 

F^(A2  +  k2)Y^  (Aj  -f  A2  +  &i  4-  k2) 


5)  The  bispectrum  estimate  of  x(k)  is  then  given  as 

£“(wi,a>2)  =  ££4(lW"0  (2.63) 

i= 1 


where  coy  =  Ai  and  w2 


Both  direct  and  indirect  bispectral  estimation  methods  have  high  vari¬ 
ance  and  therefore  need  a  large  number  of  data  samples.  The  estimates 
can  be  made  smoother  if  more  segments  are  used  but  this  will  decrease  the 
resolution  of  the  bispectral  estimates  as  a  result  of  introduced  nonstation- 
arity  problems  [11].  In  general,  however,  the  above  methods  are  easy  to 
implement  using  FFT  algorithms. 


Other  Methods 

There  a,re  other  more  elaborate  methods  for  estimating  the  bispectrum. 
Hinich  [48]  proposed  another  periodigram  bispectral  estimation  technique. 
His  technique  is  similar  to  the  direct  method  with  slight  modification. 
Hinich’s  estimation  method  is  generally  useful  for  avoiding  aliasing  in  the 
triangle  adjacent  to  the  principle  triangle  although  oversampling  is  another 
way  to  avoid  the  aliasing  problem.  Statistical  properties  of  the  Hinich  al¬ 
gorithm  can  be  found  in  [36,50]. 

Thompson,  [63],  used  a  multiwindow  bispectral  estimation  method  where 
the  information  in  a  narrow  band  is  summarized  by  the  coefficients  of  its 
expansion  in  Slepian  (discrete  prolate  spheroidal)  sequences.  The  bispec¬ 
trum  estimates  are  then  computed  by  averaging  over  time  or  over  a  three 
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dimensional  element  in  the  frequency  domain.  Details  of  this  method  can 
be  found  in  [63],  In  general,  this  method  is  computationally  more  demand¬ 
ing  than  the  above  methods.  Thompson  [63]  claims  that  frequency  domain 
averaging  over  a  three  dimensional  element  is  less  sensitive  to  nonstation- 
arity  in  the  data  than  the  conventional  methods.  The  bispectral  estimates 
derived  using  Thompson’s  multiwindow  techniques  are  consistent,  unlike 
the  direct  method  discussed  above  [63],  where  consistency  means  that  both 
the  variance  and  the  bias  of  estimate  approach  zero  as  the  number  of  data 
points  increases  [60].  It  is  shown  in  [12]  that  the  bias  in  the  estimates 
obtained  using  the  direct  method  does  not  approach  zero  as  N  — »  oo. 
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Chapter  3 

THE  BIRANGE  PROFILE 
OF  RADAR  SIGNATURES 


3.1  Introduction 

The  focus  in  this  chapter  is  on  interpreting  and  understanding  the  bis- 
pectral  features  of  radar  targets.  The  goal  is  to  provide  an  appropriate 
interpretation  of  the  bispectrum  of  radar  signatures,  and  to  identify  the 
advantages  and  disadvantages  of  this  type  of  radar  signal  processing.  To 
this  end,  we  first  investigate  the  role  of  the  bispectrum  in  achieving  a  better 
understanding  of  the  scattering  mechanisms  of  radar  targets.  Secondly,  we 
extract  features  that  characterize  the  illuminated  target.  An  investigation 
of  the  role  of  the  bispectrum  in  identifying  radar  targets  is  presented  in 
Chapter  V. 

The  type  of  signal  processing  addressed  in  this  study  is  based  on  a 
specified  target  scattering  model.  The  model  is  that  of  equation  (2.10), 
where  the  scattered  signal  is  a  combination  of  specular  scattering  terms 
from  localized  areas  on  the  target  and  multiple  interactions  terms.  The 
scattering  region  of  interest  is  the  resonance  region  (in  which  the  target  size 
is  on  the  order  of  a  wavelength).  The  bispectral  signatures  of  measured  (or 
computed)  scattering  coefficients  in  the  optical  region  are  also  examined. 

The  purpose  of  bispectral  analysis  of  radar  signals  is  to  add  to  the 
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knowledge  that  the  transient  response  already  offers.  It  is  important  to 
emphasize  that  the  bispectrum  is  not  an  alternative  to  conventional  radar 
signal  processing  (i.e.,  estimating  the  impulse  response)  but  an  addition 
that  may  be  useful  depending  on  the  prevailing  circumstances  such  as  data 
conditions,  noise  statistics,  and  scattering  region. 

It  is  known  that  the  scattering  mechanisms  of  complex  structures  such 
as  airplanes  and  ships  are  far  more  complicated  than  what  can  be  deduced 
from  the  target  transient  response.  The  bispectral  responses  of  such  compli¬ 
cated  targets  are  also  difficult  to  analyze.  For  this  reason,  it  is  constructive 
to  first  consider  synthesized  targets  and  simple  objects.  Thus,  the  bispec¬ 
tra  of  a  synthesized  target,  a  simple  generic  aircraft  model  and  a  tilted  flat 
plate  are  first  examined.  The  bispectral  responses  of  aircraft  models  are 
also  investigated. 

Bispectral  analysis  of  radar  signatures,  as  implemented  in  this  study, 
has  two  major  applications:  1)  identification  of  multiple  interactions,  and 
2)  partial  noise  suppression.  This  chapter  is  primarily  concerned  with  the 
interpretations  of  target  bispectra  in  terms  of  modeled  scattering  mecha¬ 
nisms  where  the  model  in  (2.10)  serves  as  a  basis  for  the  interpretation. 
The  effect  of  additive  noise  on  the  bispectral  response  of  radar  signatures 
is  also  investigated. 

3.2  The  Birange  Profile:  Definition  and  Prop 
erties 

Target  scattering  features  are  more  intuitive  using  time-domain  signatures  [2] 
than  using  frequency  domain  responses.  Further,  the  transient  response  is 
an  important  representation  of  time-domain  radar  signatures.  However, 
it  is  often  difficult  to  deduce  accurate  information  about  the  target  shape 
using  only  the  transient  response  because  of  the  presence  of  multiple  in- 
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teractions.  In  the  following,  the  definition  of  the  bispectrum  in  a  radar 
signature  analysis  framework  is  given.  The  bispectrum  of  radar  signals,  as 
defined  in  this  Chapter,  can  be  used  as  a  tool  for  discriminating  between 
multiple  interactions  and  specular  scattering. 

In  most  signal  processing  problems,  the  data  sequence  consists  of  sam¬ 
ples  taken  from  a  time-dependent  waveform.  Then  the  bispectrum,  as  de¬ 
fined  in  Chapter  II,  is  a  function  of  two  variables  oji  and  u>2  in  the  frequency 
domain.  The  radar  scattering  data  sequence,  however,  is  not  a  time  series 
but  recorded  in  the  frequency  domain.  The  triple  correlation  of  the  fre¬ 
quency  domain  data  is  a  two-dimensional  profile  in  the  frequency  domain 
and  is  defined  as 

i?(wi, W2)  =  E{H  {w)H (w  +  (w  +  W2)}  (3.1) 


The  bispectrum  is  then  given  as 

OO  OO 

£(<i,<2)=  £  H  i?(wi,w2)exp{-y'(w1<i +w2<2)}  (3.2) 

U>1  =  —  OO  U>2  =  —  OO 

Therefore,  the  bispectrum  is  a  profile  of  target  scattering  signatures  in 
the  two-dimensional  time  domain  and  can  be  termed  the  “bi-time”  pro¬ 
file  of  radar  targets.  Further,  the  bispectrum  can  be  expressed  as  a  two- 
dimensional  profile  in  the  range  domain  using  the  relation 


(3.3) 


where  r  is  the  range  from  the  radar  to  the  target,  t  is  the  time  needed  for 
the  signal  to  propagate  to  the  target  and  back,  and  c  is  the  speed  of  light. 
Therefore,  the  bispectrutn  as  a  function  of  range  is  given  as 


B{rur2)  =  SjR(/i,/2)exp{-i47r/c(/1r1  +  far?.)}  (3.4) 

fi  h 

where  /,•  =  a>;/27r.  This  representation  of  the  bispectrum  corresponds  to 
a  birange  profile  of  the  target.  Hence,  the  birange  profile  is  a  bispectral 
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display  of  target  signatures  in  the  two-dimensional  range  domain  ri  and  r2 
where  the  term  “range”  denotes  the  propagation  distance  along  the  target. 

The  birange  profile  of  complex  radar  signals  satisfies  only  one  symmetry 
property  defined  as 

B(rur2)  =  5(r2,ri).  (3.5) 

The  definition  of  the  birange  is  more  insightful  when  expressed  in  terms  of 
the  impulse  response  as 

B{rur2)  =  {h{ri)h(r3)h*(ri  +  r2))  (3.6) 

where  h(r )  is  the  impulse  response  as  a  function  of  range  and  ()  indicates 
the  ensemble  average.  Defining  the  birange  as  such  implies  that  a  non¬ 
zero  bispectral  response  is  the  result  of  interaction  between  the  responses 
at  rx  and  r2  which  appears  as  a  response  at  ri  +  r2.  Furthermore,  the 
bispectrum  at  (ri,r2)  is  non-zero  only  if  the  responses  at  ri  and  r2  and 
ri  4-  r2  are  correlated.  Notice  that  the  definition  of  the  birange  as  given  in 
(3.4)  corresponds  to  the  case  where  the  triple  correlation  is  computed  using 
(3.1);  i.e.  the  conjugate  is  placed  on  the  first  element  of  the  triple  product. 

3.2.1  Estimation  of  The  Birange  Profile 

The  classical  bispectrum  estimation  techniques,  summarized  in  Chapter  II 
and  slightly  modified,  can  be  used  to  estimate  the  birange  profile.  The 
modifications  are  due  to  data  conditions  such  as  stationarity  and  data- 
length.  Stationarity  requires  that  the  third  order  cumulant  of  the  measured 
data  be  invariant  with  respect  to  the  frequency  band  used.  Also,  having  a 
finite  set  of  data  imposes  some  constraints  on  segmenting  the  data  without 
reducing  the  resolution  in  the  birange.  The  following  bispectral  estimation 
algorithm  is  used  throughout  this  chapter. 

The  birange  profile  can  be  estimated  using  either  the  direct  method  or 
the  indirect  method  (based  on  the  triple  correlation  of  the  data).  Estimates 
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of  the  birange  profile  in  the  experimental  part  of  the  study  are  computed 
using  a  slightly  modified  version  of  the  indirect  classical  method.  The 
algorithm  used  for  birange  estimation  given  a  finite  set  of  data  {#(0;)},  is 
outlined  below: 

1)  Subtract  the  average  value  of  the  data.  The  purpose  of  this  step  is  to 
have  a  zero  mean  signal,  which  is  needed  to  achieve  noise  suppression. 

2)  Estimate  the  triple  correlation 

R{m , n)  =  -1  £  H*(l)H{l  +  m)H(l  +  n)  (3.7) 

v  i=»  1 

where  si  =  max  (0,  — m,  — n)  and  $2  =  min  (N  —  1,1V  —  1  —  m,  N  —  1  —  n). 
The  number  of  data  points  is  denoted  by  N,  and  —L  <  m,n  <  L. 

3)  The  birange  is  generated  using 

L  L 

B{r1,r2)=  ^(rn.nJ^m.nJexpl-jXwim  +  wjn)}  (3.8) 

m=—L n=-L 

where  W(m,n)  is  the  window  function  defined  in  (2.62)  and  where  L  is 
the  triple  correlation  lag  (usually  chosen  to  be  «  IV/ 5).  Because  no  data 
segmentation  is  employed,  high  variance  birange  profiles  are  obtained  par¬ 
ticularly  for  large  values  of  L.  The  reason  for  avoiding  segmentation  is 
that  the  data  may  not  be  stationary  (because  target  scattering  features 
are  dependent  on  the  frequency  band  used)  and  segmentation  introduces 
potential  non-station  arities  [11]  and  reduces  the  resolution  of  the  bispectral 
estimates  [60], 

The  assumption  that  the  scattering  data  are  stationary  requires  that 
the  triple  correlation  is  dependent  only  on  the  frequency  increment,  that  is 

E{X^(Ji)X(fj)X(Jk)}  =  E{XV)X(f  +  fj-fi)X(f  +  ,fk -/,•)}  V/  (3.9) 

In  radar  scattering,  this  assumption  is  not  valid  over  a  wide  frequency  range. 
Therefore,  the  birange  5(ri,r2)  is,  in  general,  dependent  on  frequency.  It  is 
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possible,  however,  to  assume  that  the  birange  does  not  change  significantly 
over  small  frequency  bands.  Further,  by  avoiding  data  segmentation  and 
large  data  bandwidths,  it  is  possible  to  obtain  an  estimate  of  the  birange 
profile  which  is  approximately  frequency  independent. 

The  definition  of  the  triple  correlation  as  given  in  (3.1)  requires  that  the 
measured  data  have  zero  mean  in  order  to  exclude  second  order  moments 
from  the  definition  of  the  third  order  cumulant.  This  constraint  is  par¬ 
ticularly  useful  to  ensure  suppression  of  additive  Gaussian  noise.  There  is 
no  guarantee,  however,  that  this  condition  is  satisfied.  One  can  avoid  this 
problem  by  simply  removing  the  mean  of  the  data.  But,  when  the  sample 
mean  is  removed  from  the  data,  the  birange  may  erroneously  indicate  the 
presence  of  a  response  at  time  zero  [60].  Furthermore,  the  scattering  data 
may  include  transient  signals  and  removal  of  sample  mean  is  undesirable  [60] 
in  this  situation.  On  the  other  hand,  if  the  average  is  not  subtracted  from 
the  data,  then  the  third  order  cumulant  differs  from  third  order  moment 
and  will  include  the  second  order  moment  of  the  data  sequence.  Further, 
the  addition  of  the  second  order  moment  to  the  cumulant  introduces  unde¬ 
sirable  noise  effects  because  the  second  moment  of  additive  Gaussian  noise 
is  not  zero  (see  details  in  Section  3.8).  Thus,  for  the  radar  application  the 
benefits  of  mean  removal  seem  to  outweigh  the  disadvantages. 

3.3  Interpretation  of  The  Bispectrum  of  Radar 
Signatures 

An  interpretation  of  the  bispectrum  of  radar  signatures  based  on  an  as¬ 
sumed  scattering  model  is  given  in  this  section.  The  assumed  model  corre¬ 
sponds  to  scattering  in  the  resonance  region,  and  it  is  assumed  the  scattered 
signal  is  composed  of  two  major  components.  The  first  component  is  due 
to  specular  direct  scattering  from  localized  positions  on  the  target.  The 
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scattered  signal  due  to  the  first  component  is  given  as 

m 

=  £  ^PM  exP  +  w<p)}  (3.io) 

p=i 

where  m  is  the  number  of  scattering  centers  and  Ap  is  the  amplitude  of 
the  scattered  signal  from  the  pth  scattering  center.  The  electric  phase  of 
the  pth  scatterer  is  denoted  by  and  tp  is  the  propagation  time  between 
the  radar  and  the  pth  scattering  subcomponent  and  back.  Notice  that  the 
inverse  Fourier  transform  of  the  above  returned  signal  is  given  by 

m 

hi(t)  =  J2AP  exP  i-M’pW  -  tp)  (3>n) 

P=1 

assuming  that  the  amplitude  Ap  and  phase  tJjp  are  independent  of  frequency. 
The  impulse  response  h\(t)  thus  consists  of  m  impulses  displayed  at  times 
fp,Vp  =  l,...,m.  Therefore,  each  peak  in  the  impulse  response  is  an  in¬ 
dication  of  specular  scattering  at  time  tp.  If  either  the  amplitude  Ap  or 
the  phase  is  a  function  of  frequency  then  the  corresponding  peak  may 
be  displaced  from  tp  and  may  have  nonzero  width;  this  phenomenon  is  a 
characteristic  of  frequency  dispersive  scattering. 

The  second  major  component  that  contributes  to  target  scattering  is 
the  effect  of  multiple  interactions.  This  type  of  scattering  occurs  in  sev¬ 
eral  ways.  The  interactions  may  include  multiple  propagation  between  two 
or  three  scattering  centers.  Multiple  interactions  may  also  occur  between 
scatterers  that  are  not  collinear  with  the  radar  axis.  In  general,  there  is  no 
defined  rule  for  multiple  interactions  in  terms  of  direction  or  magnitude. 

The  multiple  interactions  can  be  modeled  as 

m  n 

=  ^SyiP9(a;)eXP{-i(V’p  +  V,q+w(fp  +  f5))}  (3.12) 

p~i 5-i 

m  n 

+  £  E  A'n(u)  exP  {-rt^P  +  V’q  +  W<w)} 

p=l 9=1 
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where  Apq(c«;),  i4pg(u>)  indicate  the  amount  of  coupling  between  the  pth  and 
qth  scattering  points.  The  distinction  between  Apq  and  A '  is  needed  be¬ 
cause  some  interactions  of  magnitude  A ^  are  not  necessarily  scattered  at 
tp  -f-  tq.  The  inverse  Fourier  transform  of  H2(u)  is  given  by 

mm 

M<)  =  E  E  eXP  {— i(V»p  +  ^q)}  {Apq6(t  -  {tp  +  tq))  +  A'pqS{t  ~  tpq)) 

P=  1  9=1 

(3.13) 

and  consists  of  impulses  at  times  tp  +  tq,  and  and  p,  q  =  1, . . .  ,m  assum¬ 
ing  that  Apq,  if}p,  and  are  independent  of  frequency.  Furthermore,  the 
relations  between  Apq  or  A'pq  and  the  pair  (.Ap,  Aq)  are  not  usually  known. 
Notice  that  while  some  interactions  (indicated  by  multiple  propagation  of 
the  transmitted  radar  signal  between  scatterers)  are  received  by  the  radar 
at  tp  +  tq,  other  types  of  interactions  (mainly  due  to  single  propagation 
between  scatterers)  are  received  at  time  tpq  ^  tp  +  tq.  Both  forms  of  inter¬ 
actions  are  modeled  in  H2{w). 

A  scattering  model  which  includes  both  specular  and  interaction  terms 
is  a  combination  of  H i(w)  and  H2(w).  Hence,  the  returned  backscattered 
radar  signal  can  be  modeled  as 

H{ ;«)  =  Hx(w)  +  tf2(u>)  (3.14) 

m 

=  E^pMexP{-:^p  +  w*p)}  (3-15) 

P=1 

m  n 

+  E  E  eXP  {-i(VV  +  i>q+  v{tp  +  tq)} 

P= 1 9=1 
m  n 

+  E  E  Apg(W)  eXP  {-J'(^P  +  ^9  +  “tpq) 

p=l g=l 

Based  on  the  above  model,  which  is  more  realistic  than  using  Hi(lo) 
only,  it  is  clear  that  conventional  spectral  analysis  does  not  provide  the 
necessary  distinction  between  specular  scattering  and  multiple  interactions. 
Also,  the  Fourier  transform  or  modern  spectral  estimation  techniques  of 
the  radar  backscatter  data  cannot  identify  the  coupling  between  interacting 
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subcomponents.  Using  the  definition  of  the  bispectrum  given  in  Section  I  of 
this  Chapter,  the  bispectral  response  of  a  radar  target  whose  backscattered 
signal  satisfies  (3.15)  is 

B{tut2)  =  {h{t1)h{t2)h‘{tl+t2))  (3.16) 

m  m 

=  —  ^2  -  tq)  4-  ^(<1  —  tq)  t2  —  tp)} 

p=lg=l 

that  is,  the  birange  profile  of  H(u>)  consists  of  impulses  at  the  pairs  (tp,tq), 
p,  q  =  1, . . .  ,m  with  magnitudes  ApAqAp+q.  Hence,  the  bispectrum  of  radar 
signals  will  detect  the  interactions  with  Apq  terms,  but  will  not  detect  inter¬ 
actions  with  A'pq  terms.  In  otherwords,  the  birange  profile  detects  multiple 
interaction  terms  when  the  time  delay  is  the  sum  of  the  two  time  delays  of 
the  pth  and  qth  specular  scattering  terms.  Based  on  the  above  modeling  of 
radar  signatures,  it  is  reasonable  to  say  that  peaks  in  the  birange  (or  bis¬ 
pectrum)  of  radar  signals  are  indications  of  multiple  interactions  between 
scattering  subcomponents.  Furthermore,  these  interactions  occur  between 
subcomponents  whose  distance  to  the  radar  is  determined  by  the  coordi¬ 
nates  of  the  peaks  in  the  birange. 

The  intuition  behind  the  identification  of  these  interactions  is  that 
the  birange  profile  B(ri,r2)  (or  bi-time  B(li,t2))  is  related  to  the  im¬ 
pulse  response  using  the  Cramer  spectral  representation  (2.41)  B(ti,t2)  oc 
(h(ti)h(t2)h*(ti  +  i2))-  Therefore,  the  birange  exhibits  peaks  at  (U,ij)  when 
the  impulse  response  has  peaks  at  f,-,  tj ,  and  U+tj,  and  these  peaks  are  cor¬ 
related.  There  are  two  situations  in  which  the  impulse  response  may  peak 
at  these  points  simultaneously:  1)  when  a  multiple  interaction  of  the  form 
given  in  (3.15)  occurs,  or  2)  by  coincidence  when  there  are  specular  scat¬ 
tering  components  at  all  three  time  instants.  A  method  for  discriminating 
between  these  two  possibilities  is  discussed  in  Section  (3.3.3)  below. 

The  above  interpretation  is  valid  for  the  interior  regions  of  the  birange 
but  not  for  the  responses  at  the  axes.  In  fact  the  responses  on  the  time  (or 
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range)  axes  may  not  correspond  to  any  interactions.  The  bitime  profile  on 
the  time  axes  is  given  by 

B(0,i,)  =  Bfo.O)  =  {h(0)\h(t)\2)  (3.17) 

Notice  that  the  responses  on  the  axes  of  the  birange  are  proportional  to 
the  power  spectral  density  of  the  target  signatures.  Therefore,  the  infor¬ 
mation  provided  by  the  bi-time  profile  is  twofold:  1)  multiple  interactions 
are  displayed  in  the  internal  region  between  the  axes,  and  2)  the  square  of 
the  magnitude  of  the  impulse  response  is  displayed  along  the  time  axes.  If 
the  birange  is  estimated  using  only  the  real  component  of  the  radar  signal 
then  the  magnitude  squared  of  the  impulse  response  will  be  also  displayed 
along  the  line  <i  +  <2  =  0. 

The  interpretation  of  the  bispectrum  of  radar  signatures  as  given  above 
is  clarified  by  the  following  two  examples.  The  identification  of  multiple 
interactions  of  real  radar  targets  using  the  bispectrum  is  also  demonstrated 
using  several  examples  of  experimental  radar  scattering  data. 

3.3.1  Example  1:  Synthesized  Target 

The  first  example  consists  of  scattering  from  a  hypothetical  three  point- 
scatterer  hypothetical  target.  The  purpose  of  this  example  is  to  present  an 
intuitive  interpretation  of  the  role  of  the  bispectrum  in  identifying  multiple 
interactions  as  discussed  in  the  previous  section.  For  conceptual  purposes 
only,  the  first  two  point-scatterers  correspond  to  the  leading  and  the  trail¬ 
ing  edges  of  a  blade,  and  the  third  scatterer  corresponds  to  the  specular 
scattering  mechanism  from  a  sphere. 

Consider,  as  shown  in  Figure  3.1,  a  blade  of  width  T\  —  r0  and  a  sphere 
at  a  distance  r2  —  r\  from  the  trailing  edge  of  the  blade.  If  one  considers 
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only  four  ray  paths,  the  backscattered  signal  (assuming  r0  =  0)  is 

St(fj)  =  Ay  +  ^.2  exp 

+Ai  exp  j47r(ri  + 

We  have  assumed  that  the  amplitudes  A^,A2,A3,  and  A.4  (see  Figure  3.1 
for  the  values  of  A.i,.A2,.A3,  and  Aj)  are  frequency  independent  and  that 
the  sphere  creeping  wave  terms  are  negligible. 

The  scattered  signal  from  this  target  satisfies  the  complex  exponential 
model  of  (3.15).  The  transient  response  of  the  blade  and  sphere  target  as 
computed  using  Fourier  transform  with  a  Hanning  window  over  a  frequency 
band  from  1.5  GHz  to  12  GHz  is  shown  in  Figure  3.2.  Note  the  one  to 
one  correspondence  of  the  first  three  responses  with  the  geometry  of  the 
target.  The  fourth  response,  however,  is  an  interaction  term  and  does  not 
correspond  to  any  target  components.  The  birange  of  the  above  signal  in 
the  first  symmetry  region  (defined  by  tx,t2  >  0, +  i2  =  1/2A /  where  A / 
is  the  sampling  period)  is 

Bst{Ri,Rz)  =  \{A316(RuR2)  +  A1A226(R1-r1,R2)  (3.19) 

+J41^.3^(i?1  —  r2,  R2)  4-  A\  A\8(Ri  —  (rj  +  r2),  R2) 
-\rA2A2A<i8(Ri  —  r2,  R2  —  rj)}. 

Figure  3.2  shows  this  birange  profile  (with  no  windowing).  Notice  the 
response  at  (fi,f2)  =  (1.6, 2.4)  meters  which  is  due  to  interaction  between 
the  trailing  edge  of  the  blade  and  the  sphere  and  is  represented  by  the 
term  exp  j47r(r;  +r2  —  r0)^).  Further,  the  response  on  the  range  axes 
rq  =  0,  r2  =  0  and  also  along  the  line  rx  +  r2  =  0  is  proportional  to  the 
magnitude  squared  of  the  transient  response.  Therefore,  the  birange  of 
the  blade-sphere  target  includes  information  about  multiple  interactions  at 
(ri,r2)  as  well  as  the  magnitude  squared  of  the  transient  response. 


^  j  +  A3  exp 


f  f.y 

—j4nr2— 

,  c  t 


r2) 


(3.18) 
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Tapering  the  triple  correlation  function  of  the  frequency  response  using  a 
two  dimensional-window  (to  avoid  Gibbs  sidelobe  phenomenon)  [14]  results 
in  the  bispectral  response  shown  in  Figure  3.3.  The  window  function  used  is 
the  Minimum  Bias  Supremum  given  in  (11).  Note  that  each  impulse  in  the 
original  bispectrum  is  spread  into  an  ellipse  corresponding  to  the  contour 
of  the  window  response  in  the  time  domain  [14]. 


3.3.2  Example  2:  Thin-Blade  &  Sphere  Target 


The  scattering  model  in  this  case  is  similar  to  the  previous  example.  The 
width  of  the  blade  r j  —  r0,  (or  the  distance  between  the  first  two  point- 
scatterers),  is  reduced  such  that 

2(7^1)  (3'20) 

where  —  B  which  is  the  measurement  bandwidth.  The  resolution 

A t  of  the  transient  response  when  estimated  using  Fourier  transform  tech¬ 
niques  is 

A(  =  i  (3.21) 

which  is  equivalent  to  range  resolution 


(3.22) 


Therefore,  due  to  this  choice  of  blade  width,  the  first  two  peaks  of  the  tran¬ 
sient  response  cannot  be  separated  by  Fc  *  rier  processing  because  of  resolu¬ 
tion  limitations.  The  transient  response  and  the  birange  of  this  target  are 
shown  in  Figure  3.4.  The  first  peak  in  the  transient  response  corresponds 
to  a  superposition  of  direct  scattering  from  both  leading  and  trailing  edges 
of  the  thin  blade.  The  second  peak  corresponds  to  a  superposition  of  direct 


scattering  from  the  sphere  and  the  interaction  term.  Thus,  due  to  the  lim¬ 
ited  resolution  of  the  transient  response,  we  are  unable  to  detect  the  two 


edges  of  the  blade  or  the  interaction  term. 
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Now  consider  the  birange  profile  of  this  target  (see  Figure  3.4).  Notice 
that  the  birange  exhibits  a  peak  at  (r!,r2)  as  an  indication  of  multiple  in¬ 
teraction.  Also,  notice  that  the  magnitude  square  of  the  transient  response 
is  no  longer  displayed  along  the  axes.  The  reason  is  that  the  magnitude 
of  the  birange  at  the  range  axis  i?(0,r)  or  B(r,  0)  is  proportional  to  the 
transient  response  at  t  =  0,  and  h(0)  =  0. 

Two  conclusions  can  be  drawn  from  this  example:  1)  the  birange  may 
detect  interactions  that  may  not  be  displayed  in  the  impulse  response,  per¬ 
haps  for  resolution  reasons,  2)  the  magnitude  squared  of  the  transient  re¬ 
sponse  is  displayed  along  the  range  axes  only  if  h(0)  ^  0  (if  the  data  is 
segmented,  then  h(0)  ^  0  is  not  sufficient  to  obtain  the  spectral  density 
along  the  axes  of  the  birange).  Hence,  this  example  shows  that  neither 
the  transient  response  nor  the  birange  is  sufficient  to  analyze  scattering 
mechanisms.  A  combination  of  both  techniques  is  the  key  to  time-domain 
scattering  analysis. 

3.3.3  A  Method  for  Detecting  False  Interactions 

A  peak  at  r*  +  r2  in  the  transient  response  is  possibly  an  indication  of 
interactions  between  two  scattering  subcomponents  at  r\  and  r2  along  the 
target.  An  alternative  scenario  however,  is  that  the  peak  at  ra  +  r2  is  due 
to  two  components;  one  is  specular  which  is  coincidently  located  at  r\  -f  r2 
and  the  other  is  the  result  of  multiple  interaction  between  the  scatterers  at 
r\  and  r2.  The  problem  is,  then,  how  to  determine  whether  the  response 
at  ri  +  r2  is  specular,  or  is  the  result  of  an  interaction,  or  is  a  combination 
of  the  two.  This  can  be  done  by  considering  the  bicoherence  at  the  peak. 
The  bicoherence  at  (ri,r2)  is  simply  the  birange  at  (ri,r2),  normalized  with 
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Figure  3.3:  Birange  profile  of  blade-sphere  target  (windowed  triple  correla¬ 
tion) 
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respect  to  the  spectral  density  at  r1}r2  and  rx  -f  r2; 

bicH(rur2)  =  BH(rXir2)  -  (3.23) 

\JSH(ri)SH(r2)SH(ri  +  r2) 

where  SH(r )  is  the  spectral  density  of  the  target  response  H(f )  at  r.  If 
the  bicoherence  at  (ri,r2)  is  equal  to  one,  then  the  response  at  ri  +  r2 
is  totally  caused  by  an  interaction.  If  the  bicoherence  is  less  than  one, 
then  the  response  at  ri  +  r2  is  due  to  a  partially  correlated  interaction,  or 
due  to  both  a  multiple  interaction  and  a  specular  scattering  term.  If  the 
bicoherence  is  zero  then  there  is  no  interaction  at  rj  +  r2  and  the  response 
is  due  to  specular  scattering. 

As  an  example,  consider  the  specular  scatterers  at  rp,  rq  and  rp  +  rq. 
Assume  that  there  is  a  multiple  interaction  between  rp  and  rq  that  produces 
a  response  at  rp+rq  in  the  transient  response.  Then,  the  frequency  response 
is  given  as 

H{u>)  =  Ap  exp  {— j‘2— (rp)  +  ^P}  +  Aq  exp  {-j2-(rq)  +  if>q} 

c  c 

+Apq  exp  {— j2— (rp  +  rq)  +  ipp  +  V><,} 

+Ap+q  exp  {—j2—(rp  +  rq)  +  ^p+,}  (3.24) 

The  bispectrum  of  H(u)  is  then  given  as 

BH{rur2)  =  ApAqApq  (£(ra  -  rp,r2  -  rq)  +  -  rq,r2  -  rp ))  (3.25) 


and  the  spectral  density  of  H(f)  is 

SH(r)  ■-  A2p6(r  -  rp)  +  A2qS(r  -  rq) 

+  (^p9  +  ^p+q)^(r  _  (rP  +  rg )) 


(3.26) 
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(3.27) 
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where  Apq  is  the  magnitude  of  the  interaction  between  Ap  and  Aq,  while 
Ap+q  is  the  magnitude  of  direct  specular  scattering  at  rp  +  rq.  Notice  that 
if  Ap+q  =  0,  then,  the  bicoherence  is  equal  to  one.  If  there  is  no  interaction 
between  scatterers  at  rp  and  r9,  then  Apq  =  0  and  the  bicoherence  is  equal 
to  zero. 


3.3.4  Example  3: 


Consider  the  canonical  scatterer  shown  in  Figure  3.5.  Assume  that  the 
radar  returns  for  this  target  satisfy  (3.15),  (i.e.,  they  can  be  modeled  as 
a  sum  of  complex  exponentials).  Further,  assume  that  the  backscatter 
signals  of  magnitudes  A  and  B  correspond  to  scattering  from  the  first  and 
the  second  edge  of  the  arc,  respectively.  Also,  assume  that  the  backscatter 
signal  of  magnitude  C  represents  an  interaction  term  and  corresponds  to 
propagation  to  the  first  edge  then  to  the  second  edge  and  back  to  the  radar. 
The  returned  signal  of  magnitude  D  corresponds  to  specular  scattering  from 
the  tip  of  the  flat  section  connected  to  the  arc.  The  geometry  of  this  target 
is  chosen  such  that  the  time  of  arrival  of  the  multiple  interaction  coincides 
with  the  backscatter  from  the  tip  of  the  flat  edge  (see  Figure  3.5  where 
the  dimensions  are  given  in  units  of  propagation  time).  The  total  returned 
signal  is  given  as  (where  u>  =  2i r/,  and  /  is  in  GHz) 

H(u>)  =  Aexp{-j(i>i  +  4u>)}  +  Bexp{-j(rj)2  +  16u?)} 

+C  exp  +i>2  +  20u>)}  +  Dexp{-j(ipz  +  20u>)} 


The  transient  response  of  this  target  has  three  peaks;  the  first  peak  at  4  ns 
of  magnitude  |A|,  the  second  peak  at  16  ns  of  magnitude  |Z?|,  and  the  third 
peak  at  20  ns  of  magnitude  | C  +  D\.  The  birange  of  this  target  shows  a 
peak  at  (4,18),  and  (18,4)  of  magnitude  \ABC\.  The  peak  in  the  birange 
is  an  indication  of  an  interaction  between  the  two  edges  of  the  arc.  This 
does  not  imply,  however,  that  the  peak  at  20  ns  in  the  transient  response 
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Figure  3.5:  A  canonical  target  whose  specular  scattering  component  of  mag¬ 
nitude  D  is  received  by  the  radar  in  the  same  time  instant  as  the  interaction 
of  magnitude  C. 
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is  totally  due  to  multiple  interactions.  Moreover,  the  bicoherence  of  the 
peak  at  (4,16)  ns  in  the  birange  profile  is  G/\/G2  +  D 2  which  indicates 
the  relative  magnitude  of  the  response  due  to  the  interaction  compared  to 
the  total  magnitude  of  the  response  at  20  ns  as  displayed  in  the  transient 
response.  This  example  shows  that  the  bicoherences  of  the  peaks  in  the 
birange  may  be  used  to  separate  the  responses  due  to  multiple  interactions 
from  the  responses  due  to  specular  scattering  when  they  coincide  in  the 
transient  response. 

Using  the  bicoherence  to  separate  between  interactions  and  specular 
scattering  is  not  as  accurate  as  discussed  above.  The  estimate  of  the  bico¬ 
herence  is  accurate  if  the  scatterin;  '  •'m  a  radar  target  satisfies  the  model 
of  equation  (3.15).  It  is  often  the  case,  however,  that  the  scattering  is 
dispersive,  so  that  the  amplitude  of  specular  scattering  or  interactions  is 
dependent  on  frequency.  Then  the  bicoherence  of  a  pure  interaction  at 
( rp,rq )  will  not  equal  one.  Furthermore,  if  neither  interaction  nor  specular 
scattering  occurs  at  (rp,rq),  the  bicoherence  may  not  equal  zero,  because 
of  residual  values  of  the  birange  and  spectral  density  at  rp  +  rT  Neverthe¬ 
less,  the  bicoherence  of  radar  targets  can  be  used  to  quantify  the  level  of 
interactions  between  scatterers. 

3.4  The  Birange  of  Dispersive  Scatterers 

The  scattering  model  given  in  (3.15),  on  which  the  interpretation  of  the  bi¬ 
range  is  based,  does  not  include  scattering  from  dispersive  subcomponents. 
In  the  following  section  we  discuss  the  effect  of  scattering  from  dispersive 
targets  on  the  interpretation  of  the  birange. 

Physically,  dispersion  occurs  when  the  interrogating  signal  (transmitted 
by  the  radar)  loses  energy  gradually  while  propagating  along  the  target.  A 
non-dispersive  target  scatters  energy  strongly  at  short  time  intervals  rela- 
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tive  to  the  propagation  time  along  the  target.  A  dispersive  target  scatters 
energy  with  lower  magnitude  but  over  longer  periods  of  time. 

An  example  of  dispersive  scattering  behaviour  is  when  the  amplitude 
A(iv)  changes  with  frequency  as 

A(to)  -  A0exp{-j~}  (3.28) 

where  p  is  a  constant.  The  measured  frequency  response  in  this  case  is 
Hd(u)  =  exp  {—j—}H(u>) 

U) 

m 

=  2  Av  exp  {-jw(~J  +  ip)  -  jtpp}  (3.29) 

p=i  w 

by  ignoring  the  higher  order  terms.  Notice  that,  with  this  form  of  dis¬ 
persion,  the  signal  that  scatters  at  time  tp  for  the  non-dispersive  case  will 
scatter  gradually  as  a  function  of  w.  For  high  frequencies,  scattering  oc¬ 
curs  at  time  tp  approximately,  but  as  the  frequency  decreases  the  scattering 
center  is  displaced  away  from  tp.  This  corresponds  to  gradual  energy  loss 
along  the  target  when  illuminated  by  the  radar.  This  type  of  dispersion  has 
similar  effects  on  both  the  transient  response  and  the  birange  profile.  Peaks 
in  the  transient  response  and  the  birange  become  less  localized  because  of 
the  gradual  energy  loss.  In  a  special  case,  however,  if  an  interaction  oc¬ 
curs  between  tp  and  tq  where  the  the  dispersion  coefficients  are  such  that 
pp  =  —  pq,  then  the  positive  displacement  at  tp  is  compensated  by  the  neg¬ 
ative  displacement  at  tq ,  without  displacing  the  position  of  the  interaction 
term  at  (tp,tg). 

Windowing  the  measured  scattering  data  and  the  triple  correlation  of 
the  data  may  reduce  the  effect  of  this  type  of  dispersion  on  the  transient 
response  and  the  birange  profile.  By  windowing  the  frequency  response 
of  the  target,  the  peaks  corresponding  to  non-dispersive  scatterers  are  dis¬ 
played  over  W  seconds  (where  W  is  the  width  of  the  main-lobe  of  the 
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window  used).  The  peaks  corresponding  to  dispersive  scatterers  (of  disper¬ 
sion  width  8  seconds),  however,  are  displayed  over  W  +  8  seconds.  Thus, 
by  choosing  W  >  8,  it  is  possible  to  reduce  the  disagreement  between  the 
signature  of  a  dispersive  scatterer  and  that  of  a  non-dispersive  one. 

As  another  example,  consider  the  scattering  from  a  small  conducting 
sphere  as  compared  to  scattering  from  a  point  scatterer.  The  amplitude 
A(u;)  of  the  scattered  signal  from  a  sphere,  when  measured  in  the  resonance 
region,  changes  with  frequency  as  [45] 

j4(o>)  =  A0  (1  +  a  cos(/3u;))  (3.30) 


where  a,/3,  and  Aq  are  constants.  Such  a  behaviour  is  a  result  of  what 
is  commonly  known  as  the  creeping  wave.  The  transient  response  of  the 
sphere  peaks  in  two  positions.  The  first  peak  corresponds  to  direct  specular 
scattering  from  the  sphere;  the  second  peak  is  caused  by  the  creeping  wave. 
Therefore,  the  dispersive  behaviour  of  the  sphere  as  compared  to  a  point 
target  is  characterized  by  an  additional  peak  in  the  transient  response. 

Assume  that  the  three-point  scatterer  target  in  Example  1  is  replaced  by 
a  three-sphere  target  where  the  amplitudes  A(jw)  are  frequency  dependent 
as  given  in  (3.30).  The  transient  response  and  the  birange  of  this  target  are 
shown  in  Figure  3.6  for  a  =  .6,(3  =  .01,  and  A0  =  1  .  Notice  that  there  are 


more  peaks  in  the  transient  response  due  to  creeping  waves  around  each  of 
the  three  spheres.  The  birange,  however,  is  unaffected  by  such  dispersive 
behaviour  except  for  changes  along  the  axes  (where  the  magnitude  squared 
of  the  transient  response  is  displayed).  Hence,  while  the  transient  response 
of  the  dispersive  target  has  changed,  the  number  of  multiple  interactions 
in  the  birange  remains  fixed.  The  reason  is  that,  in  this  example,  the 
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of  the  target  for  the  creeping  waves  to  introduce,  by  coincidence,  false 
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responses  in  the  birange.  In  general,  unmodeled  dispersive  behaviour  cause 
changes  in  both  the  transient  response  and  the  birange  profile.  The  birange, 
however,  is  less  sensitive  to  these  changes  than  the  transient  response. 

3.5  Examples  of  the  Birange  of  Experimen¬ 
tal  Targets 

Examples  1  and  2  above  provide  an  intuitive  demonstration  of  the  role  of  the 
bispectrum  in  identifying  multiple  interactions.  The  interpretation  of  the 
birange  is  also  valid  when  scattering  from  a  target  incorporates  unmodeled 
dynamics  such  as  frequency  dependent  amplitudes  and  dispersive  behavior. 
In  the  following,  examples  of  the  birange  of  experimental  radar  data  are 
given.  The  purpose  of  these  examples  is  to  examine  the  birange  of  real 
radar  signatures  in  terms  of  the  interpretations  given  above. 

This  section  considers  some  real  radar  targets  whose  frequency  response 
is  either  measured  experimentally  or  derived  theoretically.  The  reasons  for 
considering  the  following  examples  are:  1)  To  determine  what  can  be  ex¬ 
tracted  from  the  birange  profile  of  canonical  targets  with  known  scattering 
behavior.  2)  To  study  the  effect  of  various  parts  of  a  generic  aircraft  on 
the  birange.  3)  To  identify  some  of  the  scattering  mechanisms  from  com¬ 
plex  objects  (such  as  commercial  aircraft).  4)  To  consider  targets  whose 
backscattered  data  are  measured  in  the  S-band  where  the  term  Apq  in  the 
scattering  model  of  (3.15)  has  low  magnitude. 

The  birange  profiles  that  are  displayed  in  the  following  examples  only 
show  those  responses  with  bicoherences  >  40%.  The  reason  for  displaying 
the  birange  profiles  as  such  is  to  simplify  the  interpretation  of  the  responses 
in  terms  of  the  target  substructures.  The  triple  correlations  of  the  scatter¬ 
ing  data  used  to  compute  the  birange  profiles  are  all  windowed  using  the 
Minimum  Bias  Supremum  window  function  given  in  equation  (2.45). 
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The  following  suggestions  can  be  followed  when  examining  the  birange 
of  real  radar  targets. 

•  Ignore  any  responses  outside  the  prin^-?  2  triangle  of  each  sector. 
The  principle  triangle  of  the  first  sector  is  defined  by  ri,r2  >  0  and 
r\  +  r2  <  2^7  where  A /  is  the  sampling  period  of  the  measured  fre¬ 
quency  response.  For  properly  sampled  radar  signals  (i.e.,  satisfying 
the  Nyquisit  sampling  criterion)  no  responses  should  appear  outside 
the  principle  triangle. 

•  It  helps  to  display  the  target  along  each  axis  and  rotated  according 
to  its  orientation  angle  with  respect  to  the  radar. 

•  A  peak  at  (ra,rb)  in  the  birange  profile  is  an  indication  of  interaction 
between  a  scattering  point  on  the  target  at  a  distance  ra  from  the 
origin  and  another  scattering  point  on  the  same  target  at  a  distance 
tf,  from  the  origin. 

3.5.1  The  Birange  of  a  Canonical  Scatterer 

Canonical  scattering  objects  are  usually  considered  for  mechanism  iden¬ 
tification  problems,  and  are  useful  for  compact  range  studies.  Usually, 
scattering  from  canonical  objects  is  examined  in  the  time  domain  using  the 
transient  response.  Because  canonical  scatterers  have  a  small  number  of 
scattering  terms,  their  scattering  responses  are  relatively  simple  and  can 
be  studied  in  detail. 

The  Tilted  Flat  Plate 

The  canonical  target  considered  is  a  2ft  x  2fv  flat  plate  whose  trailing  edge 
is  tilted  by  46°  above  the  horizontal.  In  addition  to  specular  scattering  from 
the  edges  of  the  plate,  a  few  multiple  interactions  are  also  displayed  in  the 
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transient  response  (see  Figure  3.8).  A  possible  scenario  of  an  interaction 
of  this  flat  plate  is  (see  Figure  3.7)  radar,  to  first  corner,  to  second  corner, 
to  radar  where  the  first  and  second  corners  are  either  on  a  common  edge 
(denoted  by  g )  or  diagonally  opposite  (denoted  by  c).  Another  possible 
scenario  of  interactions  of  this  plate  is  radar,  to  first  edge,  to  second  edge, 
to  first  edge,  to  second  edge,  to  radar  (denoted  by  d).  Notice  that  the 
time  delay  of  the  interactions  that  follows  scenario  g  will  be  dependent 
on  the  width  of  the  plate  and  not  on  the  delay  to  each  of  the  interacting 
scatterers  because  the  signal  that  encounters  the  interaction  travels  once 
along  the  plate.  However,  the  time  delay  of  the  interactions  that  follow 
the  second  scenario  is  dependent  on  the  propagation  time  to  each  of  the 
interacting  scatterers.  Therefore,  it  is  expected  that  the  birange  will  detect 
the  interactions  that  follow  the  second  scenario. 

The  transient  response  and  the  birange  profile  of  the  flat  plate  at  5° 
azimuth  are  shown  in  Figure  3.8.  The  frequency  band  used  is  2  —  18  GHz 
with  an  increment  of  20  GHz.  Notice  that  in  both  figures  only  one  interac¬ 
tion  is  detected  in  the  birange  which  can  be  interpreted,  using  the  second 
scenario,  as  a  double  bounce  between  the  leading  and  trailing  edges  of  the 
plate  (denoted  by  d  in  Figure  3.7).  This  interpretation  can  be  further  jus¬ 
tified  by  the  fact  that  the  birange  profile  of  the  plate  at  0°  azimuth  shows 
the  same  interaction. 


3.5.2  The  Birange  of  a  Generic  Aircraft  Model 


We  first  examine  the  scattering  from  a  generic  aircraft  model  with  remov¬ 
able  parts  (denoted  by  FWST  (75]  in  Figure  3.9).  The  advantages  of  study¬ 
ing  a  target  with  removable  parts  are:  1)  to  examine  the  contribution  of 
each  part  of  the  aircraft  to  the  total  number  of  scattering  mechanisms  that 
can  be  extracted,  2)  to  follow  the  changes  in  the  transient  responses  or  the 
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Figure  3.8:  Transient  response  and  birange  profile  for  tilted  flat  pla.te  at  5° 
azimuth. 
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birange  profiles  as  more  paxts  are  added  to  (or  removed  from)  the  target. 

The  frequency  band  of  the  data  used  in  this  example  is  2  GHz  to  18  GHz 
with  a  frequency  increment  of  100  MHz.  The  FWST  model  is  6  inches  long 
with  a  wing  span  of  4  inches.  Figures  3.10-3.12  show  the  birange  profile 
of  different  configurations  of  the  aircraft  parts.  Notice  that  for  the  case 
where  scattering  is  from  the  fuselage  only  (see  Figure  3.10),  there  is  a 
single  response  in  the  birange  which  could  be  due  to  scattering  from  the 
nose  and  then  the  trailing  edge  of  the  fuselage. 

If  the  horizontal  stabilizer  is  added  to  the  fuselage  (see  Figure  3.11), 
then  a  multiple  interaction  between  the  leading  edge  of  the  stabilizer  and 
the  nose  (denoted  by  a  in  Figure  3.11)  is  detected.  Also  detected  in  the 
birange,  is  an  interaction  (denoted  by  6  in  Figure  3.11)  between  the  leading 
edge  of  the  stabilizer  and  the  opening  in  the  fuselage  where  the  wing  part 
of  the  aircraft  is  inserted.  Similarly  an  interaction  (denoted  by  c  in  the 
same  Figure)  is  detected  between  the  trailing  edge  of  the  stabilizer  and 
the  same  opening  in  the  fuselage.  The  responses  on  the  diagonal  are  not 
well  understood  but  may  be  due  to  scattering  from  the  corners  between  the 
front  edge  of  stabilizer  and  the  fuselage. 

Figure  3.12  shows  the  birange  profile  and  the  transient  response  when 
all  parts  of  the  aircraft  are  present.  Notice  that  the  response  denoted  by  a  is 
an  interaction  between  the  trailing  edge  of  the  wing  and  the  corner  between 
the  wing  leading  edge  and  the  fuselage.  A  similar  interaction  (denoted  by 
b)  between  the  leading  edge  of  the  stabilizer  and  the  wing-fuselage  corner 
is  also  detected.  The  response  on  the  diagonal  (denoted  by  c)  is  possibly 
due  to  interaction  between  both  ends  of  the  wing  (corresponds  to  scattering 
from  one  end  of  the  wing  then  scattering  from  the  other  end).  Notice  that 
the  interactions  between  the  stabilizer  and  the  opening  in  the  fuselage  that 
are  shown  in  Figure  3.12  (denoted  by  6  and  c)  do  not  appear  in  Figure  3.12 
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after  inserting  the  wing  into  the  fuselage. 

The  above  examples  of  the  birange  profiles  of  the  generic  FWST  aircraft 
show  that  interactions  between  the  edges  of  the  stabilizer  and  the  edges  of 
the  wing  are  likely  to  appear  in  the  birange.  Further,  the  addition  of  parts 
of  large  size  relative  to  size  of  the  fuselage  (such  as  the  stabilizer  or  the 
wing)  yields  significant  changes  in  the  birange  profile  of  the  generic  target. 

3.5.3  The  Birange  of  Real  Aircraft  Models 

In  the  following,  we  study  the  birange  profiles  of  real  aircraft  models.  The 
birange  profiles  of  the  generic  aircraft  model  studied  in  the  previous  section 
may  help  examining  the  birange  of  real  aircraft  models.  The  major  differ¬ 
ence  between  scattering  from  the  generic  aircraft  and  a  real  target  model 
is  scattering  from  the  engines.  The  engines  are  often  major  contributors  to 
scattering  from  real  radar  targets  [75].  Therefore,  in  addition  to  scattering 
from  the  wing  and  stabilizers,  interactions  between  the  engines  and  other 
scattering  subcomponents  of  a  real  target  are  likely  to  be  detected  in  the 
birange  profile. 

There  are  two  reasons  for  examining  the  birange  profile  of  real  aircraft 
models:  1)  such  targets  are  likely  radar  targets,  and  2)  their  structures  are 
very  complex  and  therefore  consist  of  a  broad  spectrum  of  scattering  mech¬ 
anisms  including  multiple  interactions.  The  goals  of  considering  the  birange 
of  real  aircraft  models  are:  1)  to  attempt  to  identify  some  of  the  scattering 
mechanisms  in  the  birange,  and  2)  to  discriminate  between  targets  based  on 
their  birange  profiles.  The  discrimination  problem  is  considered  in  Chapter 
V. 

The  birange  profiles  for  scaie  models  of  five  commercial  aircraft  models 
are  shown  in  Figures  3.13-3.17.  Details  of  the  radar  measurements  of  these 
targets  can  be  found  in  [68] .  The  backscattered  signals  are  in  the  frequency 
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Figure  3.9:  Photo  of  commercial  aircraft  used  in  this  study. 
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Figure  3.12:  Impulse  response  and  birange  profile  for  generic  aircraft  model 
(complete  set). 


range  (1.5  —  12)  GHz  which  corresponds  to  the  resonance  scattering  region 
for  full  scale  targets.  The  polarization  mode  used  is  Horizontal  transmit 
Horizontal  receive  (HH).  Notice  that  the  dominant  responses  in  each  figure 
correspond  to  scattering  from  the  wing  and  engines.  Additional  responses 
around  the  stabilizer-tail  region  are  also  recovered. 

Figure  3.13  shows  the  birange  profile  and  the  transient  response  for  the 
Boeing  707  at  0°.  The  birange  of  this  target  shows  interactions  between  the 
engines  and  the  stabilizer  (denoted  by  a),  between  both  edges  of  the  stabi¬ 
lizer  (denoted  by  b),  and  between  the  fuselage-tail  corner  and  the  leading 
edge  of  the  wing  ^denoted  by  c).  Figure  3.13  also  shows  other  interactions 
of  lower  magnitudes.  Figure  3.14  shows  the  transient  response  and  the 
birange  profile  for  the  Boeing  727.  The  interactions  detected  in  this  case 
are  between  the  engines  and  the  wing-fuselage  corner  (denoted  by  a)  and 
between  the  wing  tip  and  the  front  edge  of  the  fuselage  (denoted  by  b). 
Some  of  the  responses,  however,  cannot  be  explained  easily  (such  as  the 
one  denoted  by  c)  and  are  perhaps  due  to  interactions  between  the  win¬ 
dows  of  the  airplane  and  other  parts  of  the  aircraft.  Figure  3.14  also  shows 
other  interactions  of  lower  magnitudes.  Figure  3.15  shows  the  interactions 
detected  in  the  birange  profile  of  the  Concord.  These  interactions  are  be¬ 
tween  the  engines  and  fuselage  front  edge  (denoted  by  a),  between  the  tail 
and  the  wing-fuselage  corner  (denoted  by  6),  and  between  the  tail  and  the 
wing  tip  (denoted  by  c).  The  birange  profile  of  the  Boeing  747  (see  Fig¬ 
ure  3.16)  shows  interactions  between  the  wing  tips  and  the  engines  (denoted 
by  a)  and  between  the  stabilizer  and  the  wing-fuselage  corner  (denoted  by 
b).  Figure  3.17  shows  the  birange  and  the  impulse  response  for  the  DC10 
target  model.  Three  interactions  are  detected  in  the  birange  of  this  target; 
the  first  is  between  the  engine  and  front  edge  of  the  fuselage  (denoted  by 
a),  the  second  is  between  the  stabilizer  and  wing-fuselage  corner  (denoted 
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by  b),  and  the  third  is  perhaps  between  a  window  of  the  airplane  and  the 
leading  edge  of  the  wing  (denoted  by  c). 

It  is  clear  from  the  above  figures  that  some  aircraft  models  have  more 
significant  interactions  than  other  models.  Hence,  the  birange  is  more  infor¬ 
mative  for  some  targets  than  others.  Also,  it  is  clear  from  these  figures  that 
different  targets  have  significantly  different  birange  signatures.  Therefore, 
the  birange  profile  is  a  potential  candidate  for  target  classification. 

After  detecting  the  interaction  terms  (with  bi coherence  «  1)  in  the 
birange  profiles  of  the  above  target  models,  the  corresponding  terms  in  the 
transient  responses  were  located.  By  examining  the  relative  amplitude  of 
the  interaction  terms  of  real  target  models  as  they  appear  in  the  transient 
response,  it  is  possible  to  conclude  that  these  interactions  have  relatively 


low  magnitude  compared  to  the  magnitude  of  the  responses  that  correspond 
to  specular  scattering.  To  quantify  this  result  let  SIR  denote  the  ratio  of 
specular  scattering  to  interactions  in  dB  as  they  appear  in  the  transient 
response.  The  specular  scattering  component  is  computed  as  the  average 
of  the  magnitudes  of  the  interacting  scatterers  (|h(ri)|  +  |/i(r2)|)/2.  The 
magnitude  of  the  interaction  is  given  as  |h(7,i-fr2)|  provided  that  no  specular 


scattering  appears  at  rj  +  r2  (bicoherence  %  1).  SIR  is  then  given  as 

^=2»'-[JtFrS!]  (331) 

After  estimating  SIR  for  several  interaction  terms  in  the  above  examples, 


we  observe  that 


1  dB  <  SIR  <  20  dB. 


(3.32) 


This  result  indicates  that  the  interaction  terms  are  not  always  easily  rec¬ 
ognized  in  the  transient  response  of  real  target  models.  Therefore,  it  is 


reasonable  to  assume  that  most 


responses  of  comparable  magnitudes  in  the 


transient  response  of  a  real  target  model  correspond  to  specular  scatter¬ 


ing.  If  the  goal,  however,  is  to  search  for  interactions  and  the  interest  is 
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in  the  type  of  scatterers  that  produced  these  interactions,  then  the  above 
result  shows  the  importance  of  the  birange  profile  in  detecting  multiple 
interactions,  especially  those  with  low  magnitudes. 

3.5.4  The  Birange  of  S-Band  Scattering  Data 

In  this  section  we  examine  the  birange  profile  of  an  aircraft  measured  with 
an  S-band  radar.  The  target  size  is  about  200  times  the  operating  wave¬ 
length.  The  measurement  frequency  band  is  3.1  —  3.4  GHz  for  a  full  scale 
aircraft.  The  data  consist  of  64  frequency  measurements  corresponding  to 
a  frequency  increment  of  4.1  MHz.  The  measurements  were  taken  while  the 
aircraft  was  in  operation  and  are  therefore  corrupted  with  noise.  Details 
on  the  calibration  and  motion  compensation  of  thL  data  are  given  in  [69]. 

Scattering  at  such  high  frequencies  is  very  localized  and  highly  spec¬ 
ular  [75].  In  general,  due  to  the  large  target-size  to  wavelength  ratio, 
the  specular  part  of  the  signal  is  the  most  dominant  in  the  transient  re¬ 
sponse  [75]. 

Figure  3.18  shows  the  transient  response  and  the  birange  of  the  A4 
whose  backscatter  data  were  measured  at  S-band.  It  is  clear  that  the  bi¬ 
range  of  such  a  target  reveals  the  complex  scattering  mechanisms  between 
various  subcomponents  of  the  target  scattering  elements.  Several  inter¬ 
actions  are  detected  particularly  between  the  engines  and  the  stabilizer 
(denoted  by  a),  between  the  tail  and  the  wing  tip  denoted  by  6),  and  be¬ 
tween  the  fuel  tank  and  the  wing  (denoted  by  c).  A  common  feature  of 
all  detected  interactions  is  that  the  SIR  is  higher  than  that  obtained  for 
scattering  in  the  resonance  region.  The  reason  is  that,  while  scattering  in 
the  S-band  is  more  localized  and  can  be  more  easily  resolved,  compare  d 
to  scattering  in  the  resonance  region,  the  propagating  short  wave  between 
interacting  scatterers  becomes  weaker  due  to  diminished  electromagnetic 
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Figure  3.15:  Birange  profile  for  the  Boeing  747  model  at  0°  azimuth 
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coupling  [71]. 

It  is  possible  to  extract  many  details  about  the  target  scattering  behav¬ 
ior  in  the  S-band  frequency  range  from  the  birange  profile.  These  interac¬ 
tions,  however,  are  sensitive  to  changes  in  the  target  aspect  angle  because 
of  the  large  size  of  the  target  relative  to  the  wavelength. 

3.6  Extraction  of  Scattering  Features  From 
The  Transient  Response  and  The  Birange 

The  above  examples  of  the  birange  of  real  aircraft  models  have  shown  that 
some  of  the  interactions  can  be  detected  using  the  birange  concept.  The 
coordinates  of  the  responses  in  the  birange  indicate  the  location  of  the 
interacting  scattering  subcomponents.  The  accuracy  of  these  coordinates 
depends  on  the  fast  Fourier  transform  algorithms  and  the  number  of  pixels 
used  to  display  the  birange  profile.  It  is  also  clear  from  these  figures  that 
it  is  possible  to  characterize  a  radar  target  by  a  small  number  of  features 
extracted  from  both  the  transient  response  and  birange  profile.  Therefore, 
it  is  possible  to  reduce  the  data  storage  requirement  of  the  birange  pro¬ 
file  if  the  interactions  shown  in  these  examples  are  extracted  and  labeled 
in  terms  of  their  coordinates  and  bicoherences.  An  algorithm  for  extract¬ 
ing  scattering  mechanisms  from  the  birange  profile  when  estimated  using 
classical  bispectral  estimation  techniques  is  given  in  this  section. 

It  is  possible  to  use  the  extracted  features  in  a  structural  pattern  recog¬ 
nition  machine  that  discriminates  between  targets  using  relational  graph 
matching  techniques  [76].  Also,  estimates  of  the  locations  of  scattering 
mechanisms  (which  are  more  accurate  than  the  estimates  given  by  the 
FFT  algorithm)  can  be  obtained  using  the  following  feature  extraction  al¬ 
gorithm.  The  following  algorithm  is  a  generalization  of  the  peak  search 
algorithm  used  in  [77]  for  estimation  the  location  of  the  peak  in  the  tran- 
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Figure  3.18:  Transient  response  and  birange  profile  of  the  A4  aircraft 
(S-band  scattering). 
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sient  response. 


Feature  Extraction  Algorithm: 

1-  Estimate  the  birange  profile  and  the  transient  response  as  given  in  Section 
3.1.  (Better  results  were  obtained  when  the  window  function  used  is  Gaus¬ 
sian  [77],  possibly  because  the  Fourier  transform  of  a  Gaussian  function  is 
also  Gaussian). 

2-  Search  for  local  maxima  in  the  birange  where  each  selected  peak  occurs 
among  six  adjacent  points  as  shown  in  Figure  3.19.  If  (tu^2j)  indicates  the 
location  of  the  selected  peak,  then 

|-B(tii,t2i)|  “  max{|i?(ti(i_i), t2(j— 1))|,  ^2(j+i))|>  (3.33) 

W«,ia(/-i))|>  l-®(*i*>*2(j+i))|>  |J9(<i«<2i)| 

|#(tl(i+l)>t2(j-l))|>  |-#(*l(t+l)>*2j)|}- 

Also,  search  for  the  peaks  in  the  transient  response  using  the  algorithm 
given  in  [75]. 

3-  Fit  a  Gaussian  function  to  each  set  of  data  consisting  of  the  selected 
peak  in  the  birange  and  the  neighboring  six  minima.  Let  <lc  and  t2c  denote 
the  coordinates  of  a  local  maximum.  Then  the  Gaussian  window  function 
is  given  as 

w(tut2)  =  Aexp{-a(ti  -  tlc)2  - /3(t2  -t2c)2}  (3.34) 

Given  the  coordinates  and  the  magnitudes  of  seven  bispectral  samples  (see 
Figure  3.20)  including  the  maximum  given  using  the  FFT  algorithm,  we 
need  to  solve  for  a,  j3 ,  t\c,  t2 c,  and  A.  Therefore,  estimates  of  these  pa¬ 
rameters  can  be  obtained  using  local  least-squares  fit  on  all  seven  samples. 
The  estimates  of  tlc  and  Uc  jn.  the  bispectral  response  can  then  be  obtained 
by  the  following 

u.  =  ^4^  <3-36> 
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he  — 


Xs  +  $4 


where 


_ (ri  ~  ~  ^i(t-i)) _ 

2ti(^i(*— i)  —  tu)  —  2(<i(i_i)  —  hi) 

{r2  -  l)(t2u  ~  tlg-p) 

2ri(h(i-i)  -  hi)  ~  2(<i(i_i)  —  hi) 
_ (r3  ~^){t\i~  ^l(j-i)) _ 

2r2(i2(j-i)  —  hj)  ~  2(t2(j-i)  —  t2j) 

_ (r4  ~  l)(*ii  ~  *l(i-l)) _ 

2ra{h(i-i)  ~  hi)  ~  2(<i(*_i)  —  hi) 


(3.36) 


ln(|^(^it,  ^2j)|)  —  ln(l^(<i(t+1),/2j-)|) 

ln(\B(hi,  t2j)\)  -  ln(|£(/1(;_i),  <2j)|) 
ln(\B(huhj)\)  -ln(l^(<1f^2(j+i))|) 
ln( \B(tu ,i2,-)|)  “  ln(|5(<li,<2(i-t))|) 
ln(|B(<i(,--i),<afj+i))|)  -  ln(l^(<i(»-x),<2j)l) 

ln(|.B(<i(i+i),<2j)|)  _  ln(|P(<i(,-+i),i2(j-i))|) 

_  ln(|-P(<i(»-i), <2(j+i))|)  -  ln(|i?(<it,f2(j+i))|) 

H\B(h(i+i),h(j-i))\)  ~  ln(|5(<n,<2(i_1))|) 


(3.37) 

(3.38) 

(3.39) 

(3.40) 


4-  Fit  a  Gaussian  window  to  each  set  of  data  consisting  of  the  selected  peak 
and  the  neighboring  minima  in  the  transient  response  [77].  The  location  tc 
of  the  peak  in  the  transient  response  can  then  be  estimated  using 


=  *ln  -*g  +  r(*fc  +  *Li) 

2  r(tk  —  tk-i)  +  2<k+i  +  2  tk 


(3.41) 


where  tk  and  tk+i,tk-i  denote  the  location  of  the  selected  peak  as  given 
using  the  Fourier  transform  and  the  location  of  the  adjacent  values.  The 
parameter  r  is  given  as 


ln(|/t(fe)|)  —  ln(|h(fc  +  !)|) 
ln(|/i(fc)|)  —  ln(|h(fc  —  1)|) 


(3.42) 
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where  \h(k)\  is  the  magnitude  of  the  transient  response  at  time  k. 

5- All  detected  maxima  must  satisfy  the  symmetry  properties  of  the  birange. 
This  condition  is  needed  to  detect  any  possible  failures  of  this  peak  search 
algorithm.  If  some  peaks  do  not  satisfy  the  symmetry  properties  of  the 
bispectrum  then  repeat  the  search  starting  from  a  different  corner  of  the 
bispectral  hexagon. 

6- If  the  measured  data  are  noisy  then,  to  eliminate  noise  peaks,  ignore  any 
local  maximum  in  the  birange  profile  where  the  magnitude  of  the  birange  at 
this  maximum  is  such  that  |B(r!,r2)|  <  K<r3  where  cr  is  the  noise  standard 
deviation  and  if  is  a  constant  (if  =  1  is  used  in  this  study).  Also,  ignore 
any  local  maximum  in  the  transient  response  where  the  amplitude  of  the 
response  at  this  maximum  is  such  that  |h(r)|  <  a 

Experimental  Results 

Table  3.1  shows  two  examples  of  the  estimates  of  the  coordinates  of  the 
interaction  term  of  the  blade-sphere  target  using  the  above  algorithm.  The 
true  values  of  the  coordinates  and  those  given  by  the  Fourier  transform 
algorithm  are  also  shown.  The  range  resolution  is  6r  =  0.015  meters,  cor¬ 
responding  to  10  GHz  frequency  band.  Notice  that  the  estimates  obtained 
using  the  above  algorithm  are  closer  to  the  true  values  compared  to  those 
given  by  the  fast  Fourier  transform. 

The  features  extracted  from  the  birange  profiles  of  five  commercial  air¬ 
craft  using  the  above  algorithms  are  listed  in  Tables  3.2-3. 3.  Notice  that 
neither  the  number  of  scattering  centers  (extracted  from  the  transient  re¬ 
sponse)  nor  the  number  of  interactions  (extracted  from  the  birange)  is  the 
same  for  both  targets. 


Figure  3.19:  Feature  extraction  technique;  locating  a  maximum  surrounded 
by  six  minima 
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Table  3.1:  Comparison  between  bi-range  estimates  obtained  using  Gaussian 
fitting  technique,  FFT  algorithm  and  the  true  values. 


I 

I 

I 


Blade-Sphere  Target  Sr  =  .015  m 

Location  of  Interaction  Term 

Examples 

True 

FFT 

FIT 

Example  #1, 

0.1789 

0.1715 

0.1745 

Example  #1,  r2 

0.3111 

0.3185 

0.3168 

Example  #2,  ri 

0.0872 

0.0725 

0.0726 

Example  #2,  r2 

0.3965 

0.3869 

0.4105 

Table  3.2:  Features  extracted  from  birange  profile  and  transient  response 
for  a  Boeing  707 


B707  at  SNR  =  10  dB,  0  deg.  Azimuth 

Model  length  «  .32  m,  Sr  «  .015  m 

Scatt.  Centers 

Range  [meters] 

Amplitude  [dB] 

No.  1 

.0165 

1.15 

No.  2 

.1213 

2.14 

No.  3 

.1901 

2.36 

No.  4 

.2372 

2.03 

Interactions 

Birange  (ra,r2) 

Bicoherence 

No.  1 

(.1209,-2822) 

47  % 

No.  2 

(.1933, .3185) 

53  % 
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Table  3.3:  Features  extracted  from  birange  profile  and  transient  response 
for  a  Boe\ag  727 


B72’7  at  SNR  =  10  dB,  0  deg.  Azimuth 

Model  length  «  .235  m,  6r  ~  .015  m 

Scatt.  Centers 

Range  [meters] 

Amplitude  [dB] 

No.  1 

.0143 

1.13 

No.  2 

.1036 

1.69 

No.  3 

.1775 

1.94 

No.  4 

.2249 

0.90 

Interactions 

Birange  (r1}r2) 

Bicoherence 

No.  1 

(.0673,-1106) 

87  % 

No.  2 

(.0252,1765) 

44  % 

No.  3 

(.1841,1093) 

52  % 

3.7  The  Effect  of  Noise  on  The  Birange  Pro¬ 
file 

Radar  scattering  signatures  are  usually  corrupted  with  additive  noise  and 
clutter.  In  this  section  we  investigate  the  effect  on  the  birange  profile  of 
additive  disturbances  to  the  frequency  domain  scattering  data.  The  distur¬ 
bance  model  used  for  noise  and  clutter  is  a  stochastic  model  [45].  The  noise 
model  is  Gaussian  (colored  or  white).  Ksienski  and  White  [22],  described 
the  noise  model  of  a  radar  system  as  being  both  additive  and  multiplicative, 
where  the  additive  part  is  due  to  thermal  and  shot  noise,  and  the  multi¬ 
plicative  part  is  due  to  measurement  errors  and  is  dependent  on  the  signal. 
According  to  [22],  if  none  of  the  noise  sources  is  dominant  and  by  virtue 
of  the  central  limit  theorem,  the  assumed  noise  model  is  zero  mean  white 
Gaussian.  The  clutter  model  is  either  Gaussian  (due  to  migrating  birds  and 
weather,  [17])  or  non-Gaussian  (due  to  ground  clutter  [17]).  Combining  the 
two  models,  it  is  possible  to  classify  disturbances  into  two  groups:  Gaussian 
and  non-Gaussian. 
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3.7.1  Additive  Gaussian  Noise 

The  third  moment  of  any  zero  mean  random  variable  with  symmetric  dis¬ 
tribution  is  zero.  Hence,  the  triple  correlation  of  a  zero  mean  random 
variable  with  symmetric  probability  density  function  is  zero,  which  implies 
that  the  bispectrum  of  a  pure  noise  signal  is  zero.  Let  H'(f)  denote  the 
noise  contaminated  scattered  radar  signal  at  frequency  /. 

II'(f)  =  H(f)  +  n(f)  (3-43) 

where  n(/)  is  a  zero  mean  Gaussian  random  variable,  and  H(f)  is  the  noise¬ 
less  backscattered  signal.  Also,  assume  that  the  sample  mean  is  removed 
from  both  the  noisy  data  H'(f )  and  the  noiseless  catalog  data  H(f),  so 
that  E{H'(f)}  =  0.  Then  the  third  order  cumulant  of  H'(f)  is 

R"‘(fuh)  =  f  +  +  /2)}  =  RH(fi,fi)  (3.44) 

=>  BH'(tuh)  =  BH(tut2). 

Therefore,  the  birange  of  noisy  radar  signals  should  be  unaffected  by  ad¬ 
ditive  Gaussian  noise  provided  that  an  infinite  number  of  data  samples  is 
available.  In  practice,  however,  only  a  finite  number  of  data  samples  is 
available.  Hence,  additive  Gaussian  noise  is  not  totally  suppressed  in  the 
birange  profile. 

Figure  3.21  shows  the  transient  response  and  the  birange  profile  of  the 
hypothetical  target  of  Example  1  where  the  simulated  radar  returns  are 
corrupted  with  additive  white  Gaussian  noise  at  -10  dB  signal-to-noise  ra¬ 
tio  (SNR).  Compare  Figure  3.21  with  the  corresponding  noiseless  profiles 
shown  in  Figure  3.2.  Notice  that,  while  it  is  difficult  to  identify  some  of 
the  scattering  mechanisms  in  the  noisy  transient  response,  it  is  possible  to 
do  so  in  the  noisy  birange  profile  as  shown  in  Figure  3.21.  Thus,  birange 
processing  seems  to  more  effectively  suppress  the  additive  Gaussian  noise 
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than  does  transient  response  processing.  This  does  not  necessarily  imply 
that  noisy  radar  targets  can  be  classified  with  a  lower  probability  of  error 
using  the  birange  (see  Chapter  V). 

Suppression  of  additive  Gaussian  noise  is  dependent  on  the  number 
of  data  samples  used  to  estimate  the  triple  correlation  of  the  data.  To 
study  the  effect  of  noise  on  the  birange  profile  when  it  is  estimated  using 
either  finite  or  infinite  data  sequences,  it  is  appropriate  to  examine  the 
estimates  of  the  triple  correlations  of  noisy  radar  returns.  The  expected 
value  of  the  triple  correlation  of  noisy  backscatter  data  denoted  by 

RH\hif2 )  and  estimated  using  (3.7)  is  related  to  expected  value  of  the 
triple  correlation  of  the  noiseless  data  RH(fi,f2)  by  the  following; 


RH\h,h)  =  R!IUuh)  (3.45) 

+4  EW'  +  AHA  +  A)  +  £  i>‘(AWA  +  AHA  +  A) 


1  N 


»=1 

N 


+ Jf  T,m  +  AW  A  +  AK(A)  +  Tf  £jr(A)»(A.+  AHA  +  A) 


i'=l 


N  ■ 

i=i 


+4  2>(  A  +  AH(AMA  +  A)  +  4  i>(A  +  AK(AMA  +  A) 

iv  t=l  iv  »=1 

+  4  XX(/i)n(/i  +  /l)»( fi  +  f*)- 

•/V  i=l 

As  IV  — >  oo  the  second,  third,  and  fourth  terms  on  the  right  hand  side  of  the 
above  equation  converge  to  [57]  {rH(fi)  +  rH(f2 )  +  rH(fi  -  f2)}E{n}  =  0 
where  rH(fk )  =  E{H*(f)H(f  +  /*)},  and  E{n}  =  0.  Also,  the  fifth,  sixth, 
and  seventh  terms  converge  to  Z(t2E{H }  =  0,  where  <t 2  is  the  noise  variance 
and  E{H}  =  0  (due  to  subtracting  the  mean  from  the  data).  The  last  term 
can  be  written  as  E{n3}  =  0.  Segmentation  or  averaging  over  independent 
records  may  be  needed  for  convergence  of  these  terms  [57].  Therefore,  as 
N  — »  oo  the  estimate  of  the  noisy  triple  correlation  approaches  the  estimate 
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of  the  noiseless  triple  correlation,  thus  producing  noise-free  birange  profiles. 


lim  RH'(fi,h)-*Ra(h,h)  (3.46) 

N—*oo 


An  expression  for  the  expected  value  of  the  error  between  the  estimate 
of  the  triple  correlation  using  a  finite  sequence  of  noisy  radar  data  and  the 
estimate  using  noiseless  data  is  given  below.  This  expression  represents  the 
effect  of  noise  on  the  birange  profile  when  estimated  using  a  finite  number 
of  data  samples  N.  Let  a  =  rff(0).  Then  after  algebraic  manipulation  of 
products  of  the  terms  in  (3.45),  the  following  is  obtained 


E{RH\fuf2)  -  RH{fuh)Y  =  [9aV  +  45acr4  +  15<r°]  (3.47) 


where  cr2  is  the  noise  variance.  This  means  that  the  mean  square  error 
between  the  noisy  triple  correlation  RH'  and  the  noiseless  triple  correlation 
RH  is  dependent  on  the  noise  power  and  increases  as  a  function  of  cr6. 
Equation  (3.47)  can  be  regarded  as  a  quantitative  measure  of  the  effect  of 
noise  on  the  birange  when  estimated  using  finite  data  records. 


3.7.2  Additive  Non-Gaussian  Noise 

By  virtue  of  the  central  limit  theorem  [22],  the  additive  noise  in  the  radar 
scattering  data  are  assumed  to  be  Gaussian.  Generally,  the  additive  noise 
is  not  Gaussian  for  two  reasons:  1)  The  noise  is  partially  multiplicative, 
and  so  is  only  asymptotically  Gaussian.  2)  Clutter  is  often  modeled  as  a 
non-Gaussian  disturbance  added  to  the  scattering  data. 

In  this  section  we  study  the  effect  of  additive  non-Gaussian  noise  on  the 
birange.  Consider  the  following, 

H'U)  =  «(f)  +  0(f)  (3.48) 

where  H'(f )  is  the  backscattered  frequency  response  of  a  certain  target 
and  Gf  represents  non-Gaussian  noise.  The  triple  correlation  of  the  above 
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Figure  3.21:  Transient  response  and  birange  profile  of  the  three-point  syn¬ 
thesized  target  of  example  #  1  at  SNR=  -10  dB. 
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signal  (for  infinite  sequence  of  data)  is  given  as 

BPUuh)  =  +  h)H'(!  +  /2)} 

=  Rtt(fuf2)  +  R°(f„h). 

Thus,  the  effect  of  non-Gaussian  noise  on  the  bispectral  response  of  a  zero 
mean  signal  is  represented  by  Rc(fufa),  which  is  dependent  on  the  third 
moment  of  G.  Non-Gaussian  disturbances  (such  as  ground  clutter)  often 
have  non-zero  first  and  third  order  moments  [17].  For  this  reason,  additive 
non-Gaussian  noise  is  not  suppressed  in  the  birange  even  for  an  infinite 
number  of  samples. 

For  a  finite  data  sequence,  following  the  same  derivations  as  in  the 
Gaussian  noise  case,  the  effect  of  additive  non-Gaussian  noise  on  the  triple 
correlation  can  be  quantified  as 

B{R^'Uuf2)-RB{fuh)y-  (3.49) 

i  \!)a2B{C}2  +  9 aE{0“}  +  6a£{C}£{C3}  +  W£{C6}] 

where  a  =  rff(0).  Notice  that  if  N  — >  oo  the  mean  square  error  between 
noisy  and  noiseless  triple  correlations  approaches  E{C6}.  Also,  because 
E{C}  ^  0,  and  JE?{C73}  ^  0,  the  term  6aE{C}E{C3}  contributes  to  the 
mean  square  error.  However,  this  result  does  not  imply  that  classification 
of  radar  targets  in  noisy  environment  (where  the  noise  is  non-Gaussian) 
based  on  features  extracted  from  the  birange  produces  poor  performance 
compared  to  using  other  target  representations.  The  reason  is  that  non- 
Gaussian  noise  has  similar  effects  on  other  types  of  display  of  target  fea¬ 
tures. 


3.8  Summary 

The  birange  profile  of  radar  signatures  was  introduced  in  this  chapter. 
The  definition  of  the  birange  as  a  display  of  target  signatures  in  the  two- 
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dimensional  time  domain  is  based  on  the  bispectrum.  The  introduction  of 
the  birange  to  radar  signature  analysis  is  motivated  by  the  interest  in  iden¬ 
tifying  multiple  interactions,  thus  achieving  a  more  accurate  description  of 
the  target  scattering  subcomponents  and  enhancing  the  target  identifica¬ 
tion  process.  The  birange  of  radar  signatures,  as  defined  in  this  chapter,  is 
based  on  a  scattering  model  platform  that  includes  specular  scattering  and 
multiple  interactions.  The  model  used  is  a  sum  of  complex  exponentials  at 
frequencies  related  to  the  locations  of  scattering  centers  along  the  target. 

The  estimation  of  the  birange  is  slightly  different  from  conventional 
bispectral  estimation  techniques.  Details  of  an  estimation  algorithm  are 
given  in  this  chapter.  Also,  the  interpretation  of  the  responses  in  the  birange 
is  presented  using  two  examples  of  synthesized  targets.  Means  for  detecting 
false  interactions  in  the  birange  profile  are  also  discussed.  Further,  the 
effect  of  unmodeled  dyn?mics  in  the  backscattered  signal  on  the  birange  is 
investigated. 

The  application  of  the  birange  profile  as  a  radar  signature  analysis  tech¬ 
nique  is  first  examined  using  experimental  radar  returns  for  a  canonical  ob¬ 
ject  with  somewhat  known  scattering  behavior.  Then  the  birange  profiles 
of  experimental  radar  returns  from  a  generic  aircraft  model  (with  remov¬ 
able  parts)  are  investigated.  The  detected  responses  in  the  birange  profile 
of  several  configurations  of  the  parts  of  this  aircraft  are  associated  with  the 
geometry  of  the  aircraft. 

The  birange  profiles  of  experimental  radar  returns  from  real  aircraft 
models  at  different  frequency  bands  are  also  examined  in  terms  of  the  de¬ 
tected  interactions  and  the  geometry  of  the  targets.  Also,  a  feature  extrac¬ 
tion  algorithm  to  estimate  the  location  of  the  interactions  in  the  birange 
profile  is  given.  This  algorithm  is  used  to  provide  tabulated  features  ex¬ 
tracted  from  both  the  transient  response  and  the  birange.  The  extracted 
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features  can  be  used  to  identify  radar  targets  using  structural  pattern  recog¬ 
nition  algorithms  without  requiring  large  memory  space.  Finally,  the  effect 
of  additive  Gaussian  and  non-Gaussian  noise  on  the  birange,  when  esti¬ 
mated  using  either  finite  or  infinite  data  sequences  is  investigated. 


Chapter  4 

BIRANGE  ESTIMATION 
USING  AUTOREGRESSIVE 
MODELING 


4.1  Introduction 

This  chapter  is  con.  rned  with  the  estimation  of  high  resolution  birange 
profiles.  The  birange  p.  le  of  radar  signatures  as  introduced  in  the  previ¬ 
ous  chapter  is  estimated  using  classical  Fourier  transform  techniques.  The 
resolution  of  the  estimated  profiles  is  limited  by  the  measurement  band¬ 
width.  Furthermore,  the  two-dimensional  window  function  used  in  Chapter 
III  has  a  wide  mainlobe,  which  further  limits  the  resolution  in  the  birange. 
If  the  birange  is  estimated  without  windowing  the  triple  correlation,  then 
two-dimensional  sidelobes  will  appear  along  each  range  axis.  Therefore, 
with  finite  sets  of  radar  scattering  data,  resolution  problems  in  the  birange 
when  estimated  using  Fourier  techniques  are  inevitable. 

The  key  to  solving  the  resolution  problem  is  in  avoiding  the  limita¬ 
tions  imposed  by  Fourier-based  processing  of  finite  length  data  records. 
Instead,  parametric  modeling  of  the  scattering  data  is  employed.  “Para¬ 
metric”  means  that  scattering  as  a  function  of  frequency  satisfies  a  certain 
model  whose  parameters  are  estimated  from  the  measured  data.  The  main 
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advantage  of  parametrization  is  the  removal  of  the  effect  of  limited  reso¬ 
lution,  thus  providing  higher  resolution  than  can  be  obtained  with  Fourier 
transform  based  methods.  The  limitations  of  parametrization  arise  because: 
1)  electromagnetic  scattering  at  frequencies  outside  the  measurements  can¬ 
not  be  easily  predicted  from  the  measured  data,  but  depends  on  many 
factors  such  as  dispersion,  scattering  region,  etc.,  and  2)  model  parameters 
are  often  estimated  using  nonlinear  algorithms,  and  thus,  are  sensitive  to 
noise  in  the  data.  If  the  measured  backscatter  has  a  high  signal-to-noise 
ratio  and  the  target  has  a  slowly  varying  frequency  response,  then  these 
limitations  are  small,  and  significant  resolution  improvements  over  Fourier 
methods  can  be  achieved. 

Parametric  modeling  of  radar  signatures  has  been  treated  extensively 
in  the  literature  [5,7,8].  A  brief  summary  of  the  algorithms  used  and  the 
results  obtained  is  given  in  Chapter  II.  The  methods  used  are  based  on 
AR,  ARMA,  MUSIC,  and  Prony  spectral  estimation  techniques.  While  the 
ARMA  methods  ignore  any  prior  information  about  the  target  scattering 
data,  the  MUSIC  algorithm  and  the  Prony  method  incorporate  the  a  priori 
known  exponential  nature  of  the  radar  scattering  data  as  modeled  in  (3.15). 
The  consideration  of  these  parametric  methods  is  motivated  by  the  interest 
in  obtaining  high  resolution  impulse  responses. 

There  has  been  much  recent  interest  in  high  resolution  bispectral  esti¬ 
mation  [11].  The  routes  that  are  being  followed  to  achieve  high  resolution 
bispectra  can  be  classified  into  four  categories:  1)  maximum  likelihood  bis¬ 
pectral  estimation,  2)  maximum  entropy  bispectral  estimation  3)  AR,  MA, 
and  ARMA  modeling,  and  4)  harmonic  decomposition  methods  (MUSIC, 
etc..  ).  Partial  results  have  been  obtained  for  the  first  two  categories 
[11,24,34],  but  those  remain  open  problems.  It  is  not  yet  clear  whether 
high  resolution  bispectral  estimation  methods  that  are  based  on  a  maxi- 
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mum  likelihood  approach  or  on  maximizing  the  entropy  of  the  data  can 
be  found.  Unlike  spectral  estimation  methods,  bispectra  estimated  using 
AR  modeling  are  not  linked  to  either  maximum  entropy  or  maximum  like¬ 
lihood  bispectra  [11].  Hence,  it  is  not  possible  to  deduce  one  from  the 
other.  Most  of  the  recent  results  on  high  resolution  bispectral  estimation 
are  based  on  autoregressive  (AR),  moving  average  (MA),  and  autoregressive 
moving  average  (ARMA)  modeling  techniques,  and  on  harmonic  decompo¬ 
sition  methods.  In  [30],  a  method  based  on  the  MUSIC  algorithm  and 
incorporating  triple  correlations  is  proposed  for  high  resolution  bispectral 
estimation  purposes. 

The  estimation  of  the  model  parameters  is  in  some  cases  related  to  the 
order  of  polyspectra  under  investigation.  For  example,  when  the  spectral 
response  is  needed,  the  model  parameters  are  estimated  using  autocorrela¬ 
tions  or  both  triple  correlations  and  autocorrelations.  But  when  the  bis¬ 
pectral  response  is  needed  the  model  parameters  are  estimated  using  triple 
correlations  only. 

In  this  chapter,  an  algorithm  that  adopts  AR  parametric  modeling  ap¬ 
proach,  where  the  AR  parameters  can  be  estimated  using  third  order  cu- 
mulants,  is  devised  and  tested  using  experimental  radar  data.  The  goals 
are:  1)  to  attempt  to  improve  the  accuracy  of  the  extracted  multiple  inter¬ 
actions  from  the  birange  profile,  and  2)  to  obtain  encoded  birange  profiles 
of  the  targets  of  interest  in  terms  of  the  model  parameters.  The  purpose  of 
the  present  investigation  is  to  investigate  whether  AR.  modeling  can  be  used 
to  generate  high  resolution  birange  profiles.  The  quality  of  the  obtained 
birange  profiles  is  examined.  The  reason  for  discussing  AR  modeling  in 
particular  is  that  the  techniques  used  to  estimate  the  AR  parameters  de¬ 
pend  on  the  order  of  the  polyspectral  response  sought.  Incorporating  third 
order  cumulants  into  the  estimation  procedure  is  a  major  feature  of  the 
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estimation  algorithms. 


4.2  Birange  Estimation  Algorithm  Using  AR 
Modeling 


Autoregressive  modeling  of  radar  signatures  has  physical  basis  related  to 
the  nature  of  scattering  data.  Based  on  the  exponential  scattering  model 
of  (3.15),  the  measured  backscatter  at  frequency  /*,  denoted  by  H(fk ),  can 
be  modeled  as 

*(/*)“ -£«W*-i)  +  «(A)  (4-1) 

{=1 

where  e(fk)  is  zero  mean  white  noise  with  third  order  moment  (3 ,  and  p 
indicates  the  model  order  assumed  to  be  chosen  a  priori.  The  AR  system 
that  models  the  birange  is  assumed  to  be  driven  by  white  non-Gaussian 
noise.  The  noise  must  be  non-Gaussian  noise  because  otherwise  the  bispec¬ 
trum  of  the  output  is  zero.  The  parameters  {a;}  are  the  AR  parameters 
and  can  be  estimated  from  second  or  third  order  cumulants  of  the  data 
sequence  {H(fk)}. 

The  inverse  Z  transform  of  H(fk),  h(t )  is  the  target  impulse  response. 
Therefore,  the  impulse  response  of  a  radar  target  whose  backscatter  satisfies 
(4.1)  can  be  expressed  as  [8] 


h(t)  = 


1  +  £„=i  ^  exp  {— jn7r(l  -  2 t/T)} 
P 


(4.2) 


£■ 


ha  exP  {J’r(1  “  2i/T)}  -  Pi 

where  T  =  1/A/,  and  p,-,  Ci  are  the  partial  fraction  expansion  coefficients. 

The  estimation  of  the  birange  profile  using  AR  modeling  is  based  on 
the  definition  of  the  bispectrum  using  the  Cramer  spectral  representation 
(2.12).  The  birange  B(ti,t2)  is  given  as 


i?(ii,^2)  —  h(ti)h(t2)h*  (ti  +^2)  (4.3) 


111 


Using  (4.3),  the  birange  is  estimated  as 

B{tut2)  = 

_ d -  (4.4) 

Ek,l,m=o  akata^exp  {j  27r[fc(l  -  f-)  +  1(1  -  -  m(l  -  2{h^h})]} 

Hence,  the  birange  is  parameterized  by  the  AR  coefficients  {a;}.  These 
coefficients  are  to  be  estimated  from  the  cumulants  or  estimates  of  the 
cumulants.  The  above  method  for  estimating  the  birange  profile  using  an 
AR  model  has  the  following  disadvantage:  the  birange  as  given  in  (4.3)  is  a 
transformation  of  the  impulse  response  from  a  one-dimensional  profile  into 
a  two-dimensional  profile.  For  example,  a  single  peak  at  <*.  in  the  impulse 
response  is  mapped  into  three  straight  lines  in  the  birange  profile  where 
these  straight  lines  satisfy  {tj  + 12  =  <*,  <i  —  t2  =  tk,  t2  —  =  4}  because 

by  substituting  the  magnitude  of  the  impulse  response  at  these  points  in 
(4.3)  a  birange  peak  of  large  magnitude  (approaches  infinity)  is  obtained. 
A  true  peak  in  the  birange  profile,  however,  is  the  result  of  three  peaks  in 
the  impulse  response  at  U,tj,  and  4  =  U  +  tj,  and  all  three  peaks  have 
relatively  large  magnitudes. 

AR  Parameters  Estimation  Methods 

There  are  several  methods  that  can  be  used  for  the  estimation  of  the  AR 
parameters  from  the  cumulants  of  the  data  [27,28].  The  following  relation 
is  the  basis  for  estimation, 

p 

-  m,i  —  n)  =  0  V  m,n  (4.5) 

i=0 

where  the  triple  correlation  is  defined  as 

R{m,  n)  =  E{H*(k)H{k  +  m)H{k  +  n)}.  (4.6) 

The  above  relation  is  easily  derived  [42]  by  multiplying  both  sides  of  (4.1) 
by  Hm(fk  —  fm)H*(fk  —  /n)»  and  taking  the  expected  values.  The  relation 
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given  by  (4.5)  defines  an  oveniAternuned  set  of  equations  that  can  be  used 
to  solve  for  the  {a,}  coefficients- 

The  estimation  of  the  Ah,  coefficients  {oi}  is  a  prerequisite  for  estimating 
the  birange  profile.  The  nlgOtithms  used  in  this  study  for  estimating  the 
birange  (based  on  the  methods  §iven  in  [27,28,32,42])  are  the  following: 
Algorithm  #  1:  Estimate  the  triple  correlation  77(ra,n),  m,n  =  0,...,  of 
the  measured  data  Ii(fk)  'theb.  substitute  in 


77(0,0)  77(1,1) 

7?(-l,  -1)  77(0,0) 


[  77(— p,  — p)  7?(-p  +  l,-p  4-1) 
which  can  be  denoted  as, 


77  (p,p) 

'  1 

i?(p-l,p-  1) 

ai 

77(0,0) 

.  ap  . 

07o  =  O. 


0(4.7) 


(4.8) 


Notice  that  the  entries  in  th^  matrix  57  are  chosen  from  only  one  slice 
(77(m,m),  m  =  n)  of  the  chtfUalant  plane. 

Algorithm  #  2:  It  is  possible  to  choose  entries  that  belong  to  many  slices 


nulant  plane  as  gi vsh.  below 

77(0,0) 

77(1,1) 

7?(p,p) 

77(0, -1) 

77(l,0) 

...  77(p,p -  1) 

77(0, -p) 

l  >-  p) 

77(p,0) 

7?( — 1, 0) 

77(0,1) 

...  77(p-l,p) 

-R(-p,  -p) 

77(-ptl,-p-l-  1) 

...  77(0,0) 

1 

ax 

a„ 


=  (H.9) 


where  the  {oi}f=1  are  obtained  hy  solving  an  overdetermined  set  of  linear 
equations.  In  both  algorithm  the  estimates  of  the  birange  at  (ti,t2)  are 
computed  using  (4.5). 

Unlike  using  Yule- Walker  equations  to  solve  for  the  Alt  parameters 
(where  the  autocorrelation  iA^ttijc  is  always  full  rank),  the  triple  corre¬ 
lation  matrix  is  not  necessarily  full  rank,  particularly  for  algorithm  #1. 
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Estimates  of  the  birange  profile  using  AR  modeling  techniques  exist  if  a 
full  rank  triple  correlation  matrix  can  be  found.  However,  if  a  full  rank 
triple  correlation  matrix  exists  under  certain  conditions,  it  is  not  guaran¬ 
teed  to  remain  full  rank  if  the  data  conditions  are  changed  (e.g.,  if  the  data 
are  contaminated  by  noise)  [11]. 

4.2.1  Difference  Between  Modeling  Using  Triple  Cor¬ 
relations  and  Using  Autocorrelations 

In  the  previous  section  we  developed  algorithms  which  estimate  the  AR  co¬ 
efficients  using  third  order  cumulants.  Alternatively,  there  are  algorithms 
which  estimate  the  AR  coefficients  of  the  same  model  (4.1)  using  second  or¬ 
der  cumulants  only  [8].  These  two  approaches  represent  different  modeling 
philosophies,  and  below  we  discuss  the  differences  between  the  two. 

If  the  goal  of  the  modeling  procedure  is  to  obtain  parametric  bispectral 
estimates  then  only  the  third  order  cumulants  should  be  used  to  estimate 
the  AR  coefficients.  The  reason  for  not  including  second  order  cumulants  is 
to  avoid  any  responses  that  are  not  bicorrelated  even  when  the  wavenum¬ 
ber  of  one  component  is  the  sum  of  wavenumbers  of  two  other  components. 
If  only  second  order  cumulants  are  used  to  estimate  the  AR  parameters, 
estimates  of  the  target  impulse  response  will  be  obtained,  but  the  bicor¬ 
relation  terms  in  the  birange  profile  cannot  be  obtained.  If  the  goal  is  to 
model  the  system  including  all  time  components  (bicorrelated  or  not)  then 
one  can  use  both  second  and  third  order  cumulants  to  estimate  the  AR 
parameters.  Therefore,  the  system  model  that  is  derived  using  third  order 
cumulants  exhibits  time  components  that  are  engaged  only  in  bicorrelation, 
whereas  the  system  model  derived  using  second  and  third  order  cumulants 
combined  exhibits  all  time  components  of  the  signal.  In  this  chapter,  we 
are  concerned  with  modeling  the  birange  profile,  therefore  we  will  use  third 
order  cumulants  only. 
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h(t)  has  eight  responses  only,  i.e.,  h(t )  =  { h(ti ),  i  —  1, ...  ,8},  where  t5  = 
03 +  U,  and  ts  =  te  +  h-  Assume  also  that  the  responses  {h(t3),  h(t4),  h(U)} 
are  bicorrelated.  If  H(<v)  is  modeled  using  only  third  order  cumulants  then 
only  the  time  components  h(t 3),  and  h(t5)  are  detected.  Thus,  the 

birange  profile  will  have  a  peak  at  (t3,  £4)  and  (£4,  t3),  and  be  zero  elsewhere. 
However,  if  H(iv)  is  modeled  using  second  order  cumulants  then  all  time 
components  are  detected.  If  the  parametric  bispectrum  is  generated  using 
second  and  third  order  cumulants  then  it  will  include  peaks  at  both  ( t6 ,  t7) 
and  at  (t3,ii). 

Notice  that,  while  second  order  cumulants  cannot  be  used  to  model  the 
bispectrum,  third  order  cumulants  may  be  used  to  model  the  spectrum  (Lii 
and  Rosenblatt  [78])  given  that  full  rank  cumulant  slices  exist. 

4.2.2  Experimental  Results 

Blade-Sphere  Example: 

Algorithm  #  1  of  Section  4.2  was  used  to  estimate  the  birange  profile  of  the 
blade  sphere  target  example  of  Chapter  III.  In  this  example,  as  in  Chapter 
III,  1000  data  points  were  used  to  estimate  the  birange  profile.  Specifically, 
from  these  1000  data  points  we  estimated  the  triple  correlation  R(k,k)  for 
k  =  0,...,p  using  (4.6).  We  solved  for  the  AR  coefficients  using  (4.7). 
Finally,  the  birange  profile  was  obtained  using  (4.4).  Figure  4.1  shows 
the  birange  of  the  blade  sphere  target  when  estimated  using  an  AIt  filter 
with  complex  coefficients  of  orders  10  and  12.  This  figure  shows  that  the 
birange,  as  estimated  using  the  above  algorithm  is  sensitive  to  model  order 
as  indicated  by  the  presence  of  spurious  peaks  (denoted  by  a,  b,  and  c)  in 
the  birange  when  the  model  is  increased  from  10  to  12.  The  estimates  of  the 
coordinates  of  the  interaction  (see  Table  4.1),  however,  are  more  accurate 
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when  the  birange  is  estimated  using  the  above  algorithm  than  using  Fourier 
transform  or  even  using  Gaussian  function  fitting  on  the  Fourier  transform 
(compare  Table  4.1  with  Table  3.1).  For  this  example,  using  more  than  a 
single  slice  (as  in  Algorithm  #  2)  to  estimate  the  AR  parameters  has  an 
insignificant  effect  on  the  birange  estimate  of  this  target. 

Figure  4.2  shows  both  the  impulse  response  and  the  birange  profile 
of  the  blade  sphere  target  using  AR  model  of  order  12.  The  purpose  of 
Figure  4.2  is  to  show  the  correspondence  between  the  impulse  response  and 
the  birange  of  the  blade-sphere  target  when  both  are  obtained  using  AR 
modeling.  Note  that  all  the  responses  of  the  target  are  modeled  in  the 
impulse  response. 

Scaled  Model  Target  Examples: 

The  triple  correlation  matrix  given  in  Algorithm  #  2  was  used  (in  a  least 
squares  sense)  to  estimate  the  AR  parameters  and  the  birange  profile  for 
each  of  the  five  model  aircraft  discussed  in  Section  3.6.5.  The  data  were 
decimated  so  that  an  unambiguous  window  of  6  nanoseconds  is  obtained, 
and  the  triple  correlations  used  to  obtain  the  AR  parameters  were  estimated 
from  the  decimated  data.  The  decimation  rule  followed  is 

Mh)  =  \Wk-i)  +  2  H{fk)  +  H(fM)) 

where  Hd(fk)  is  the  decimated  response  at  frequency  /&.  A  triple  correlation 
lag  of  21  data  points  was  used.  The  data  belong  to  the  1.5  —  12  GHz 
frequency  band  with  a  frequency  increment  of  50  MHz.  Algorithm  #  1  was 
also  used  but  found  to  give  inferior  results,  and  are  not  presented. 

Figures  4.3-4. 5  show  the  birange  profiles  for  three  commercial  aircraft 
obtained  using  AR  models  of  order  14  and  using  the  Fourier  transform:  Fig¬ 
ure  4.6  shows  the  impulse  responses  for  three  aircraft  obtained  using  Fourier 
transform;  and  AR  modeling.  Notice  that  both  the  impulse  responses  and 
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the  birange  profiles  obtained  using  AR  models  do  not  correlate  well  with 
those  obtained  using  Fourier  transform. 

The  difference  in  the  impulse  responses  may  be  due  to  the  fact  that 
AR  parameters  estimated  using  triple  correlations  model  the  bispectrum 
(or  birange)  of  the  target  and  not  its  spectrum  (or  impulse  response).  We 
have  tried  several  experiments  using  different  model  orders  and  using  dif¬ 
ferent  triple  correlation  lags;  we  have  also  tried  incorporating  singular  value 
decomposition  into  the  AR  parameter  estimation  process.  The  results  in 
4. 3-4. 5  represent  the  best  results  we  could  obtain  using  this  algorithm.  Fig¬ 
ure  4.5  represents  the  most  compatible  and  Figure  4.4  the  least  compatible 
fit  between  AR  and  Fourier  birange  profiles  of  the  five  targets.  It  is  also 
clear  from  Figures  4. 1-4.5  that  the  AR  modeling  algorithm  given  in  this 
chapter  can  reproduce  the  Fourier-based  birange  profiles  using  synthesized 
data  but  has  difficulty  resolving  closely  spaced  birange  peaks  using  real 
radar  data. 

Figures  4.7  shows  how  the  impulse  response,  that  models  the  birange 
profile  of  the  DC-10  aircraft,  changes  as  a  function  of  model  order.  Notice 
that,  as  the  model  order  is  increased,  spurious  peaks  are  obtained  as  shown 
in  Figure  4.7.  However,  a  spurious  peak  at  time  t(  in  the  impulse  response 
does  not  always  produce  a  spurious  peak  in  the  birange  profile  unless  there 
are  two  other  peaks  (not  necessarily  spurious)  at  tj,  U  —  tj,  and  tj  —  U. 
Figure  4.8  shows  the  birange  profile  of  the  DC-10  aircraft  using  the  same 
model  orders  as  in  Figure  4.7.  Notice  that,  unlike  the  impulse  response 
shown  in  Figure  4.7,  increasing  the  model  order  does  not  always  introduce 
spurious  birange  peaks  (see  the  upper  left  plot  in  Figure  4.8  for  AR(20)). 

The  parametric  birange  profiles  obtained  using  third  order  cumulants 
assume  that  the  data  sequence  is  random  with  hidden  bi correlations.  If  we 
assume  that  the  data  is  deterministic,  then  it  is  possible  to  model  the  signal 
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using  second  order  cumulants  or  autocorrelations.  In  this  case  the  Fourier- 
based  birange  profile  can  be  obtained  using  jS(<i,t2)  =  h(<1)/i(t2)h*(<i  +  <2) 
and  the  AR-based  birange  profiles  can  be  obtained  using  Yule- Walker  equa¬ 
tions  [60].  Figures  4.9-4.10  show  the  birange  profile  and  the  impulse  re¬ 
sponse  for  the  aircraft  DC-10  obtained  using  Fourier  transform  and  AR 
modeling  and  using  the  same  data  and  the  same  decimation  procedure  as 
in  the  above  experiments.  These  Figures  show  that  when  the  data  is  as¬ 
sumed  to  be  deterministic,  the  birange  profiles  obtained  using  the  Fourier 
transform  can  be  matched  with  higher  resolution  fidelity  by  AR  based  esti¬ 
mates  using  second  order  cumulants.  The  reason  for  the  agreement  between 
the  Fourier-based  birange  profile  and  the  AR  birange  profile  is  that  both 
are  directly  estimated  from  from  the  impulse  responses  using  (4.3)  and 
there  is  an  agreement  between  the  impulse  response  obtained  using  Fourier 
transform  and  using  the  AR  parameters.  This  experiment  suggests  that 
the  difficulty  in  modeling  the  birange  may  be  related  to  the  fact  that  we 
are  using  the  triple  correlation  of  complex-valued  data  to  estimate  the  AR 
parameters. 

4.3  Discussion 

The  results  obtained  on  modeling  scaled  model  aircraft  using  the  AR-based 
algorithm  proposed  in  Section  4.2  are  not  encouraging  in  the  sense  that  no 
closely  spaced  peaks  in  the  birange  are  resolved  and  little  compatibility 
between  Fourier-based  birange  profiles  and  AR-based  birange  profiles  is 
observed.  This  lead  us  to  investigate  certain  modeling  aspects  related  to 
the  type  of  scattering  data  that  we  are  trying  to  model. 

One  of  the  key  issues  of  the  AR  modeling  procedure  adopted  in  this 
Chapter  is  that  the  radar  scattering  data  are  complex- valued.  Unlike  mod¬ 
eling  the  spectral  responses  of  complex-valued  signals,  there  seems  to  be 
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a  significant  difference  between  modeling  real-valued  and  complex-valued 
data.  This  difference  is  explained  by  the  following  theorems  where  com¬ 
plex  exponential  signals  with  implicit  bicorrelations  are  used  as  examples  of 
complex-valued  data.  The  following  theorems  show  that  parameters  of  an 
AR  model  of  a  certain  order  can  be  found  such  that  they  uniquely  satisfy 
equation  (4.8)  and  yet  they  cannot  be  used  to  model  the  true  bispectrum 
of  the  complex- valued  signal.  Thus,  these  theorems  suggest  other  possible 
reasons  for  discouraging  AR  modeling  results. 

4.3.1  On  AR  Modeling  of  Complex- Valued  Signals 

The  estimation  of  the  bispectrum  using  AR  modeling  requires  the  choice 
of  the  model  order.  The  model  order  is  dependent  on  the  number  of  peaks 
in  the  bispectral  response.  If  real-valued  data  is  used  and  the  birange  has 
I  peaks  in  one  symmetry  sector,  then  the  number  of  peaks  in  the  impulse 
response  that  models  the  birange  is  upper  bounded  by  3 1,  which  corre¬ 
sponds  to  a  model  order  of  67  (assuming  real  AR  coefficients).  This  result 
has  been  shown  in  (42)  when  modeling  real-valued  sinusoids  with  implicit 
bicorrelations. 

Modeling  complex- valued  signals  requires  a  different  approach  from  that 
used  to  model  real-valued  signals.  The  following  theorems  explain  the  mod¬ 
eling  aspect  of  complex-valued  signals  (complex  exponentials  used  as  an 
example  of  complex-valued  signals)  and  the  model  order  needed.  These 
theorems  show  a  fundamental  difference  between  modeling  real- valued  sig¬ 
nals  and  complex-valued  signals. 

Theorem  4.1:  Let  H(w)  be  of  the  following  form 

i 

#(w)  =  £(exP  OH™ +  M) 

i=i 

+  exp  {j(wtbi  +  0W)}  +  exp  {j(u(tai  +  tbi)  +  0ai  +  0fci)}fl4.1O) 
where  0ai,Obi  are  uniformly  distributed  random  variables  over  [0, 27r] .  Note 
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that  H(u>)  includes  I  implicit  couplings  as  indicated  by  the  third  component 
of  (4-10).  There  is  a  unique  set  of  complex  AR  coefficients  {an}n=o>  with 
ao  =  1,  such  that 

i 

ani?j/(n  —  A;,n  — /)  =  0.  (4.11) 

n=0 

However,  when  this  set  of  AR  coefficients  is  used  in  (4-%)  and  (4-3),  the 
result  is  in  general  not  the  true  birange  profile. 

Proof:  The  triple  correlation  of  H (u>)  is  given  as  (replacing  u>i ,  and  u;2  by 
k  and  l  respectively. 

/ 

Rir{k,  l )  =  £[exp  +  ltbi)}  +  exp  {j{ktbi  +  ltai)}}.  (4.12) 

i=i 

and  the  birange  profile  is 
/ 

BH(U,h)  =  £[*(<i  -  tai,t 2  -  hi)  +  S(h  -  tbi,t2  -  tai) )  (4.13) 

~  i=i 

Then, 

/ 

^2anRH(n~  k,n-  l)  = 

71=0 

I  / 

£  £k»{exp  {j{tai{n  -k)  +  tbi(n  -  l))} 

n=0 i=l 

+  exp  {j{tbi{n  -  k)  +  tai(n  -  /))}}]  (4.14) 

/ 

=  £[(exp  {-j(taik  +  tbil)} 


+  exp  {-j{tail  +  hk)})  £  on  exp  {j(tai  +  <bf)n}] 

n=0 


This  expression  is  zero  if  and  only  if  (see  [60])  the  AR  parameters  {a;} 
satisfy 


i  i 

£  anzn  =  J](z  -  exp  {j(tai  +  4,-)})  (4.15) 

n=l  i=l 

Therefore,  it  is  sufficient  to  choose  a  model  of  order  I  to  model  I  responses  in 
the  bispectrum.  The  impulse  response  obtained  using  these  AR  parameters 
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has  I  peaks  at  4t  4-  4i,  t  =  1,. . . and  is  given  as  (using  equation  (4.2)) 


^  exp  {,/7r(l  -  2t/T)}  -  exp  {j7r(l  -  2(tai  +  4.))/^} 
Then  the  birange  is  obtained  using  (4.3)  as 


(4.16) 


B{ti,  4)  (jC  exp  _  24/T)}  -  exp  {j7r(l  -  2(4;  4-  4i))/r} 


exP  OM1  -  2 4/T)}  -  exp  {j'tt(1  -  2(4;  +  4;))/r} 


\£i  exp  {— jtt(1  -  2(4  4-  t2)/T)}  -  exp  {— i?r(l  -  2(4;  +  tbi))/T}J 
The  birange  shows  a  peak  at  (4,4)  when  either  4,4,  or  4  +  4  is  equal 
to  4;  4-  4;.  For  example,  if  /i(4)  =  a,  A(4)  =  /?,  (where  ct,/3  ^  0,  for 
an  all-pole  model)  and  A(4  4*  4)  =  7  where  limf;.+£jb_4£(i;+t(>i.{7}  — >  00,  and 
(4  4-  4  =  4t  +  4;),t  =  1, . . . then  the  birange  at  (4,4)  is  given  as  such 


(4-17) 


Similarly,  (4.17)  holds  V4,4  9  4  +  4  =  4*  4-  4i,  4  ~  4  =  4i  4-  4«,  °r 
4  -  4  =  4»  +  hi-  Note  that  the  true  birange  peak  at  (4*,  4f)  belongs  to 
the  straight  line  defined  by  {4  4-  4  =  4*  4-  4;},  however  the  magnitude  of 
this  peak  is  comparable  to  the  other  peaks  along  the  same  line  and  cannot 
be  distinguished  from  the  others.  The  only  birange  profile  that  can  be 
modeled  correctly  corresponds  to  the  case  where  4;  =  4«  =  0,  since  in  this 
case  lim(,l(j-o{-#(4,4)}  — >  (oo)(oo)(oo). 

Q.E.D. 

Theorem  4.2:  Let  H(w)  be  of  the  following  form 
1 

H{w)  =  £leXP{i(W*a.-  +  0ai)} 

i-1 

4-  exp  {j(wtbi  +  0bi)}  4-  exp  {j(u;(4;  4-  hi)  4-  9ai  4-  0ti)}p-l8) 

where  9ai,9bi  are  uniform, ly  distributed  random  variables  over  (0, 27r) .  Then, 
there  is  a  unique  set  of  complex  AR  coefficients  {an}^i0,  with  a0  =  1,  such 
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that 


(4.19) 


where 


zi 

Y  anRcH(n  —  k,n  —  l)  —  0. 

n=0 


RcH(k,l)  =  E{H*(n)H(n  +  k)H(n  +  l) 

+H(n)H*(n  +  k)H{n  +  l)  +  H(n)H{n  +  k)H'(n  +  ($.20) 

These  parameters  uniquely  model  the  true  birange  profile. 

Proof:  Using  (4.20),  the  third  order  cumulant  of  the  signal  in  (4.18)  is 
zi 

RcH{k,  l )  =  ]OexP  {j(khi  +  Itbi)}  +  exp  {j(ktbi  +  ltai)} 

i=i 

+  exp  {j{k(tai  +  hi )  +  Ihi)}  +  exP  {j(khi  4-  l{hi  +  hi)} 

+  exp  {j{k(tai  +  hi)  +  Ihi)}  +  exp  {j(ktai  +  l(tai  +  4i))B-21) 

The  proof  then  follows  similar  to  that  in  [42]  for  real  sinusoidal  signals  with 
I  couplings.  The  impulse  response  is  then  given  as 

~  S^xp  iM1  ~  2t/T)}  “  exp  0’tt(1  -  2 {tai  +  hi)/T}  + 

_ Sil _ + 

exp  {jir(l  -  2 t/T)}  -  exp  {j7r(l  -  2(tai)/T} 

_ Cjz _ l 

exp{i7r(l  -  2 t/T)}  -  exp{j7r(l  -  2(ibi)/T}> 

where  c;7-  ^  0,  j  —  1,2,3.  The  birange  is  then  obtained  using  (4.3).  Using 

this  model,  all  responses  at  tai,  hi,  and  tai  +  hi,  i  =  1,  • . .  ,1  are  modeled 

in  the  impulse  response  and  each  set  of  {hi, hi, hi  +  4i}  in  the  impulse 

response  produces  a  single  peak  in  the  birange  such  as 

t,  til1?-  (°o)(co)(oo)  (4.22) 

C1  >c2*“*ca»i  lbt 

This  peak  can  be  easily  distinguished  from  other  peaks  that  are  obtained 
because  of  the  transformation  of  (4.3).  Q.E.D. 

We  have  shown  in  Theorem  4.1  (using  AR  modeling  of  complex  expo¬ 
nentials)  that  although  it  is  possible  to  obtain  AR  parameters  that  satisfy 
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(4.8),  the  derived  system  does  not  necessarily  model  the  the  true  bispec¬ 
trum.  Theorem  4.2  presents  an  alternative  approach  to  modeling  complex 
exponentials.  The  parameters  obtained  as  described  in  Theorem  4.2  model 
uniquely  the  true  bispectrum  of  complex  exponentials.  Theorem  4.2,  how¬ 
ever,  cannot  be  used  to  model  all  types  of  complex-valued  signals,  including 
radar  signals,  but  suggests  that  modeling  complex-valued  time  series,  using 
third  order  cumulants,  may  differ  significantly  from  modeling  real  time  se¬ 
ries.  The  reason  that  we  cannot  use  (4.20)  to  model  radar  signals  is  twofold: 
1)  the  electric  phase  -0  modeled  in  (3.15)  is  not  uniformly  distributed  over 
[0,27r]  (however,  we  may  assume  so  for  mathematical  reasons),  and  2)  the 
model  of  (3.15)  corresponds  to  point  target  scatterers  and  may  not  be  ap¬ 
propriate  to  model  real  aircraft  targets.  The  reason  that  the  birange  of  the 
blade-sphere  target  model  was  recovered  using  third  order  cumulants-based 
AR  modeling  is  because  scattering  from  the  model  is  deterministic  with  a 
response  at  zero  time. 

Another  aspect  of  modeling  radar  target  responses  is  related  to  the 
estimation  procedure  adopted  in  this  study.  Modeling  the  bicorrelated  re¬ 
sponses  in  the  impulse  response  and  then  transforming  the  obtained  impulse 
response  into  a  two-dimensional  profile  (the  birange  profile)  produces  ridges 
that  are  formed  by  straight  lines  in  the  birange  profile.  These  straight  lines 
do  not  belong  to  the  birange  profile  but  are  the  consequence  of  transfor¬ 
mation  from  an  all-pole  impulse  response  into  the  birange  profile.  From 
the  outcomes  of  the  experimental  study  and  using  such  an  argument,  we 
conclude  that  modeling  the  birange  profiles  of  radar  scattering  data  may 
benefit  by  using  a  two-dimensional  AR  model  that  models  third  order  cu¬ 
mulants  directly  instead  of  the  one-dimensional  modeling  procedure  (also 
based  on  third  order  cumulants)  adopted  in  this  study. 
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4.4  Summary 


The  purpose  of  this  chapter  is  to  devise  and  test  an  algorithm  that  can 
be  used  to  obtain  high  resolution  birange  profiles  by  modeling  the  mea¬ 
sured  frequency  response  using  AR  models.  It  is  shown  that  the  AR  model 
needed  to  estimate  the  birange  may  differ  from  the  AR  model  needed  to 
estimate  the  impulse  response.  Algorithms  based  on  third  order  cumulants 
of  the  measured  data  were  used  to  obtain  experimental  birange  profiles 
for  different  model  orders.  The  AR.  birange  profiles  were  then  compared 
to  the  Fourier  transform  based  birange  estimates.  The  results  obtained 
are  not  encouraging  as  indicated  by  the  incompat ability  of  Fourier-based 
birange  profiles  and  the  AR-based  birange  profiles.  This  lead  us  to  fur¬ 
ther  investigate  the  properties  of  (triple  correlation)-based  AR  modeling  of 
complex- valued  radar  signals.  The  results  of  this  investigation  indicate  that 
modeling  complex-valued  exponential  signals  using  third  order  cumulants 
does  not  result  in  the  true  birange  profile  and  a  slightly  different  model¬ 
ing  strategy  is  required.  Even  when  the  birange  is  modeled  correctly,  the 
estimated  birange  profile  from  the  impulse  response  exhibits  false  birange 
peaks  that  belong  to  straight  lines  in  the  two-dimensional  birange  plane. 
Such  a  disadvantage  may  be  avoided  by  using  a  two-dimensional  AR  model 
that  is  based  directly  on  the  third  order  cumulants.  Henceforth,  the  model¬ 
ing  aspect  of  the  two-dimensional  birange  profiles  that  are  estimated  from 
complex-valued  signals  remains  an  open  research  area. 
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Figure  4.1:  Birange  profile  for  the  blade-sphere  target  using  AR  modeling 
based  on  third  order  cumulants  for  models  of  order  10  and  12  (for  compari¬ 
son,  the  box  in  the  lower  left  shows  the  birange  profile  generated  by  Fourier 
transform,  and  the  box  in  the  upper  right  corner  shows  the  true  location 
of  birange  responses). 
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Figure  4.2:  Impulse  Response  and  Birange  profile  for  the  blade-sphere  tar 
get  using  an  AR  model  of  order  12  (based  on  third  order  cumulants). 
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Figure  4.3:  Birange  profile  for  DC-10  using  Fourier  transform  and  AR 
model  of  order  14  (based  on  third  order  cumulants). 
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Figure  4.4:  Birange  profile  for  the  Concord  aircraft  using  Fourier  transform 
and  AR  model  of  order  14  (based  on  third  order  cumulants). 


IMPULSE  RESPONSE 


Figure  4.6:  Impulse  response  for 
transform  and  AR  model  of  order 
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Figure  4.7:  Impulse  response  for  DC-10  aircraft  using  AR  models  (based 
on  third  order  cumulants)  of  different  orders. 
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Figure  4.8:  Birange  profiles  for  DC-10  aircraft  using  AR  models  (based 
third  order  cumulanls)  of  different  orders. 
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Table  4.1:  Comparison  between  AR  and  Fourier  birange  features  of  the 
synthesized  target  example. 


Three-point  target  Sr  —  .015  m 

Peaks 

AR 

Fourier 

True 

n 

.1752 

.1715 

.1789 

r2 

.3161 

.3185 

.3111 
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Chapter  5 


CLASSIFICATION  OF 
RADAR  SIGNATURES 


5.1  Introduction 

The  focus  in  this  chapter  is  on  using  the  birange  profile  as  a  basis  for  dis¬ 
criminating  between  target  signatures  and  for  classifying  unknown  airborne 
radar  targets.  A  radar  target  identification  system  is  designed  to  distin¬ 
guish  birange  profiles  of  unknown  objects  as  belonging  to  one  of  a  set  of  M 
classes,  with  each  class  corresponding  to  a  particular  target.  The  goal  is  to 
devise  and  test  efficient  algorithms  that  can  be  employed  for  classification 
purposes. 

Bispectral  processing  of  radar  signals  has  an  inherent  advantage  over 
spectral  processing;  namely  the  suppression  of  additive  Gaussian  noise  or 
any  additive  unwanted  random  signal  with  symmetric  probability  density 
function.  This  advantage  motivates  our  interest  in  using  the  birange  profile 
for  target  recognition  purposes,  particularly  if  the  backscatter  waveforms 
are  corrupted  with  disturbances  that  can  be  modeled  as  zero  mean  additive 
Gaussian  noise  [22].  Therefore  it  is  of  interest  to  assess  the  reliability  of 
target  identification  when  using  the  birange  profile  as  a  feature  source. 

Although  total  suppression  of  Gaussian  noise  is  possible  when  the  scat¬ 
tering  measurements  are  recorded  over  an  infinite  frequency  range,  in  prac- 
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tice  the  radar  measurement  frequency  range  is  limited.  Therefore,  some 
noise  components  remain  in  the  birange  profile.  Consequently,  the  percent¬ 
age  of  correct  classification  of  radar  targets,  using  the  birange  as  a  feature 
source,  is  expected  to  decrease  as  the  additive  noise  power  is  increased 
or  the  signal-to-noise  ratio  (SNR)  is  decreased.  To  analytically  study  the 
effect  of  the  noise  on  classification  performance,  we  need  to  know  the  un¬ 
derlying  distribution  of  the  birange  as  a  function  of  noise  power  cr2.  This 
distribution  is  derived  in  Section  5.2. 

The  process  of  target  classification  is  approached  in  two  ways;  paramet¬ 
ric  and  nonparametric.  The  parametric  approach  requires  the  knowledge  of 
the  underlying  distribution  of  the  birange  profile  and  the  a  priori  probabil¬ 
ity  of  occurence  of  each  target.  Such  an  approach  is  useful  for  estimating 
the  optimal  classification  performance  that  could  be  achieved  using  the 
birange  as  a  feature.  Due  to  the  lack  of  knowledge  of  the  a  priori  proba¬ 
bilities  of  possible  targets,  and  due  to  the  computational  burden  imposed 
by  parametric  classification,  it  is  often  more  appropriate  to  use  simpler 
nonparametric  classifiers. 

Using  the  birange  or  other  time  domain  profile  in  target  classification 
may  be  rewarding  when  nonparametric  recognition  methods  are  required. 
If  the  underlying  distribution  of  the  scattered  waveform  is  not  known  or  if 
the  target  probability  of  occurence  is  not  known  a  priori,  then  one  might 
consider  classification  using  processed  radar  signatures  such  as  the  birange 
profile  or  the  transient  response.  Target  classification  under  these  con¬ 
ditions  may  lack  optimality.  Therefore,  if  parametric  classification  is  not 
possible  due  to  certain  conditions  imposed  on  the  data  then  both  frequency- 
domain  and  time-domain  features  should  be  considered.  Hence,  the  birange 
profile  can  be  considered  as  a  feature  source  for  nonparametric  classifiers. 
The  focus  in  this  chapter  is  on  devising  and  testing  the  performance  of 
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some  non-parametric  classifiers.  Such  classifiers  include  cross- correlation 
and  nearest  neighbor  methods. 

5.2  The  Underlying  Distribution  of  a  Single 
Response  in  the  Birange 

The  underlying  distribution  of  the  birange  profile  of  noisy  radar  signals  is 
an  indication  of  the  effect  of  additive  Gaussian  noise  on  the  estimate  of  the 
birange.  The  derivation  of  the  likelihood  function  of  birange  estimates  is 
the  first  step  in  parametric  classification  studies.  The  likelihood  function  of 
the  estimate  of  the  birange  profile  at  a  point  given  that  this  birange 

profile  corresponds  to  a  target  of  class  i  is  denoted  as 

p(B(kyl)/i)  Vfc,Z  €{!,..., N}  (5.1) 

where  N2  is  the  number  of  pixels  in  the  birange  profile,  where  about  .75 N2 
pixels  are  used  to  display  the  birange  (the  birange  hexagon  is  displayed 
inside  the  square  defined  by  N2  pixels).  Therefore,  the  likelihood  function 
of  the  birange  profile  is  simply  the  joint  likelihood  function  for  all  pixels. 

Assume  that  the  estimate  of  the  birange  is  obtained  using  the  direct 
approach  without  segmenting  the  data.  Then, 

B,(rk,r,)  =  KhTMh^,)  (5.2) 

where  h'ik  =  h{k  +  n,*  is  the  real-valued  noisy  impulse  response  of  target 
i  evaluated  using  a  discrete  Fourier  transform  of  the  band-limited  target 
frequency  response  and  K  is  a  constant.  If  the  impulse  response  h'ik  is  com¬ 
puted  without  zero-padding  the  measured  frequency  response,  then  it  can 
be  shown  [75]  that  the  transient  response  at  range  rk  has  a  Gaussian  distri¬ 
bution  with  mean  h{k  and  variance  a2.  Note  that  are  independent 

identically  distributed  (i.i.d.)  random  variables.  Therefore,  the  likelihood 


function  of  the  pixel  (rfc,rj)  of  the  birange  profile  is  just  the  likelihood 
function  of  the  product  of  three  independent  identically  distributed  (iid) 
samples  of  the  transient  response  (as  given  in  (5.2)  ).  This  can  be  done  in 
two  ways;  using  a  transformation  of  random  variables,  or  using  the  Mellin 
transform.  Unfortunately,  in  either  case  it  is  difficult  to  derive  a  closed 
form  expression  of  the  underlying  distribution.  These  two  approaches  are 
investigated  below. 


A-Transformation  of  random  variables: 


It  is  possible,  using  the  following  transformation  of  the  random  variables 
{hik}  (where  /i,-*  is  the  impulse  response  of  target  i  at  range  r &)  to  a  new  set 
of  random  variables  {zjt},  to  derive  the  likelihood  function  of  the  product 
in  (5.2)  as  the  marginal  distribution  of  the  new  random  variable  23  = 
hikhuhi(k+l)- 


hik  =  (5.3) 

hu  =  z2 
hi(k+l)huhik  —  23 


The  Jacobian  for  the  transformation  from  hik  to  Zk  is  |J|  =  Then, 
the  likelihood  function  of  the  birange  at  (r^,r{)  is  simply  the  marginal 
distribution  of  23  given  as: 


O°o  1 


exp 


222i 

■(*.  -  hiky  -  (z,  -  h„y  hi 


(5.4) 


( k+l ) 


2cr2 


dz\dz2 


where  hik,hn,h^k+i)  are  samples  from  the  noise  free  transient  response  at 
ranges 

rktriiTk  +  ri  for  a  target  of  class  i.  It  is  possible  to  solve  numerically  for 
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B-Using  the  Mellin  transform 

The  distribution  for  the  product  of  n  independent  nonstandardized  random 
variables  can  be  found  as  an  infinite  summation  using  the  Mellin  transform 
and  the  complex  residue  theorem  [19] .  For  n  =  3,  the  likelihood  function 
of  a  single  pixel  in  the  birange  is 


p(B(k,l)/i)  =  cpq 

p= 0 q=Q 

k=  1/2  9=0 

where  ±  corresponds  to  >  0,  and  where 


B2' 

p 

8<r6 

8<r«J 

1 2-9 


£  <5-5> 


’  B2  ’ 

k 

[-1,-1 

8a\ 

8<r» 

1  2 -q 


B 

4<r3 


At  = 


exp 


(2t r)3/2cr3 


At  =  2^4  iW 

<t3  j= [ 

The  coefficients  Cjq ,  <4,  can  be  obtained  using 

q  /  \  /  3  \ 

dk,  = 


(q-m) 


uP 


b  - 


r(s  +  b  +  p ) 


*-( _ 

ds *  \(s  +  b  —  l)(s  +  b) . . .  (s  +  b  —  2  +  p); 

a  = -p,  p  =  0,1,2,... 

5  *’  *  2’2’2’"' 


(5.6) 

(5.7) 


(5.8) 

(5.9) 

(5.10) 

(5.11) 


*•«  = 

(5.12) 

,  , ,  /2s  +  1  3  h?-  \ 

=  n  2  '1*2^; 

(5.13) 

,,  .  "  a  z  a(a  + 1)  z2 

$(a;c;,)  =  >  +  cl,  +  c(c  +  x)  2,  +  - 

(5.14) 
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where 


where  /i,j  is  the  noise  free  transient  response  at  range  rj,  1  <  j  <  N  for 
target  i.  Also,  the  ( q  —  m)th  derivative  of  the  function  $  can  be  evaluated 

using  Leibnitz’s  rule  as 
r  n  i 


n$°w 


(5.15) 


q-m  Ji  Jn-1 


=  £  E  £  (,7m)  - 

jt=0j3=0  jn=0  \  /  \  Jn  J 

where  jc,  c  =  1,2  must  satisfy  j1  +  j2  +  . . .  +  jn  —  9  —  m. 

The  derivative  of  the  function  can  be  obtained  using  the  following 


•  •  •  P j  $ir_il)(i) . . . 


where 


y6(0(a)  =  s¥l\s  +  b  +  p)  +  Y,(-i)l+1 

i=0 

x(f  4-  l)!(s  +  b  —  1  +  z)-*-1 


(5.16) 


(5.17) 


and  where  ^(^(s  +  b  +  p)  is  the  digamma  function,  and  $W(s  -f  b  +  p)  is 
the  polygamma  function  given  as  for  l  >  1 


$(0(s  +  b  +  p)  =  (~l)l+H\  £(/+!  +  a)-(,+b+p) 


(5.18) 


The  above  is  the  likelihood  function  for  the  birange  profile  at  a  point 
(M),  i  <  k,l,<  N.  In  order  to  obtain  the  likelihood  function  for  all  points 
k,l  we  need  to  compute  the  joint  distribution  for  all  pixels  in  the  birange 
domain.  The  birange  pixels,  however,  are  correlated  (along  vertical  and 
horizontal  lines)  i.e., 


#  £{Si(Jfe>J)}jB{3i(lfe>»)} 


(5.19) 

(5.20) 


where  Brfk.l )  denotes  the  birange  of  target  z  at  (rfc,r/).  Also  notice  that 
the  variance  of  a  single  response  in  the  birange  is 


var (Bi{k,  l ))  =  o’6  +  cr4  [h]k  +  h\  +  h\k+l) 


(5.21) 
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+  a2  [h2kh2i2  +  h2{lh2( fc+{)  +  h2kh2^k+ ,)] 

This  indicates  that  the  variance  is  proportional  to  the  third  power  of  the 
noise  variance  and  is  also  dependent  on  the  data. 

Figure  5.1  shows  an  example  of  the  likelihood  function  (or  conditional 
probability  density  function)  of  one  pixel  in  the  birange  profile  of  noisy  radar 
signals.  In  this  figure,  the  values  hy  for  the  noise  free  impulse  response  of 
target  i  at  ranges  rj,  j  =  1,2,3  are  given  as  fin  =  1,  ha  =  1,  ha  =  .4. 

Notice,  that  the  distribution  is  not  symmetric,  and  its  skewness  depends 
on  the  noise  free  values  of  the  transient  response  at  ri,r2,  and  ri  +  tv  The 
distribution  in  Figure  5.1  corresponds  to  cr  =  1;  therefore,  cr4  =  cr6  =  1.  If 
<r  >  1  then  the  variance  of  the  above  distribution  will  increase  rapidly  with 
a. 

Using  the  above  likelihood  function,  it  is  possible  to  redefine  the  bi¬ 
range  as  a  deterministic  component,  determined  by  the  expected  value  of 
B ,  added  to  a  random  variable  with  a  probability  density  function  as  given 
in  (5.4).  This  definition  of  the  birange  may  be  helpful  in  experimental  tar¬ 
get  recognition  studies  where  it  is  possible  to  obtain  an  estimate  of  the  noisy 
birange  by  adding  a  randomly  generated  number  whose  probability  density 
function  is  as  shown  in  Figure  5.1  to  the  triple  product  of  transient  re¬ 
sponse  samples.  It  is  also  possible,  by  assuming  that  the  birange  pixels  are 
uncorrelated,  to  obtain  an  approximate  evaluation  of  a  parametric  target 
classification  algorithm  based  on  the  birange.  Furthermore,  a  more  accu¬ 
rate  expression  of  the  probability  of  detection  using  the  bispectrum  can  be 
obtained  by  using  the  above  distribution  instead  of  the  normal  (asymptotic 
approximation)  distribution  used  in  [56,57}. 
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likelihood  function 


Likelihood  function  for  one  bispectra!  sample 


*4-3-2*10  12  34 

bispectrum 


Figure  5.1:  Likelihood  function  of  a  single  response  in  the  birange 
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5.3  Classification  Techniques 


The  classification  of  radar  signals  using  the  birange  profile  as  a  feature 
source  is  developed  in  this  section.  The  focus  is  on  nonparametric  classifica¬ 
tion  methods.  Although  it  is  possible  to  design  a  parametric  classifier  such 
as  Bayesian  classifier  using  the  likelihood  function  derived  in  the  previous 
section,  it  is  very  computationally  inefficient  to  do  so.  Further,  there  seems 
to  be  no  significant  reasons  for  classifying  targets  parametrically  using  pro¬ 
cessed  radar  data  since  optimal  parametric  classification  performance  can 
be  achieved  using  the  measured  data  directly. 

Three  non-parametric  birange-based  classification  methods  are  discussed 
in  this  section.  These  classifiers  are:  the  birange  nearest  neighbor  classifier, 
and  the  birange  cross-correlation  classifier.  Employing  any  of  these  classifi¬ 
cation  methods  is  dependent  on  the  available  information  about  the  target, 
and  in  particular  on  the  zero-time  reference. 

5.3.1  The  Birange  Nearest-Neighbor  Classifier 

The  distribution  of  birange  iV2-tuples  in  the  feature  space  indicates  that  the 
classes  are  distributed  on  convolved  non-planer  surfaces  which  cannot  be 
separated  by  hyperplanes  except  possibly  in  very  localized  regions.  There¬ 
fore,  linear  discriminants  cannot  be  used  to  separate  classes  of  radar  targets 
in  terms  of  birange  profiles.  The  nearest  neighbor  (NN)  decision  rule  is  a 
less  constrained  classification  method  that  can  be  easily  employed  to  dis¬ 
criminate  between  targets  based  on  the  birange  profile.  The  NN  rule  has 
performed  satisfactorily  when  used  for  classification  of  frequency- domain 
target  signatures  [26,37]. 

The  nearest  neighbor  decision  rule  can  be  used  to  classify  unknown 
radar  targets  if  the  zero-time  reference  of  the  target  with  respect  to  the 
radar  is  known  a  priori.  Thus,  the  NN  rule  is  used  when  absolute  phase 
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information  is  available  (zero-time  reference  is  known).  The  NN  rule  in  this 
case  is  given  as, 

Assign  the  unknown  target  u  to  class  ( i )  if 
di  =  min  { dj },  j  =  1, . . . ,  M  where 

j 

di  =  ,/£EW  jf  (5-22) 

V  tl  t2 

where  is  the  birange  profile  of  target  x.  Notice  that  this  decision 

rule  is  based  on  selecting  the  class  whose  Euclidean  distance  to  the  unknown 
target  prototype  is  the  smallest.  Therefore,  it  is  necessary  for  the  birange 
profile  of  both,  the  unknown  target  and  the  catalog  target  to  be  displayed 
exactly  in  the  same  position  in  the  (ti,t2)  plane. 

The  lack  of  information  about  the  target  zero-time  reference  suggests 
the  use  of  a  classification  algorithm  that  takes  into  consideration  all  possible 
positions  of  the  target  in  the  time  domain.  An  alternative  to  using  NN  rule 
as  discussed  above  is  to  include  birange  profiles  of  all  possible  positions  of 
the  target  in  the  catalog  of  each  class.  The  decision  rule,  then,  is  based  on 
the  nearest  neighbor  among  all  targets  and  among  all  possible  positions  of 
each  target.  Thus,  this  algorithm  assigns  the  unknown  target  to  class  i  if 

di  =  min  {nun  {djr}}  j  =  l,...,M,  0  <  r  <  (5.23) 

where  A /  is  the  frequency  increment  of  the  measured  scattering  data.  The 
birange  of  each  target  catalog  must  be  known  at  all  possible  choices  of 
r.  This  classification  algorithm  is  computationally  inefficient  because  it 
requires  the  computation  of  the  birange  profile  of  the  catalog  (or  the  un¬ 
known)  for  each  possible  r. 

When  the  target  zero-time  reference  is  unknown,  an  alternative  ap¬ 
proach  to  NN  classification  is  to  apply  a  linear  shift  r  to  the  transient 
response.  This  can  be  done  either  by  shifting  the  transient  response  of  the 


145 


unknown  target,  or  by  applying  a  ramp  shift  (as  a  function  of  frequency) 
with  slope  t  to  the  phase  of  the  measured  frequency  response  H(f)  as 

H'{f)  =  H(f)exv{2nfT}  Vr  <  (5.24) 

Which  corresponds  to  shifting  the  birange  profile  along  the  ti  =  t2  axis.  In 
fact,  applying  a  linear  phase  shift  to  the  data  corresponds  to  uni-directional 
correlation  between  the  two-dimensional  birange  profile  and  the  catalog. 
This  form  of  cross-correlation  considers  all  possible  positions  the  zero-time 
reference  along  the  radar-target  axis.  Computationally,  however,  it  is  ineffi¬ 
cient  to  compute  the  birange  profile  of  the  unknown  target  and  the  catalog 
targets  and  then  correlate  in  the  time  domain. 

5.3.2  Cross-correlation  Classifier 

The  classification  of  birange  profiles  using  cross-correlation  is  developed  in 
this  section.  This  classifier  assumes  that  the  target  zero-phase  reference  is 
known  a  priori.  The  focus  in  this  section  is  on  the  definition  of  the  cross¬ 
correlation  classifier.  The  goal  of  this  algorithm  is  to  find  a  catalog  element 
(i)  that  matches  the  bispectral  response  of  an  unknown  target  Bu(rl)r2). 
That  is  one  wishes  to  minimize 

min{  /  [  (Bi{rur2)  -  Bu(rur2))2  drxdr2}  =  (5.25) 

*  Jr\  Jr 2 

min {J  J  B?(r1,r2)dr1dr2  +  J  J  5^(r1,r2)dr1d7’2  (5.26) 

-/  J  Bi(ri,r2)Bu(ri,r2)dridr2} 

If  the  first  two  terms  in  (5.29)  are  fixed,  this  entails  maximizing 

/  /  Bi(r\ ,  r2)Bu(r\ ,  r2)dr\dr2  (5.27) 

Jr\  Jr 2 

and  since  the  target  zero-phase  reference  is  assumed  to  be  known,  then  this 
is  equivalent  to  maximizing 

riu(0,0)=  (5.28) 
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_ Jr.  Sn  Bi{r1,r2)Bu{rx,r2)  dn  dr2 _ 

[Jr.  Sn  \Bu{rur2)\2drxdr2y  [/P1  /rj  |.B;(r1,r2)|2dr1dr2] 2 

where  I\u(0,0)  is  the  normalized  cross  correlation  of  the  catalog  target 
birange  profiles  for  i  =  1,2, ...,M  and  the  test  target  response  u.  Using 
Fourier  transform  identities  and  Parseval’s  theorem,  we  find  that  the  cross¬ 
correlation  can  be  written  for  the  discrete  frequency  case  as 

riu(o,o) _ ™mHfuh)KUum _ T  (529) 

[Ea  Ea  \Hfuh)\'Y  [Ea  Ea  W/..A)I2]’ 

The  test  target  is  classified  as  target  c  if 

rcu(0,0)  =  max{riu(0,0)}  i  =  1,. . .  ,M.  (5.30) 

5.4  Experimental  Study 

A  comparison  between  the  performance  of  the  cross-correlation  classifier 
using  the  birange  and  the  performance  of  other  optimal  and  suboptimal 
classifiers  is  given  below.  The  comparison  includes  classification  in  additive 
white  Gaussian  noise,  additive  colored  Gaussian  noise,  and  additive  non- 
Gaussian  noise.  Classification  with  azimuth  ambiguity  is  also  investigated. 
Also,  the  effect  of  extraneous  uncatalogued  scatterers  on  the  classification 
performance  is  examined. 

The  probabilities  of  target  misclassification  at  different  signal-to-noise 
ratios  are  estimated  using  Monte-Carlo  simulations.  When  the  a  priori 
probability  of  a  target  hypothesis  is  not  known,  then  N  test  samples  are 
drawn  randomly  and  then  used  to  determine  whether  the  classifier  gives 
the  correct  decisions  for  these  samples  or  not.  If  the  a  priori  probabilities 
Pi,  i  =  1, . . . ,  M  ( M  is  the  number  of  hypotheses)  are  known  then  n,  =  P{N 
samples  are  drawn  from  class  (i),  to  which  noise  is  added  and  then  these 
samples  are  tested  by  the  classifier. 
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The  classifiers  simulated  in  this  study  are  the  following: 

1-  Nearest  Neighbor  rule  based  on  the  backscatter  data  without  employ¬ 
ing  any  signal  processing.  Given  that  the  measured  backscatter  is  Hu  = 
[Hu( u>i), . . . ,  Hu(uk)}  (where  K  is  the  number  of  frequencies  used)  then  the 
NN  rule  classify  the  unknown  as  target  i  if 


=  j  =  1 . M.  (5.31) 


2-  Cross-correlation  based  on  its  one-dimensional  time-domain  response. 
For  this  algorithm,  target  {i}  is  chosen  such  that 


r[u  =  max 


£n  |fc,(n)|» 


) j  —  !)•••) M 


(5.32) 


Note  that  this  is  equivalent  to  maximizing  the  down-range  cross-correlation 
between  the  unknown  target  impulse  response  and  the  catalog  impulse  re¬ 
sponse. 

3-  Cross-correlation  of  birange  profiles  as  discussed  in  Section  5.3.3. 

The  performance  of  the  cross-correlation  classifier  using  the  birange  is 
dependent  on  the  birange  estimation  procedure  (or  the  estimation  of  triple 
correlation).  The  amount  of  segmentation  used  in  computing  the  triple  cor¬ 
relation  had  a  significant  effect  on  the  performance  of  the  classifier.  Over¬ 
segmenting  the  data  may  result  in  a  degraded  classification  performance. 
The  triple  correlation  lag  used  and  the  number  of  data  points  also  influence 
the  performance  of  the  classifier.  Finally,  removing  the  average  from  both 
the  unknown  and  the  catalog  improves  the  classifier  performance. 

A  Monte-Carlo  simulation  was  performed  to  evaluate  the  performance 
of  the  classifiers  mentioned  above.  The  data  base  used  in  these  examples 
has  been  frequently  used  in  radar  target  identification  studies  [6,7,8].  The 
data  base  consists  of  experimental  measurements  in  the  frequency  band 
from  1  —  12  GHz  of  scale  models  of  commercial  aircraft.  A  photo  with 


148 


corresponding  names  of  these  aircraft  is  shown  in  Figure  3.9.  The  scaled 
data  corresponds  to  measurements  of  the  radar  cross  section  (RCS)  of  full 
scale  aircraft  in  the  HF/VHF  frequency  band  (8-58  MHz).  Details  of  these 
measurements  can  be  found  in  [64]. 

Decision  statistics  for  each  target  are  computed  at  a  fixed  noise  level, 
and  total  statistics  of  classification  error  for  all  targets  are  obtained.  One 
hundred  experiments  were  performed  for  each  test  target  (for  a  total  of  500 
experiments).  The  entire  test  is  repeated  at  different  noise  levels.  Finally 
the  misclassification  (error)  percentage  is  plotted  versus  signal  to  noise  ra¬ 
tio.  When  classifying  with  the  birange,  the  best  performance  was  obtained 
when  the  data  is  segmented  into  K  =  5  records  of  M  =  21  samples  each 
with  a  correlation  lag  L  =  10  (see  Section  2.3.4). 

Figure  5.2  shows  the  classification  performance  for  five  commercial  air¬ 
craft  with  complete  azimuth  information  using  additive  white  Gaussian 
noise.  The  catalog  consists  of  scattering  data  for  five  commercial  aircraft 
at  0°,  10°,  and  20°  azimuth.  The  performance  of  the  nearest  neighbor  (NN) 
algorithm  is  optimal  for  this  case.  The  birange  classifier  is  outperformed 
by  the  impulse  response  classifier.  We  also  re-iterate  that  the  performance 
of  the  birange  classifier  is  very  sensitive  to  the  triple  correlation  estimation 
procedure.  The  performance  of  the  birange  classifier  may  degrade  if  the 
number  of  correlation  lags  is  altered  or  if  the  window  function  is  changed 
(the  window  used  in  this  study  is  given  in  (2.60)).  This  performance  figure 
shows  that  birange  profile  estimates  can  be  used  effectively  in  radar  tar¬ 
get  identification.  Increasing  the  number  of  data  samples  and  employing 
an  optimized  classification  scheme  may  improve  the  performance  of  the  bi¬ 
range  classifier  and  reduce  its  sensitivity  to  the  triple  correlation  estimation 
procedure. 

Figures  5.3  and  5.4  show  the  classification  performance  when  additive 
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colored  noise  generated  by  passing  white  Gaussian  noise  through  an  AR  or 
MA  filter  respectively.  The  AR  filter  coefficients  {a»},-=0  are 
[.5, .6, .7, .8, .7, .6, .5, 0,0, .5, .6, .7, .8, .7, .6,. 5]  (see  [61]  for  the  autocorrela¬ 
tion  of  the  noise  generated  using  this  filter),  and  the  MA  filter  has  coeffi¬ 
cients  [1 ,  .8] .  The  target  azimuth  is  assumed  to  be  completely  known.  The 
NN  rule  (which  may  be  suboptimal  in  this  case)  applied  to  the  frequency- 
domain  data  outperforms  the  time-domain  classifiers.  Also,  the  perfor¬ 
mance  of  the  birange  classifier  compares  favorably  with  the  performance 
of  the  impulse  response  classifier.  The  degradation  of  the  time-domain 
classifiers  compared  to  the  nearest  neighbor  rule  is  slightly  less  for  the  ad¬ 
ditive  white  noise  case,  which  may  suggest  a  comparable  performance  of  all 
classifiers  under  other  colored  noise  conditions. 

Figure  5.5  shows  the  classification  performance  when  additive  non- 
Gaussian  noise  is  used  (the  square  root  of  a  Weibull  distributed  random 
variable  added  to  both  the  in-phase  and  quadrature  components  of  the 
data).  The  azimuth  is  assumed  to  be  completely  known.  The  performance 
of  the  birange  classifier  is  improved  and  is  comparable  to  the  other  classi¬ 
fiers,  which  may  indicate  a  significant  role  for  the  birange  in  classification 
of  unknown  targets  in  a  non-Gaussian  noise  environment. 

Figures  5.6  and  5.7  show  the  classification  performance  when  the  az¬ 
imuth  is  known  only  to  within  ±10°  (Figure  5.6)  and  ±20°  (Figure  5.7). 
Although  the  nearest  neighbor  rule  is  not  optimal  in  this  case,  it  outper¬ 
forms  the  time-domain  classifiers.  Further,  the  performance  of  the  birange 
classifier  degrades  compared  to  the  impulse  response  classifier  when  the 
azimuth  ambiguity  range  increases.  In  fact,  if  the  target  is  assumed  to  be 
known  within  ±30°  (not  shown  in  the  Figures)  then  the  performance  of  the 
birange  classifier  degrades  significantly  compared  to  the  impulse  response 
classifier.  These  figures  show  that  the  birange  is  sensitive  to  changes  in 
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target  azimuth  position.  This  sensitivity  may  be  explained  by  the  fact  that 
changing  the  azimuth  may  introduce  additional  multiple  interactions  and 
delete  others.  Although  these  interactions  also  affect  the  impulse  response, 
they  more  strongly  affect  the  birange. 

Figure  5.8  shows  the  classification  performance  when  the  classifier  is 
mis-informed  about  the  target  azimuth  position  with  an  error  of  ±10°. 
That  is,  a  target  at  azimuth  A°  is  compared  with  the  catalog  targets  at  az¬ 
imuth  A0  ±  10°.  This  type  of  mismatch  in  design  specifications  affects  the 
classification  performance  of  all  classifiers  including  the  birange  classifier. 
It  is  clear  from  Figure  5.8  that  time-domain  classification  techniques  out¬ 
perform  the  nearest  neighbor  classifier  in  this  case.  Although  the  birange 
classifier  is  sensitive  to  change  ’n  target  azimuth,  Figure  5.8  shows  that  the 
birange  is  less  affected  by  inaccurate  azimuth  information  than  the  nearest 
neighbor  classifier. 

5.4.1  Effect  of  Extraneous  Scatterers 

In  the  following,  we  discuss  the  effect  of  adding  extraneous  scatterers  to  the 
unknown  target  frequency  response  (scatterers  not  included  in  the  catalog) 
on  the  performance  of  the  classifiers  under  investigation.  This  type  of  clas¬ 
sification  problem  may  occur  when  clutter  (in  the  form  of  point  scatterers 
with  or  without  interactions)  is  detected,  or  when  the  catalog  target  model 
is  inaccurate  or  incomplete. 

The  effect  of  this  type  of  disturbance  is  simulated  as  a  set  of  additive 
point  scatterers  with  variable  scattering  magnitude  located  at  different  po¬ 
sitions  with  respect  to  the  target.  The  ratio  of  scattering  from  the  target 
plus  extraneous  scatterers  to  scattering  from  the  extraneous  scatterers  only 
is  denoted  by  SER  (Signal-to- Extraneous  signal  Ratio). 

Tables  5.1-5.10  show  classification  results  obtained  for  different  scenar- 
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Nearest  Neighbor 

Impulse  Response  (cross-con.) 

Bsraoge  Profile  (cross-con.) 


Figure  5.2:  Classification  performance  of  five  commercial  aircraft  with 
known  azimuth  and  additive  white  Gaussian  noise 
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Nesrert  Neighbor 

Impel*  Response  (crow-con.) 

Bimge  Profile  (ciow-con.) 


Figure  5.3:  Classification  performance  of  five  commercial  aircraft  with 
known  azimuth  and  additive  colored  Gaussian  noise  generated  by  an  AR 
filter 
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Figure  5.4:  Classification  performance  of  five  commercial  aircraft  with 
known  azimuth  and  additive  colored  Gaussian  noise  generated  by  an  MA 
filter 
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Nearert  Neighbor 

inpolse  Respoote  (crciKorc.) 

linage  Profile  (croii-corr.) 


Figure  5.5:  Classification  performance  of  five  commercial  aircraft  with 
known  azimuth  and  additive  non-Gaussian  noise 
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Signd-To-Noise  Ratio  [dB] 


Figure  5.6:  Glassification  performance  of  five  commercial  aircraft  with 
known  azimuth  within  ±10°  and  additive  white  Gaussian  noise 
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Neueit  Neighbor 

bnpolje  Respooae  (croiworr.) 

Mange  Profile  (croci-con.) 


Figure  5.7:  Classification  performance  of  five  commercial  aircraft  with 
known  azimuth  within  ±20°  and  additive  white  Gaussian  noise 
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-f-  Neirest  Neighbor 
_  IupttUe  Reroute  (croir-corr.) 

_  Binge  Profile  (crou-corr.) 


ios  of  extraneous  scatterers  added  at  different  positions  with  respect  to  the 
target.  The  error  rates  shown  in  these  tables  represent  the  percentage  of 
misclassification  when  scattering  from  extraneous  point  scatterers  (these 
scatterers  are  described  in  the  captions)  is  added  to  the  frequency  response 
of  all  five  targets.  Also,  scattering  from  unknown  targets  is  assumed  noise- 
free  since  the  purpose  of  these  experiments  is  to  evaluate  the  performance 
of  the  birange  classifier  when  backscattered  signals  include  uncatalogued 
responses.  The  extraneous  scatterers  used  to  generate  the  results  in  Ta¬ 
bles  5.7  and  5.9  include  an  interaction  whose  response  appears  at  —.8  ns, 
and  those  used  to  generate  the  results  in  Tables  5.8  and  5.10  include  an 
interaction  whose  response  appears  at  0  ns.  The  bicoherences  for  both 
interactions  are  assumed  fixed  with  a  value  of  40%. 

The  following  conclusions  can  be  drawn  from  those  experimental  results: 
First,  classification  with  nearest  neighbor  rule  is  the  most  sensitive  to  the 
presence  of  extraneous  uncatalogued  scatterers.  Second,  classification  with 
impulse  response  is  less  sensitive  to  extraneous  scatterers  than  the  NN 
classifier  but  unaffected  by  changing  the  position  of  extraneous  scatterers. 
Finally,  we  see  that  the  birange  is  less  sensitive  to  the  presence  of  extraneous 
scatterers  than  the  nearest  neighbor  classifiers.  The  birange  classifier  is 
even  less  sensitive  to  uncatalogued  scatterers  if  the  responses  from  such 
extraneous  scatterers  do  not  coincide  with  the  response  of  the  unknown 
target.  This  insensitivity  may  be  due  to  the  fact  that  while  the  presence  of 
extraneous  scatterers  produces  additional  peaks  in  the  impulse  response,  it 
takes  three  extraneous  scatterers  with  interaction  to  produce  a  peak  in  the 
birange  profile  of  the  unknown  target. 

Figures  5.9-5.10  show  a  comparison  between  classification  results  (plot¬ 
ted  as  probability  of  misclassification  versus  the  SER).  It  is  clear  from  these 
figures  that  classification  using  the  birange  profiles  outperforms  classifica- 
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tion  using  the  radar  cross  section  measurements  (using  nearest  neighbor) 
even  when  the  azimuth  of  the  unknown  target  is  partially  known  (or  known 
within  a  certain  azimuth  range).  The  extraneous  signatures  used  in  these 
figures  are  three  point  scatterers  where  the  location  of  the  response  due 
to  each  scatterer  is  random  and  uniformly  distributed  over  [-f  ,f],  where 
T  =  A_,  (where  A /  is  the  frequency  increment  of  the  measured  data). 
The  locations  of  the  responses  of  the  extraneous  scatterers  are  independent 
identically  distributed.  Two  hundred  experiments  including  five  targets  (a 
total  of  1000  experiments)  are  simulated  for  each  SER  ratio. 

Figure  5.11  shows  classification  results  obtained  when  the  frequency  re¬ 
sponses  of  nine  extraneous  scatterers  is  added  to  the  frequency  responses 
of  the  five  aircraft.  The  location  of  the  response  due  to  each  scatterer  is 
random  and  uniformly  distributed  over  [— f  ,f).  Figures  5.9  and  5.11  show 
that  classification  using  the  impulse  response  outperforms  classification  us¬ 
ing  the  birange  when  the  number  of  extraneous  scatterers  is  increased  to 
nine.  This  result  can  be  explained  since  by  adding  the  responses  of  nine 
extraneous  scatterers  to  the  data,  the  likelihood  of  having  a  birange  peak 
due  to  extraneous  scatterers  is  higher  than  that  when  only  three  extraneous 
scatterers  are  present.  Further,  each  of  the  five  targets  has  less  than  nine 
scattering  centers,  which  indicates  that  if  the  number  of  the  extraneous 
scatterers  exceeds  the  number  of  scattering  centers  then  the  impulse  re¬ 
sponse  classifier  outperforms  the  birange  classifier  by  a  significant  margin. 
Figures  5.9  and  5.11  also  show  that  the  performance  of  the  nearest  neigh¬ 
bor  classifier  is  inadequate  whether  the  number  of  extraneous  scatterers 
exceeds  the  number  of  scattering  centers  or  not. 
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5.5  Summary 


This  chapter  is  concerned  with  studying  the  classification  aspect  of  the  bi¬ 
range  profile  when  used  as  a  target  recognition  tool.  The  purpose  of  this 
study  is  to  investigate  the  feasibility  of  employing  the  birange  as  a  feature 
source  in  target  recognition  systems  and  to  examine  the  conditions  un¬ 
der  which  the  birange  classifier  performs  satisfactorily.  The  classification 
schemes  suggested  in  this  chapter  are  suboptimal  in  the  sense  of  minimizing 
classification  error.  However,  using  these  classification  schemes,  it  is  possi¬ 
ble  to  identify  a  framework  of  application  for  the  birange  classifier  in  terms 
of  noise  conditions  and  availability  of  a  priori  information  about  aspect 

I 

angle.  The  comparisons  made  in  this  study  include  the  performance  of  the 
nearest  neighbor  classifier,  the  impulse  response  cross- correlation  classifier, 
the  birange  cross-correlation  classifier.  The  conditions  investigated  include 
additive  white  Gaussian  noise,  additive  colored  Gaussian  noise,  additive 
non-Gaussian  noise,  additive  uncatalogued  scatterers,  and  error  in  aspect 
information. 

The  classification  results  in  this  chapter  can  be  summarized  as  follows: 
The  nearest  neighbor  classifier  outperforms  both  the  impulse  response  and 
the  birange  classifiers  except  for  the  case  where  the  unknown  target  is  mis¬ 
represented  in  the  catalog.  An  unknown  target  is  misrepresented  in  the 
catalog  if  the  response  due  to  extraneous  scatterers  is  added  to  its  fre¬ 
quency  response  or  if  its  aspect  angle  is  not  included  in  the  catalog.  The 
impulse  response  classifier  outperforms  the  birange  classifier  except  for  the 
case  of  non-Gaussian  additive  noise.  When  scattering  from  extraneous  un¬ 
catalogued  scatterers  is  added  to  the  frequency  response  of  the  unknown 
target  then  classification  with  the  birange  outperforms  classification  with 
impulse  response  if  the  signal-to-extraneous  scatterers  ratio  is  low  and  if 
the  number  of  extraneous  scatterers  does  not  exceed  the  number  of  scatter- 
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ing  centers  along  the  target.  Finally,  by  considering  the  outcomes  of  this 
study,  we  conclude  that  classification  with  the  impulse  response  is  more  re¬ 
liable  than  other  classifiers  (all  circumstances  combined).  However,  it  may 
be  rewarding  to  use  a  joint  decision  rule  that  combines  both  the  impulse 
response  and  the  birange  classifiers. 
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Table  5.1:  Classification  error  rate  when  a  single  extraneous  point  target  is 
placed  at  -2  ns  (within  the  nose  area). 


SER  [dB] 

NN 

TIME 

BIRANGE 

1 

0 

0 

0 

-4 

20 

0 

20 

-9 

40 

20 

40 

-12 

40 

20 

40 

-15 

60 

40 

40 

Table  5.2:  Classification  error  rate  when  a  single  extraneous  point  target  is 
placed  at  0  ns  (within  the  wing  area). 


SER  [dB] 

NN 

TIME 

BIRANGE 

1 

0 

0 

0 

-4 

0 

0 

0 

-9 

0 

20 

20 

-12 

m 

60 

20 

-15 

K9 

60 

40 

Table  5.3:  Classification  error  rate  when  a  single  extraneous  point  target  is 
placed  at  1.2  ns  (within  the  tail  area). 


SER  [dB] 

NN 

TIME 

BIRANGE 

1 

0 

0 

0 

-4 

0 

0 

0 

-9 

■a 

40 

0 

-12 

n 

40 

0 

-15 

60 

40 

20 

I 
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Figure  5.9:  Comparison  between  classification  using  measured  data  (us¬ 
ing  nearest  neighbor),  impulse  response,  and  birange  profile  as  the  sig- 
nal-to-extraneous  scatterers  ratio  (SER)  changes  with  complete  azimuth 
information  (three  extraneous  scatterers  used). 
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Signd-To-Extroneous  Scatter ers  [dB] 


Figure  5.10:  Comparison  between  classification  using  measured  data  (us¬ 
ing  nearest  neighbor),  impulse  response,  and  birange  profile  as  the  sig- 
nal-to-extraneous  scatterers  ratio  (SER)  changes  with  partial  azimuth  in- 
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Figure  5.11:  Comparison  between  classification  using  measured  data  (us¬ 
ing  nearest  neighbor),  impulse  response,  and  birange  profile  as  the  sig- 
nal-to-extraneous  scatterers  ratio  (SER)  changes  with  complete  azimuth 
information  (nine  extraneous  scatterers  used). 
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Table  5.4:  Classification  error  rate  when  two  extraneous  point  targets  are 
placed  at  -2  ns  and  0  ns  respectively. 


SER  [dB] 

NN 

TIME 

BIRANGE 

-2 

0 

0 

0 

-8 

0 

20 

20 

-12 

El 

40 

40 

-14 

El 

60 

40 

-16 

60 

40 

60 

Table  5.5:  Classification  error  rate  when  two  extraneous  point  targets  are 
placed  at  -2  ns  and  1.2  ns  respectively. 


SER  (dB] 

NN 

TIME 

BIRANGE 

-2 

0 

0 

0 

-8 

■a 

mm 

20 

-12 

Kb 

20 

-14 

n 

40 

40 

-16 

80 

40 

40 

Table  5.6:  Classification  error  rate  when  two  extraneous  point  targets  are 
placed  at  0  ns  and  1.2  ns  respectively. 


SER  [dB] 

NN 

TIME 

BIRANGE 

-2 

0 

0 

0 

-8 

0 

mm 

20 

-12 

El 

BE 

40 

-14 

n 

80 

40 

-16 

40 

80 

40 

Table  5.7:  Classification  error  rate  when  three  extraneous  point  targets  are 
placed  at  -2  ns,  0  ns,  and  1.2  ns  respectively. 


SER  [dB] 

NN 

TIME 

BIRANGE 

-4 

0 

0 

0 

El 

20 

20 

-13 

19 

20 

20 

-16 

80 

40 

60 

-18 

80 

60 

60 
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Table  5.8:  Classification  error  rate  when  three  extraneous  point  targets  are 
placed  at  -3  ns,  2  ns,  and  3  ns  respectively. 


SER  [dBj 

NN 

TIME 

BIRANGE 

0 

0 

0 

I 

El 

20 

mm 

Mr- 

m 

40 

msM 

1 

40 

60 

40 

-18 

60 

60 

40 

Table  5.9:  Classification  error  rate  when  three  extraneous  point  targets 
are  placed  at  -2  ns,  0  ns,  and  1.2  ns  respectively  with  partial  azimuth 
information  (within  20°). 


SER  [dB] 

NN 

TIME 

BIRANGE 

-4 

0 

0 

0 

-10 

20 

30 

20 

-13 

40 

60 

40 

-16 

60 

60 

60 

-18 

60 

60 

80 

Table  5.10:  Classification  error  rate  when  three  extraneous  point  targets  are 
placed  at  -3  ns,  2  ns,  and  3  ns  respectively  with  partial  azimuth  information 
(within  20°). 


SER  [dB] 

NN 

TIME 

BIRANGE 

-4 

0 

10 

0 

msm 

H 

20 

0 

-13 

n 

30 

20 

-16 

60 

40 

40 

-18 

60 

40 

40 
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Chapter  6 

CONCLUSIONS 

6.1  Summary  of  Results 

This  report  has  investigated  the  use  of  bispectral  processing  methods  for 
radar  backscatter  processing  and  for  radar  target  identification.  Bispectral 
processing  methods  were  adapted  to  the  radar  signature  analysis  problem, 
and  resulted  in  the  so-called  birange  profile.  Analysis  of  the  birange  pro¬ 
file  showed  that  it  was  able  to  detect  interactive  scattering  terms  in  the 
backseat tered  signal.  This  result  led  us  to  investigate  whether  this  addi¬ 
tional  information  could  be  effectively  used  as  a  discriminant  for  target 
identification. 

The  state  of  the  art  in  time-domain  radar  signal  processing  was  briefly 
described  in  Chapter  II.  Research  on  time-domain  signature  analysis  is 
focused  on  the  extraction  of  the  major  scattering  components  from  the 
impulse  response  of  a  radar  target.  High  resolution  spectral  estimation 
techniques  are  also  being  employed  to  realize  this  goal.  The  focus  in  these 
studies  is  on  extracting  important  scattering  characteristics  from  the  im¬ 
pulse  response.  The  present  consideration  of  the  bispectrum  serves  as  a 
tool  that  enhances  the  feature  extraction  process. 

In  Chapter  III  we  studied  the  birange  profile  as  a  tool  for  detecting  mul¬ 
tiple  interactions.  A  definition  for  the  bispectrum  of  complex  radar  signals 
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was  introduced.  Also,  a  birange  estimation  technique  based  on  the  “indi¬ 
rect  bispectral  estimation  method”  with  slight  modifications  and  suitable 
for  radar  signals  was  proposed  and  used  in  this  report.  An  interpretation  of 
the  birange  based  on  a  specified  scattering  model  which  is  basically  a  sum 
of  complex  exponentials  with  implicit  couplings,  is  presented.  Examples  on 
the  birange  of  hypothetical  targets  were  also  given.  The  concept  of  birange 
processing  was  tested  using  experimental  measurements  of  scattering  from 
a  tilted  flat  plate  and  a  generic  aircraft  model  with  removable  parts,  then 
using  real  aircraft  models  such  as  Boeing  707,  727,  747,  DC  10,  and  the 
Concord.  Finally,  a  feature  extraction  algorithm  based  on  Gaussian  fitting 
was  used  to  locate  the  interactions  in  the  birange  profiles. 

The  extracted  features  were  compared  with  the  physical  structure  of  the 
corresponding  target  model  and  they  support  the  claim  that  the  bispectrum 
may  be  used  to  detect  multiple  interactions.  It  is  shown  that  it  is  possible 
to  detect  many  responses  in  the  birange  profile  using  scattering  data  at  the 
resonance  and  the  optical  scattering  regions.  Also,  it  is  possible  to  detect 
responses  in  the  birange  that  are  undetectable  in  the  transient  response. 
However,  not  all  interactions  can  be  recovered  from  the  bispectrum.  It  was 
also  shown  that  by  using  the  “bicoherence”,  one  may  obtain  an  estimate 
of  the  magnitude  of  interactions;  thus  the  bicoherence  is  useful  to  separate 
interactions  from  specular  responses  in  the  impulse  response. 

A  birange  profile  estimation  algorithm  based  on  AR  modeling  of  radar 
signatures  was  introduced  in  Chapter  IV.  This  algorithm  was  tested  using 
both  synthesized  scattering  data  and  real  data.  The  Fourier-based  birange 
profile  of  the  synthesized  target  was  reproduced  using  this  algorithm,  but 
the  Fourier-based  birange  profiles  of  real  data  were  slightly  incompatible 
with  those  obtained  using  this  algorithm.  The  reasons  for  incompatability 
of  the  AR-based  birange  estimates  with  the  Fourier-based  estimates  were 
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investigated.  It  was  shown  that  AR  modeling  of  complex- valued  signals  may 
produce  birange  estimates  that  are  entirely  different  from  those  obtained 
using  Fourier  transform.  Furthermore,  the  transformation  involved  in  the 
birange  estimation  procedure  (from  a  one-dimensional  all-pole  response  to  a 
two-dimensional  profile),  produces  ridges  that  do  not  belong  to  the  birange 
of  the  target  in  the  birange  plane.  Therefore,  we  conclude  that  parametric 
modeling  of  the  birange  may  benefit  from  a  two-dimensional  AR  modeling 
approach. 

Chapter  V  considered  the  effectiveness  of  birange  profiles  for  target 
classification.  Tests  were  conducted  using  non-Gaussian  noise,  azimuth 
ambiguity,  and  extraneous  scatterers.  In  general,  classification  results  using 
the  birange  profile  only  were  inferior  to  those  using  the  impulse  response 
or  nearest  neighbor  classifiers.  There  were  some  cases  where  the  birange 
classifier  outperformed  other  methods.  In  the  first  case  the  additive  noise 
was  non-Gaussian  and  in  the  second  case  the  responses  of  three  extraneous 
point  scatterers  were  added  to  the  measured  target  frequency  response. 

In  general,  performance  gains  using  the  birange  profile  were  slight.  Also, 
the  computational  needs  of  the  birange  profile  are  much  higher  than  any 
of  the  other  methods  considered.  Based  on  this  argument,  these  studies 
suggest  that  the  birange  profile  can  provide  modest  classification  perfor¬ 
mance  gain  at  the  expense  of  a  large  increase  in  computation.  However, 
these  results  are  based  on  initial  studies  using  suboptimum  classifiers,  and 
may  not  be  true  under  other  circumstances. 

The  birange  as  a  time-domain  profile  is  not  an  alternative  to  the  im¬ 
pulse  response  but  an  augmentative  display  of  target  signatures.  If  the 
additional  knowledge,  provided  by  the  birange,  is  weighted  properly,  then 
it  will  have  oositive  impact  on  the  radar  target  recognition  problem.  One 
may  argue  that,  because  bcth  the  transient  response  and  the  birange  pro- 


171 


file  are  derived  from  the  same  data  source,  then  neither  one  will  enhance 
the  identification  process  beyond  what  is  already  achieved  using  the  data 
prior  to  any  signal  processing.  The  validity  of  this  argument  diminishes 
when  the  parameters  of  the  underlying  probability  density  function  are  not 
known  a  priori.  Therefore,  features  extracted  from  the  birange  profile  will 
potentially  enhance  the  identification  process  under  special  circumstances. 

6.2  Continuation  of  Research  in  This  Area 


The  following  research  objectives  are  likely  candidates  for  continuing  re¬ 
search  in  this  branch  of  radar  signal  processing. 

Notice  that  the  triple  correlation  for  any  limited  data  sequence  has 
a  hexagonal  shape  inside  a  squared  region  ( N  X  N ).  Furthermore,  the 
hexagonal  form  of  the  triple  correlation  is  admissable  in  the  (<i, <2)  plane, 
where  “admissable”  means  that  a  replica  of  the  triple  correlation  can  fill 
the  entire  triple  correlation  plane  without  any  gaps.  Hence,  it  is  possible 
to  compute  the  birange  using  hexagonal  fast  Fourier  transform  methods 
(HFFT).  By  using  the  HFFT  or  hexagonal  discrete  Fourier  transform  of 
the  triple  correlation  to  compute  the  birange,  the  number  of  data  storage 
and  computations  is  reduced  by  about  25%  compared  to  rectangular  Fourier 
transform  computation  techniques. 

The  assumed  scattering  model  in  this  report  does  not  include  disper¬ 
sive  targets.  As  in  the  transient  response,  many  forms  of  dispersion  may 
produce  gradually  decaying  interactions  as  a  function  of  time  in  the  bi¬ 


range  profile.  If  the  returned  signal  decays  as  a  function  of  frequency  or 
loses  energy  at  high  frequencies,  then,  tne  triple  correlation  decays  even 


faster.  Similarly,  if  the  returned  signal  is  gaining  strength  at  high  frequen¬ 


cies,  then  the  triple  correlation  increases  significantly  when  including  high 


frequency  components.  Therefore,  it  is  important  to  adjust  the  lag  of  the 


triple  correlation  estimate  such  that  no  significant  information  is  lost  at 
higher  frequencies. 

It  may  be  interesting  to  investigate  the  birange  profile  of  a  class  of 
canonical  scatterers  to  gain  additional  insight  about  the  nature  of  multiple 
interactions.  The  outcomes  of  such  an  investigation  can  then  be  used  to 
better  understand  scattering  from  real  radar  targets. 

A  common  problem  in  bispectral  processing  is  the  storage  of  large  ar¬ 
rays  of  bispectral  signatures.  This  problem  can  be  avoided  by  employing 
sophisticated  bispectral  feature  extraction  methods.  This  can  be  accom¬ 
plished  by  a  thorough  investigation  of  the  parametric  modeling  aspect  of 
the  birange  profiles.  The  extracted  features  can  then  be  presented  to  a 
structural  pattern  recognition  system  to  identify  radar  targets. 

This  report  has  focused  on  the  extraction  of  second  order  multiple  in¬ 
teractions  using  the  bispectrum.  Higher  order  multiple  interactions  can  be 
detected  using  higher  order  spectra  such  as  trispectrum.  An  investigation  of 
this  nature  may  have  significant  impact  on  target  recognition  studies,  par¬ 
ticularly  because  fourth  order  moments  do  not  have  the  same  shortcomings 
as  third  order  moments. 

Radar  signatures  that  are  used  in  this  study  correspond  to  the  hori¬ 
zontal  transmit  horizontal  receive  polarization.  An  interesting  problem  is 
to  incorporate  other  polarization  into  the  birange  profile.  Parameters  such 
as  ellipticity  and  tilt  of  multiple  interactions  may  be  significant  to  radar 
scattering  analysis  [75]. 


173 


Appendix  A 

Generation  of  Noise 


The  experiments  performed  in  this  study  require  the  simulation  of  additive 
white  Gaussian  noise,  colored  Gaussian  noise,  and  non-Gaussian  noise.  The 
white  Gaussian  noise  is  simulated  by  adding  twelve  uniformly  distributed 
random  variables  x ;  where  x,-  is  uniform  over  [—.5,  .5].  The  distribution 
of  the  total  sum  of  these  random  variables  is  the  twelfth  convolution  of  a 
rectangular  pulse  with  itself  which  approximates  a  Gaussian  function. 

The  colored  Gaussian  noise  is  simulated  as  the  output  of  a  stable  linear 
AR  or  MA  filter  whose  input  is  white  Gaussian  noise.  The  coefficients 
of  the  AR  (or  MA)  filter  define  the  form  of  coloring  applied  to  the  white 
noise  signal  (these  coefficients  define  the  spectra  of  the  generated  colored 
gaussian  noise). 

The  distribution  of  a  clutter  random  variable  C  can  be  modeled  by  a 
Weibull  distribution,  where  the  random  variable  D  =  {(C  —  rj) / a}7  has  the 
exponential  distribution 

Vd{D)  =  exp(-D).  (A.l) 

Then  the  probability  density  function  of  C  is 

pc(C)  =  7 / a{{C  -  tj)/ a}7-1  exp  [~ [{C  -  7?)/a}7]  (-4.2) 

The  standard  Weibull  distribution  is  obtained  for  rj  =  0  and  a  =  1  In  this 
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case 


E{C}  =  1  -  .57722/7  +  .98905/72  (A.3) 

VAR{C}  =  (tt2/672)  «  1.64493/7* 

In  order  to  generate  a  random  variable  that  has  a  Weibull  distribu¬ 
tion  we  need  to  generate  an  exponentially  distributed  random  variable 
D,po{D)  =  exp(— D)  then  the  clutter  C  can  be  computed  from  D  - 
{{C  —  7/)/a }7,  where  77,7,  and  a  are  adjustable  parameters. 

The  method  used  to  generate  an  exponentially  distributed  random  vari¬ 
able  is  the  following  [16].  Let  A';,  i  —  0, 1, . . .  be  a  sequence  of  independent 
random  variables  uniformly  distributed  over  [0,1].  Define  a  random  vari¬ 
able  N  for  N  >  0,  and  N  odd,  such  that  [16] 

X,  <  A0,  £  X,  <  Xo,  •  •  • ,  £  Xj  <  Xo,  £  Xj  >  X0  (AA) 

j=i  j- 1  j=i 

Then  let  T  be  the  number  of  sequences  rejected  before  an  odd  N  satisfies 
the  above  conditions  (T  =  0,1,...)  then  the  random  variable  D  =  T  +  X0 
is  exponentially  distributed  with  a  standard  density  exp(— D)  [16].  The 
random  variable  C  can  then  be  generated  from  D  using  D  =  {(C  —  rj) / a}y . 
Finally  notice  that  the  parameter  7  plays  the  role  of  a  (the  standard  de¬ 
viation)  for  a  Gaussian  distribution;  thus  reducing  7  results  in  increasing 
the  energy  of  the  clutter  signal  and  decreasing  the  signal  to  clutter  ratio 

(SCR). 
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